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MEETING OF 1957 MARCH 8 
Dr W. H. Steavenson, President, in the Chair 


The President announced the death of Henry Norris Russell, an Associate 
of the Society, and paid a tribute to his memory, the Fellows standing. 


The election by the Council of the following Fellows was duly confirmed :— 

George Anderson, 54 Winton Road, Ashburton, Victoria, Australia (pro- 
posed by W. H. Vale) ; 

William Bowen Bonnor, 8 Aigburth Drive, Liverpool 17 (proposed by 
J. Kershaw); 

Michael Harvey Briggs, 9 Meadway, Chadderton, Lancs. (proposed by 
F. W. Land); 

Peter John Grant, Uplands, Burley, Ringwood, Hants. (proposed by R. W. 
Cartmel); 

*Michael Guest, 122 Perrywood Road, Great Barr, Birmingham (proposed 
by P. Moore); 

Kenneth William Harrison, 61 Worthing Avenue, Gosport, Hants. (proposed 
by W. H. Day); 

Michael James Hendrie, 25 The Crossways, Westclitf-on-Sea, Essex (pro- 
posed by R. L. Waterfield); 

James Lawrence, Stonyhurst College Observatory, Lancs. (proposed by 
P. J. Treanor); 

John Vetterlein, g1 Ingrave Road, Brentwood, Essex (proposed by W. H. 
Julian); 

*Roderick Vernon Willstrop, The Observatories, Cambridge (proposed by 
R. O. Redman); and 

Raymond William Henry Wright, University College, Ibadan, Nigeria 
(proposed by M. W. Ovenden). 


The election by the Council of the following Junior Members was duly 
confirmed :— 


Muriel Hodgkinson, 78 Hardy Lane, Chorlton-c-Hardy, Manchester 21 
(proposed by Z. Kopal); and 
Richard Anthony James, Downing College, Cambridge (proposed by 
A, W. Vince). 
* Transferred from Junior Membership. 
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One hundred and twelve presents were announced as having been received 
since the last meeting, including :— 

A. G. Meissner, Astronomischer Hand-Atlas zu Riidiger’s Kenntniss des 
Himmels (presented by Flt Lt F. Whalley); 

Melbourne Astrographic Catalogue, Vol. V (presented by Sydney Observatory 
and the International Astronomical Union); 

Harold Jeffreys, Scientific Inference, 2nd edition (presented by the author) ; 

F. Link, Die Mondfinsternisse (presented by Messrs Geest and Portig); 

J. Dufay, Galactic nebulae and interstellar matter (presented by Messrs 
Hutchinson & Co.); and 

G. C. MeVittie, Distance and Time in Cosmology: the Observational Data 
(presented by the author). 


MEETING OF 1957 APRIL 12 
Dr W. H. Steavenson, President, in the Chair 


The President announced that the Council had elected the following 
Associates of the Society :— 
Guillermo Haro, Director of the National Astrophysical Observatory, 
Tonanzintla, Mexico; 
Milton La Salle Humason, Mount Wilson Observatory, California, U.S.A. ; 
and 
Miss Inge Lehmann, Chief Seismologist, Copenhagen (retired). 


The election by the Council of the following Fellows was duly confirmed :— 


John Brian Alexander, 106 Pevensey Road, Eastbourne, Sussex (proposed 
by R. v. d. R. Woolley); 

Winter Horace Barnett, 6 The Bungalows, Pottery Lane, Whittington 
Moor, Chesterfield, Derbyshire (proposed by R. R. S. Cox); 

Thomas James Edgar Palmer, 124 Kenyon Mansions, Queen’s Club Gardens, 
London, W.14 (proposed by F. M. Holborn); 

Michael Douglas Scott-Scott, Ferres, Swanage, Dorset (proposed by H. 
Bondi); 

William Smith, 60 Hillcrest Road, Southlands, Weymouth, Dorset (pro- 
posed by A. F. Alexander); 

George Eric Burdett Stephenson, 124 Duchy Road, Harrogate, Yorkshire 
(proposed by G. Fielder); and 

Harry Wilde, 363 London Road, Ewell, Surrey (proposed by J. Jackson). 


The election by the Council of the following Junior Members was duly 
confirmed :— : 
Vinod Kumar Gaur, Geophysics Department, Imperial College, London, 
S.W.7 (proposed by J. M. Bruckshaw); and 
Brian Richard May, 3 St George’s Mews, Regent’s Park Road, London, 
N.W.1. (proposed by C. W. Allen). 


THE VARIATION OF IONIZATION ALONG A METEOR TRAIL 
J. S. Greenhow and E. L. Neufeld 


(Communicated by A. C. B. Lovell) 


(Received 1957 March 29) 


Summary 

The ionized trails of faint meteors have been investigated by observation 
of radio echoes at two spaced receiving stations. ‘The mean ionization curve 
is much shorter than predicted by present theory, and in particular the rise to 
maximum electron line density is more rapid than expected. Occasional 
meteors appear almost instantaneously with maximum electron density, and 
in many cases irregularities in the ionization considerably distort the diffraction 
effects during trail formation. Evidence for the fragmentation of meteors is 
also presented. These effects are considered in relation to Jacchia’s photo- 
graphic observations of the light curves of faint visual meteors. 





1. Introduction.—A meteor in its passage through the upper atmosphere is 
subject to impacts with air molecules, and the heat liberated is sufficient to 
vaporize completely the meteoric material. Light and ionization are produced 
by collisions between air molecules and evaporated meteor atoms, which are 
moving with relative velocities of 10~70 km sec~! (Sparrow, 1926; Opik, 1933; 
Hoppe, 1937). ‘The variation of light and ionization along a meteor trail 
has been discussed by Whipple (1943) and Herlofson (1948), and the light 
curves have been investigated photographically by Jacchia (1949). Results for 
bright meteors of 0 to —7 magnitude are in general agreement with theory, 
and by analogy it has been assumed that the ionization curves are also of the same 
form. ‘The ionization curves are of considerable interest in studies of radio 
echoes from meteor trails, particularly in atmospheric density and temperature 
measurements, and in the determination of meteor mass distributions (Kaiser, 
1954; Evans, 1954). 

Recently the photographic observations have been extended to faint visual 
meteors of approximately +3 mag, and serious anomalies in the light curves 
have been discovered (Jacchia, 1955). As the majority of radio echoes are 
obtained from meteors even fainter than this, it is important to determine to 
what extent the ionization curves agree with those derived theoretically. ‘This 
paper describes experiments in which the variation of electron line density along 
a meteor trail has been investigated by observation of the radio echoes at two 
widely-spaced stations. 

2. Technique 

(a) Measurement of electron line density.—'Vhe linear electron density « in the 
section of a meteor trail from which the radio echo is obtained may be determined 
from the echo characteristics, together with a knowledge of the equipment 
parameters. ‘Thus for «~10™ electrons cm™ the initial echo amplitude 
AyxaG7!?G;)?/R®? (Lovell and Clegg, 1948). ‘lhe amplitude decays exponen- 
tially and the decay time to 1/e of the initial amplitude is A*/167*D (Herlofson, 
1948). G, and Gp, are the gains of the transmitting and receiving aerials, and R 
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the range of the echo. A is the radio wave-length and D the diffusion coefficient. 
When #210" cm-!, A,ca!4G,?G,'2/R%?. In this case echo amplitude varies 
only as the fourth root of «, and is a rather inaccurate measure of the electron 
density. On the other hand, echo duration 7’ increases with «, where 
T= — (5) (Greenhow, 1952). ¢, m are the electronic charge and mass, c 
7 c 
the checker light. Thus if Dis known, T gives a measure of the linear electron 
density. tae 
(b) Transmitter and receiver.—In order to investigate the shape of the ionization 
curve it is necessary to obtain radio reflections simultaneously from two or more 
points along a given meteor trail. In the present experiments a single transmitter 
generating 300 pulses/sec with a peak power of 50 kw was used. I he pulse length 
was 30s and wave-length 8-27 m. ‘The sensitivity of the equipment was such 
that most of the echoes used in the analysis were from meteors of +6 to + 8 visual 
magnitude. ‘Two receiving stations were set up, one on the same site as the 
transmitter at Jodrell Bank (O,, Fig. 1), and the other 20 km to the east. 
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Fic. 1.—Geometry of reflection, showing separation of reflecting points along a meteor trail with 
Spaced receiving stations. 





Because of the specular reflecting properties of a meteor trail, echoes are 
received at O, and O, from points S, and S,. If V is the meteor velocity and ¢ 
the time taken for the meteor to travel from S, to S,, then S,S,=Vt. If observa- 
tions are restricted to meteors occurring at approximately go° azimuths, then for 
a trail AB at any inclination to the plane of the paper the difference in height 
between S, and S, is given by Ah~(V2)* cot E/l. The largest permissible value 
of Ah is 5 km, and all values less than this are possible for meteor trails suitably 
orientated. The meteor velocity is estimated approximately from the rate of 
increase in echo amplitude as the meteor passes through the first Fresnel zone, 
or, in the case of very slow meteors, from the diffraction fluctuations (Davies and 
Ellyett, 1949). ¢ is given by the time delay in appearance of the echoes at O, 
and QO,, and E from the height and range of the meteor echo, 
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(c) (d) 
Fic, 2.—Film records. Upper trace, echo at O,. Lower trace, echo at Os. 150 pulses'sec 50 c/s 
pp 1 2 : : 
calibration pulses. 
(a) Echo at O, > echo at Oy. Echoes are from near end of trail and no diffraction fluctuations are 
observed at 5,. 


(+) Echo at O, > echo at O,. Echoes are from near beginning of trail and well marked diffraction 
effects are visible at Oy. 


(c) Example of irregular ionization between S, and S,, producing diffraction fluctuations before 
the lower specular reflecting point. 


d) Example of “ decay” echo with x 10'? cm! at S,, giving place to a “‘ long duration’’ echo 
Z - ~ 1 
with x > 108 cm~! at Sy (dh~2-5 km). 


J. S. Greenhow and E. L. Neufeld, The variation of ionization along a meteor trail. 
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(a) 





(d) 


Fic. 7.—Observed diffraction patterns produced by a meteor which rises to maximum electron 
density in a trail length 800 m (a) the echo at O, corresponds to the theoretical curve of Fig. 6(g). 


(b) the echo at O,. (Reproduced by courtesy of Mr K. R. R. Bowde n.) 


J. S. Greenhow and E. L. Neufeld, The variation of ionization along a meteor trail. 
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A radio link between O, and O, enables the signals from S, and S, to 
be displayed side by side at the same station. A switched pre-amplifier accepts 
alternate echo pulses from S, and S,, and the signals pass through the same 
receiving and amplifying channels to a switched beam oscilloscope (Fig. 2). 
Alternate echo pulses from the Jodrell Bank and distant receiver thus appear 
at 150 c/s on the upper and lower traces respectively. The slow time base of 
0-4 sec duration is triggered by the first pulse of predetermined amplitude to be 
received from either station. Ranges are recorded on a separate cathode ray 
tube. A full description of the recording system, without the beam switching, 
has been given previously (Greenhow, 1954). 

(c) Aerial system.—An array of 4 Yagi aerials inclined at 45° to the horizontal 
is used for transmission, giving a beamwidth of + 25° to the first zero of the main 
beam. Identical horizontally-polarized receiving aerials consisting of a single 
dipole and reflector are used at both stations. The transmitting beam looks 
eastward over the distant receiver. ‘Thus when a meteor trail is in the process 
of formation the echo will appear first on the Jodrell Bank receiver, followed a 
short time later by the echo at the distant station. The first echo to appear is 
always from a point higher up the trail than the later echo (Fig 1). 

(d) Calibration.—We are not directly concerned with the absolute sensitivities 
of the two receiving systems, as we only require to measure the relative values of 
a from the ratio of echo amplitudes at O, and O,._ A method which automatically 
takes account of all losses in the system is to compare echo rates at the two 
receivers. For equal sensitivities the rates at O, and O, should be equal, 
irrespective of trail length and whether echoes are coincident or not. 

A second and more direct system of calibration is to observe pairs of echoes 
with zero time displacements. When this occurs the meteor trail is orientated 
in such a way that S, and S, coincide, and the ratio of echo amplitudes at O, and O, 
gives the precise ratio in sensitivities in the direction of the trail. ‘The mean value 
of amplitude ratio for these echoes, after using the first method of correction, 
is 1°01 +004. ‘Thus two separate methods of normalizing equipment sensitivities 
give results which agree closely. 

(e) Errors in amplitude ratio—Echo amplitude is proportional to 
G,'?G,!?/R®? for given « and A. Differences in G, and R for echoes observed 
at the two stations are negligible, and errors due to differences in receiver polar 
diagrams are predominant. To reduce this effect to a minimum, identical broad- 
beam aerials are used for reception. As the transmitting aerial and timing 
measurements restrict echoes to the central parts of the dipole and reflector 
mainbeam where G,,"” varies very slowly, we are chiefly concerned with differences 
in the vertical polar diagrams caused by the effects of ground reflections. At 
both receiving sites the ground in front of the aerial is level to within +A/8 over 
the first Fresnel zone. Asa further precaution measurements have been restricted 
to echoes between elevations of 25° and 50°, where the rate of change of gain with 
E is small for an aerial at a height of 4/2. 

To test the vertical polar diagrams at the two stations the distributions of 
echoing points in elevation have been measured, and are illustrated in Fig. 3. 
The distributions are similar, showing that there are no appreciable systematic 
differences between the receiving systems. A point to point comparison of the 
polar diagrams is given by an inspection of echo pairs with zero time displacements 
(Section 2(d)). The amplitude ratios individual echo pairs are found to vary 
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by only +30 per cent from unity, this figure including the spread in results due 
to all sources of error. ue 

A rate count shows that we may expect approximately one echo pair in every 
fifty to be a random coincidence of two entirely separate meteors. 97 per cent of 
these coincidences can be eliminated from range and timing measurements at 
the two stations, and echo pairs of this type have therefore been rejected. 
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ELEVATION (DEGREES) 
Fic. 3.—Distribution of echoing points in elevation at the two receiving stations, showing similarity 
of vertical polar diagrams. 
— — — — theoretical distribution. 


3. Measurement of the mean ionization curve 

(a) Coincidence rate.—With a station separation of 20 km it might be expected 
that some echoes would appear only at one receiver, because of the finite trail 
length. Observation of the number of non-coincidences does not give any 
information about the shape of the ionization curve, as equal echo rates always 
result in equal numbers of non-coincident echoes at the two stations. 
Measurement of the coincidence rate does give some information on the lengths of 
meteor trails (Manning, Villard and Peterson, 1951), but these workers have not 
assigned a meaning to these lengths in terms of relative ionization densities. 

(6) Distribution of amplitude ratios.—Information about the shape of the 
ionization curve is obtained from the ratios of echo amplitude at two spaced 
stations. For example, with the theoretical ionization curve of Fig. 4 (a) (Herlofson, 
1948) and a maximum possible height separation of o-7H, the echo amplitude from 
S, can never exceed twice the amplitude from S, (when S, and S, lie on the part 
of the curve before %max). On the other hand it is possible for the echo from 5, 
to be infinitely greater than that from S,, if S, and S, lie below the point of maxi- 
mum electron density. The shape of the ionization curve can therefore be 
tested by observation of the ratios in echo amplitudes at the two stations. 

(i) Selection of echoes.—The smallest echoes recorded are those with amplitudes 
twice the receiver noise level, and the greatest number are observed with this 
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amplitude. In order to determine the shape of the ionization curve with 
reasonable accuracy, therefore, measurements are first restricted to those echoes 
greater than 10 times receiver noise (10 x m) on one or other of the receivers, and 
then further restricted to 20 times noise echoes. ‘The minimum detectable signal 
is assumed to be rather less than 4 x n, to allow for the possibility of the triggering 
mechanism failing to accept a small echo (although in practice the apparatus 
triggers reliably on 2 x m signals). Measurements are restricted to short duration 
“decay” echoes, for which «ox echo amplitude (Section 2(a)). Thus ratios in 
line density can be measured up to 2-8/1 or 5-5/1 for the smallest echoes in the two 
groups. ‘The ratior=A,/A, is determined for each pair of echoes, where A, and 
A, are the echo amplitudes observed from the upper and lower reflecting points 
respectively. When echoes occur at only one station, the ratio can only be 
specified within certain limits depending upon the echo amplitude. This value 
is always >5°5 or <o-18 for the second group, and >2-°8 or <o-36 for the first 
group of echoes. 








ELECTRON LINE 


DENSITY a (a) 
(ac ECHO AMPLITUDE A) 





7 T 7 . 





-4H -3H “2H -H ’ 
. HEIGHT 
es 
DIRECTION OF 
FLIGHT 
Ah Ahmox 
Ricnesmnaienenniennlbe 
O km 4 5 
ELECTRON LINE 
DENSITY (b) 
(ec ECHO AMPLITUDE A) 
-4H -3H -2H HEIGHT -H O H 


Fic. 4.—(a) Theoretical ionization curve 
% 9 Pp ee es 
- - I-- - (Herlofson, 1948). 
~max 4 Pmax 3 Pmax 
Abscissa expressed in terms of scale height H, where p/ pmax =exp [—(h—/max)/ HH). p and h are 
atmospheric pressure and height. pmax and hmax refer to point of maximum ionization. 











(5) Proposed ionization curve for faint meteors. 


In order to simplify the interpretation of the results, it is desirable to consider 
echo pairs with a fixed separation in height, rather than the full range of Ah from 
oto 5km. ‘The mean height separation of all echoes with 2>r>o0°5 is found 
to be 3:1 km, while for echoes with ratios 4>r>2 and 0-5 >r>o-25 it is 4°1 km. 
In order to maintain a constant value of Ah for all ratios, therefore, echo pairs 
with close spacings have been rejected. ‘The median value of AA for near unity 


ro7 


values of r then rises to 4.0 km. Thus we may assume a value Ah~o-6°()}/H. 
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The deviations include 75 per cent of the echoes. Although the separation of 
specular reflecting points for meteors observed at only one station cannot be 
determined, the small number of echo pairs of sufficient amplitude to enable 
very high ratios to be measured also show separations of approximately 4 km. 
Thus we may determine the ratio of electron line densities between pairs of points 
separated in height by approximately 4 km at various points along the ionization 
curve. 

(ii) Experimental results.—Distributions of the amplitude ratio r for the two 
groups of echoes are given in Fig. 5(a) and (6). The histograms are markedly 
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Fic. 5.—Distribution of ratios in echo amplitudes for pairs of reflecting points spaced by a height 
difference of 4 km. 
(a) Observed distribution (echoes> 5-6 x minimum detectable level). 
(6) Observed distribution (echoes>2°8 x minimum detectable level). 
(c) Theoretical distribution for ionization curve of Fig. 4 (b). 
(d) Theoretical distribution for ionization curve of Fig. 4(a). 
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asymmetrical as expected, but an unexpected feature is the large number of echo 
pairs in which the signal from S, is greater than 4 or 5 times that from S,. Ratios 
as low as o-1 are observed, indicating increases in electron line density by a 
factor of 10 over Ah~4 km, on the part of the ionization curve before &max. 
On this part of the theoretical curve the amplitude increases with decreasing 
height in the manner exp(h/H), and the maximum possible increase in A for 
Ah~4 km is only a factor of 1-8. It appears, therefore, that the rate of increase 
of ionization is much greater than expected, values up to five times the theoretical 
rate being observed. 

On the other hand, the large number of examples in which the echo from S, 
considerably exceeds that from S, is in agreement with theory. ‘Thus, casual 
inspection of observed ratio distributions suggests that ionized meteor trails 
are shorter than has previously been assumed. In the following paragraph an 
attempt is made to deduce qualitatively the shape of the observed ionization curve, 
using the measured amplitude ratios of Fig. 5. 

(c) Theoretical ratio distributions.--The theoretical distribution of amplitude 
ratios for a meteor trail of given «max, observed in a part of the polar diagram of 
sensitivity S=G,'*G,!"/R**, is simply obtained by sliding a pair of points 
spaced by 0-6 H along the horizontal axis in Fig. 4(a). The histogram showing 
the frequency of occurrence of different values of A,/A, is plotted and the process 
is stopped when A, and A, fall below an arbitrary minimum detectable level. 
The complete distribution is then obtained by superimposing a large number 
of such histograms for all values of S and «max, weighting each curve in relation 
to the collecting area of given sensitivity and the distribution of meteor masses. 
Sporadic meteors have been used in the present experiments, and the well- 
established inverse square law distribution of masses is assumed (Kaiser, 1953). 
The variation of G,"?G,"?/R*? is considered by constructing sensitivity contours 
in the 95 km plane above the earth in the manner discussed by Clegg (1948). 
In practice S and &max were varied in eight steps of 1-4, 98 per cent of all echoes 
being included in this sensitivity range of 64 to 1. A random distribution of 
radiants is assumed. The theoretical ratio distribution is illustrated in Fig. 5 (d), 
and this may be compared with the experimental observations. The theoretical 
distribution is asymmetrical, with quite a high probability of the echo lower 
down the trail being greater than 5-5 times the upper echo. On the other hand, 
when the echoing points are above amax, there should only be 50 echoes with 
r<o'36 in Fig. 5 (6) (after allowing for experimental errors), whereas 273 are 
observed; and 20 with r <o-36 in Fig. 5 (a), whereas 100 are observed. If the 
observations are restricted to even larger ratios in Fig. 5 (a), there should be no 
echo pairs in which the echo from S, exceeds that from S, by more than 5:5 to 1, 
yet this occurs in approximately 10 per cent of the observations. 

(d) The mean ionization curve.—An attempt has been made to fit an ionization 
curve to the observed ratio distributions, and the resultant curve is shown in 
Fig. 4(6). When the procedure of Section 3 (c) is repeated using this curve, the 
theoretical distribution shown in Fig. 5(c) is obtained. This is in reasonable 
agreement with the observations, giving low values of r before &max. This 
proposed ionization curve is much shorter than the curve of Fig. 4(a), the 
difference mainly arising in the absence of a long exponential rise before the 
point of maximum ionization. ‘The differences are not so great near &max- It 
must be emphasized that the curve shown in Fig. 4 (4) is not necessarily a unique 
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solution, but it represents the average behaviour of meteors observed with the 
present equipment. 

4. Additional evidence for anomalous ionization curves.—So far we have 
considered only trails with «~ 10" cm~, using echo amplitude as a measure of 
electron line density. However, when 2210" cm, the electron line density 
can be determined approximately from the echo duration independently of the 
equipment parameters (Section 2(a)). Thus in Fig. 2(d) the echo at the home 
station is of the short duration “‘decay”’ type, and without specifying the precise 
value of « we can say that it must be less than 10! cm~?. From the rate of decay 
we deduce that the appropriate value of D is 3 x 10°cm*sec~!. At a point whose 
height is only 2-5 km lower the electron line density has increased to such an 
extent that a “long duration” echo has developed, with 7T’>0-3 sec. Sub- 
stituting a value of D allowing for the 2-5 km height difference, we find that 
a2>10!3cm-!, Thus the electron line density has increased by more than a factor 
of 10 in a height difference of 2-5 km. Echo pairs of this type are frequently 
observed, although they have not been used in the present analysis as they do not 
satisfy the restrictions placed on echo characteristics. However, they serve to 
show that very rapid increases in « do occur before «max, without reference to 
the aerial polar diagrams and other variables. 


TABLE I 
Number of diffraction maxima observed in meteor echoes at two spaced stations (V 25 
km sec"). In (a) and (b) the average reflection point is the same for O, and Og, and is 
approximately at &max- In (d) the average reflecting point is approximately 4 km 
further along the trail than in (c). 




















Number of diffraction maxima I max 2 max 3, 4 max >5 max 
per cent | percent | percent | per cent 

(a) All echoes at O, 52 21 18 9 
(b) All echoes at O, 55 16 19 10 
() Echoes at O, 

Coincidences with Vt> 4 km 35 19 31 15 

Echoes at Og. 
(d) SS ; , 

Coincidences with Vt> 4 km 73 13 13 I 























5. Fragmentation.—Very frequently an echo is observed from the upper 
reflecting point showing pronounced diffraction fluctuations during trail formation, 
yet these fluctuations may be much smaller or entirely absent in the echo from 
the point a few km further down the trail. In Table I the percentage of echoes 
showing well-marked diffraction patterns are given for all echoes on the two 
stations, whether coincident or not (a and 6). Only meteors with VS 25 km sec"? 
are used (estimated from the rate of rise of amplitude), as greater velocity results 
in a smoothing out of the fluctuations because of the low pulse recurrence rate. 
At both stations approximately 28 per cent of all echoes show at least three 
diffraction maxima while 53 per cent show only one. However, if we consider 
only coincidences with Vt>4 km, the situation is very different. The percentage 


a 
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of echoes with three or more zones is increased to 46 per cent from the upper 
reflection point, and reduced to only 14 per cent 4 km further along the trail. 
Similarly the percentage showing no zone formation is reduced to 35 per cent 
at O,, but increased to 73 per cent at Og. 

‘Table II shows that faster meteors are less affected, as a higher percentage 
of echoes give good diffraction patterns. ‘These results were obtained using 
earlier single station observations at a high pulse recurrence rate. 

This smoothing of the diffraction pattern can be explained by the progressive 
fragmentation of meteors proposed by Jacchia (1955) from observations of faint 
photographic meteors. ‘Thus each fragment will produce its own diffraction 
pattern, and when the separation of fragments along the trail becomes comparable 
with a Fresnel zone length the fluctuations will be smoothed out. Differential 
velocities of the order of 1kmsec™! are required to smooth out all but the first 
maxima in a path length of 1 km. ‘The first Fresnel zone is comparatively large 


‘((~1 km), and as this is primarily responsible for the main ‘“‘ body”’ echo the first 


diffraction maximum will not be seriously affected. 


TaBLe II 


Number of diffraction maxima observed in meteor echoes as a function of meteor velocity. 

















Number of diffraction maxima I max 2 max 3, 4 max > 5 max 
per cent per cent per cent per cent 
>60 km sec~?! 22 2 29 27 
3 
40-60 km sec~! 28 23 24 25 
25-40 km sec~? 44 20 22 14 
<25 km sec~! 53 19 19 9 























6. Irregular ionization.—The effect of irregular ionization along a meteor 
trail is to distort the simple diffraction patterns illustrated in Fig. 6(a), (6) and (c). 
The distortion is most noticeable if the irregularity occurs before the centre of 
the first Fresnel zone, and as an example the result of increasing the electron 
line density in the gth zone before S, by 50 per cent is shown in Fig. 6(d). Deep 
diffraction fluctuations are now observed before the meteor reaches the foot of 
the normal to the trail. These are analogous to the Doppler whistles observed 
with C.W. techniques as the meteor approaches the observer, although in the 
present case the meteor trail produces its own reference phase. Because of the 
small amplitude of these fluctuations they are normally observed only in trails 
giving very large body echoes, which saturate the receiver. As a selection 
criterion echoes which rise to saturation as the meteor crosses the first Fresnel 
zone have been investigated. Of these 51 per cent show deep fluctuations before 
the specular reflecting point, a further 29 per cent fluctuate slightly, and only 
20 per cent show the smooth rise in amplitude expected for a uniformly ionized 
trail. 

An example of this type of echo where the irregularity actually occurred 
between the normals to the trail from the two observing stations is given in 
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Fic. 6.—Theoretical diffraction patterns. 
(a) For infinitely long uniformly ionized trail, with zero rate of diffusion. 
(b), (c) As in (a) showing effects of increasing the rates of diffusion (Kaiser, 1953). 
(d) Effect of increasing the electron line density in the gth zone before Sq by 50 per cent, for trail 
of (a). (So is the specular reflecting point.) 


(e), (f), (g) Effect of sudden commencement of ionization (linear rise from zero to maximum electron 
density in a length (RA/2)"!*). (e) Trail begins 0-4(RA/2)!/? before So, (f) trail begins at So, 
(g) trail begins 0-4(RA/2)'!? after So. 


Fig. 2(c). ‘The echo at O, rises smoothly as the meteor approaches S,,, but at 
O, large diffraction fluctuations are observed before So. is reached. 

Occasionally an echo appears only as a small diffracted signal at O,, whereas 
at O, a normal body echo is obtained (Fig. 7); this effect can only be explained if 
the first Fresnel zone at S, is absent, and the meteor begins to ionize almost instan- 
taneously between S, and S, with a rise to maximum electron density in a 
distance ~(RA/2)"* (800m). Theoretical examples of the diffraction patterns 
for trails with abrupt commencements are given in Fig. 6(e), (f) and (g). 

A consideration of echo amplitudes and trail length shows that if all trails 
are of this type 3 per cent of the echo pairs should be of the form illustrated 
in Fig. 7. ‘Thus we might expect 50 examples of this effect out of a total of 1700 
suitable echoes, whereas only 8 were observed. This suggests that between 10 
and 22 per cent of all meteor trails reach their maximum ionization almost 
instantaneously, with an even greater rate of rise than indicated by the ionization 
curve of Fig. 4(6). 

7- Discussion.—An investigation of echo amplitudes at two spaced stations 
has shown that the ionization curves of faint +6 to +8 magnitude meteors are 
much shorter than predicted by theory. In particular the rise to maximum 
electron line density is more rapid than expected. It must be emphasized however 
that only a mean ionization curve has been derived, and large variations may 
exist between individual meteors. For example the experimental observations 
could be explained by a mixture of long and short trails. Direct evidence that the 
linear electron density may increase even more rapidly than suggested by the 
mean curve is given by the transition from short to long duration echoes from 
points separated by only a few km. Almost instantaneous increases in « from 
zero to maximum in a trail length less than 800 m have been inferred for a small 
percentage of meteor trails. Evidence for the fragmentation of meteors and 
irregular production of ionization has also been presented. 
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These results may be compared with the photographic observations of faint 
+3m meteors (Jacchia, 1949, 1955). In general the light curves are shorter 
than expected and breaks in the photographic trail produced by the rotating 
shutter become filled in, leading Jacchia to suggest that the meteors are in the 
process of fragmentation. In addition a small number of these meteors are observed 
to reach maximum light almost instantaneously, and others show irregularities in 
their light curves. ‘Thus the various anomalies observed in the light curves of 
faint visual meteors, appear to have their analogues in the ionization curves of 
faint radar meteors. 
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THE DISTRIBUTION OF THE DIRECTIONS OF PERIHELIA 
OF LONG-PERIOD COMETS 


F. G. Tyror 
(Communicated by R. A. Lyttleton) 


(Received 1956 October 31) 


Summary 


The distribution of the directions of perihelia of 448 long-period comets 
is investigated. Standard statistical tests strongly suggest a departure from 
randomness in the distribution. It further appears that there is some 
evidence that the perihelion points exhibit a preference for lying near the 
galactic plane. It is indicated how the accretion theory of the origin of 
comets can account for the form of the distribution. 





1. Introduction.—The problem of the origin of comets has engaged the 
attention of numerous workers over a long period of time. Although much has 
been discovered about the behaviour of comets, there still remains no generally 
accepted hypothesis of their origin, each hypothesis presented having some not 
entirely satisfactory features. 

Among recent hypotheses, Lyttleton (1) envisages a process of accretion of 
dust particles by the Sun in the course of passage through interstellar dust clouds. 
It is shown that this may result in condensations of interstellar matter falling in 
towards the centre of mass of the solar system along a general direction which 
must be that of the local relative motion of cloud and solar system. ‘These 
condensations are identified as comets, many thousands of which may, under 
suitable conditions, be formed on the passage of the Sun through a single cloud. 
The hypothesis explains many of the observed features of the behaviour of 
comets and of the properties of their orbits. 

One fact, however, emerges which is peculiar to this type of theory, namely 
that the comets must have originally come from preferred positions in the sky. 
The comets formed in any one passage through an interstellar cloud would be 
expected to have their perihelion points grouped together on the celestial sphere 
(if allowances for planetary perturbations could be made), which feature could 
only be seriously departed from if there was some change in the direction of 
relative motion of cloud and Sun during the passage. The effect of perturbations 
during the first and succeeding perihelion passages would probably result in the 
diffusion of the perihelion points of such a family, and also possibly in the loss 
of many members (by expulsion, break-up, or direct infall to the Sun). However, 
it may still be possible, at least for families born during recent cloud passages, 
to detect some traces of perihelion grouping. 

The object of the present paper is to examine the distribution of perihelion 
points of long-period comets to determine whether or not evidence of any such 
grouping is present. Previous work on this subject was carried out by Pickering 
(2) in 1911 (with comments by Eddington (3) and Oppenheim (4) in 1924). 
These writers, who preferred to consider the distribution of aphelion points, 
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claim to have found some evidence of the grouping of the points, but it has been 
thought advisable to revise their results in the light of the additional data on 
cometary orbits made available since their work was completed. Of the number 
of orbits included in the present investigation, 26 per cent refer to comets 
discovered since Pickering’s work was completed. 

Short period and periodic comets are excluded from this survey on the 
grounds that the large planetary perturbations that must have been experienced 
by them to bring them to their present orbits may be expected to have eradicated 
most if not all reliable evidence of any original grouping. In actual fact such 
perturbations must have affected all orbits to a greater or lesser extent. Hence 
some effort will be made to isolate and investigate comets that may have suffered 
only relatively slight perturbations since these will then be expected to give a 
better indication of the original distributions of cometary orbits. 

The data on comet orbits have been taken in the first place from Baldet’s 
tables (5) and, for more recent orbits, from the annual list given in Monthly 
Notices of the Royal Astronomical Society. Much of this has been summarized 
by Porter (6). 

The available data almost certainly suffer from observational selection, and 
so care must be exercised in attempting to draw conclusions concerning the 
general distribution of comets. ‘The nature of the available data has in its turn 
dictated the methods employed in the analysis. ‘These have been essentially 
simple for it was not thought worthwhile to attempt delicate techniques on a 
possibly suspect sample. 

2. Distribution of perihelion points.—Plate 7 contains a plot (in ecliptic 
coordinates) of the positions on the celestial sphere of the points of perihelion 
of 448 comets with periods greater than 200 years up to and including comet 
1952f. Of these 

67 have a calculated eccentricity e>1, 
263 have an assigned eccentricity e= 1, 
118 have a calculated eccentricity e<1. 


On direct examination several fairly large regions on the sphere are seen to be 
particularly noticeable for the absence of points, e.g. the so-called Hoek’s zone (7) 
centred roughly on (350°, — 30°), and the zone centred on (310°, 55 ). On the 
other hand it also appears possible to distinguish definite compact groups of 
points. This grouping becomes more apparent in Table I (a), for the construction 
of which the surface of the sphere is divided into 432 approximately equal areas 
and the number of points per area given. ‘Taking these basic units of area in 
groups of six, the observed frequency of points per area has been compared with 
that expected ona Poisson distribution. ‘This indicates a preponderance of regions 
with many points and with few points—a feature typical of the existence of 
groups. The value of x? associated with the observed distribution corresponds 
to a probability of less than o-oo1. Although differing values of the probability 
may be expected from different arbitrary subdivisions of the celestial sphere, 
this figure may be taken as an order of magnitude estimate, and is indicative of 
a real departure from randomness. 

In order to investigate further the distribution given in Table I, a simple 
sensitive test given by Sterne (8) may be used. Essentially this gives information 
concerning the size of the regions in which occurs the departure from randomness 
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previously noted. If these are larger than the subdivisions used, a tendency 
for unit areas containing small numbers of points to lie close together may be 
expected, and similarly for areas containing large numbers; i.e. there will exist 
a positive correlation between the numbers of perihelion points in a unit area 
and the numbers in other areas close by. The probability of a numerical 
configuration occurring through random causes can therefore be estimated from 
a correlation diagram such as that given in Fig. 1. Here the number (x) of points 
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Fic. 1.—Correlation between number of perihelion points per unit area (x) and average number 
of points in adjacent areas (y). 


per area is plotted against the average number (y) in the four immediately adjacent 
areas. Thus 72 points appear in the correlation diagram. The probabilities 
obtained again varied with the choice of sub-division and were O(10-*). This 
value is very low and in spite of the uncertainty involved in using the number 
of points in a given area in two ways a real departure from randomness is strongly 
suggested. It may thus be concluded that grouping must exist in regions larger 
than the unit areas considered. 

3. Plane of preference of perihelion points.—Coordinates may be assigned to 
the perihelion points by adopting a right-handed cartesian system with origin 
at the Sun. The coordinate axes are chosen so that the x-axis and the y-axis 
lie in the plane of the ecliptic (8=0) with the x-axis along A=o and the y-axis 
along A=7/2, with the z-axis completing the system. 

The sample moments of the rst and 2nd orders are then found to be 


x= —0'0069; y= —0°0446; z= +0°1316; 
a=x?= +0°2993; b=? = +0°3931; c=2%= +0°3076; 
f=yz= —0°00916; g=xz= +0°02741; h=xy = — 000997. 


From a visual inspection of Plate 7 it appears possible to identify on the celestial 
sphere a great circle of apparent maximum density of points. This suggests a 
form for the non-randomness of the distribution noted in the previous section. 
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Let the direction cosines of the normal to the preferred great circle be (/, m,n). 
Then 


ao 2 2 
D?= —s > (lat my; + nz;) 
448 4 


must be a minimum subject to the restriction /*+m?+n?=1 where summation 
is over all the perihelion points. ‘Thus 
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M= jh bf}. 
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The smaliest eigenvalue A=A, gives the preferred great circle. 

The significance or degree of definition of the preferred great circle may be 
assessed by calculating the value of D? and comparing this with the value of 
Di=| ss 
normal distribution of D? about its mean value of }, the value of ({ — D®) may be 
compared with the standard deviation of the distribution about the mean (see 


which is obtained from a random distribution of points. Assuming a 


Table II). In the case considered here (} — D®) is found to have a value of more 


than four times the standard deviation, and the preferred great circle is thus 
reasonably well defined. 


TaBLe II 


Mean square distance of perihelion points from preferred plane and galactic latitude of normal 
to the preferred plane 


D? on D? Standard g l—cos g 
Deviation 
All comets 0'276 0°057 0°014 78° 0°022 
Group corrected 
comets 0'269 0063 0'016 67 0080 
P class 0°226 0°107 0°036 84°°5 0°004 
M class 0°262 0'o71 0°026 81°°5 o'O12 


On substituting A=A, equations (1) solve to give 
l= +0°7674, m= +0°0151, n= —0°6410. 
This direction may be compared with that of the galactic pole given by 
L= +0°8772, M = —0°0536, N = —0°4772. 

The angle between these two directions is given by the galactic co-latitude 
(7/2—g)=12°. If the relation between the direction of the galactic pole and the 
normal to the preferred great circle were random, the probability of the two 
lying at an angular distance less than (7/2 —g) apart would be 1—cosg. In the 
case considered here the probability is seen to be ~o-o2. The low value of the 
probability certainly suggests the presence of a galactic influence in the distribution 
of perihelion points. Consequently the distribution has been obtained in terms 
of galactic coordinates and is presented in Table I1(b). The decrease in density 
of points to the galactic pole is immediately apparent. 
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The centre of the perihelion points is found to lie in a direction given by 

A= 261°, B=71 
which it is noted is only some 20° removed from the apex of the Sun’s motion 
at 

A,=271', — B,=53°- 
The significance of the similarity of the two directions depends on the significance 
of the displacement of the centre of the distribution from the centre of the sphere. 
The significance of the displacement can be determined by considering each of 
the first order moments of the distribution as independent with a standard 
deviation of 1/(1/3) and applying a y?-test with three degrees of freedom. In the 
case considered the displacement of the centre of the distribution is found to 
be significant at the 1 per cent level. 

The two effects identified here—a preferred plane and a displaced centre 
for the distribution—will not, in general, be independent of each other. It may, 
however, certainly be concluded that they are not different aspects of the same 
phenomenon. For, by itself, a displaced centre, i.e. the preference of perihelion 
points to be distributed about a particular point on the celestial sphere, will 
define an infinite number of equally significant preferred planes, all passing 
through this particular point. Further, the poles of such planes would be go‘ 
removed from the displaced centre. ‘The actual value found is 49°. 

It does not seem possible to correct accurately the preferred plane for the 
displaced centre. Generally, however, it may be expected that the result of 
such a correction would be to move the nearer preferred pole towards the 
position of the displaced centre, i.e. away from the galactic pole. 

4. Observational selection—The previous two sections show that there is 
some evidence that the perihelion points of long-period comets tend to be 
distributed non-randomly. ‘This tendency may be due either to observational 
selection, to a real grouping phenomenon, or a combination of both. That 
observational selection must play a large part in the discovery of comets has 
long been realized. For example Holetschek (9) has shown that a comet stands 
the best chance of being discovered if its perihelion occurs in opposition and 
has a perihelion distance q lying in the range 0-3 <<q<1-2a.u. Some writers (10), 
too, have defined “‘ zones of avoidance ”’ having the shape of right circular cones 
with apex at the Earth and axis the line Earth-Sun and within which, for small 
semi-angle of cone, discovery of even bright comets is highly improbable. Again 
a comet is most likely to be discovered at or near perihelion when, generally 
speaking, it is likely to be close to the Earth, and at its greatest intrinsic brightness. 

The last point could well explain the fact that, of the 448 comets considered, 
259 have their perihelion points in the Northern Hemisphere where there are, 
or have been, a preponderance of observing statinns. On the other hand this 
distribution of points between the two hemispheres may be an indication of a 
real preference for points to lie in the northern half of the sphere. In this case 
the work and conclusions of Sections 2 and 3 continue to hold good. If, however, 
observational selection is regarded as explaining any part of the discrepancy in 
numbers between points in the Northern and Southern Hemispheres then some 
account must be taken of it. 

The effect of this particularly simple form of selection, if it exists, can be 
surmounted by considering the observed distribution in zones of equatorial 
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latitude, since the random theoretical distribution can be assumed symmetrical 
about the polar axis. The distribution in the zones is found to be significantly 
non-random, and thus it can be claimed that the form of observational selection 
considered here cannot account for the previously noted departure from 
randomness of the distribution of perihelion points on the celestial sphere. In 
zones of equal galactic latitude the departure from randomness is far less marked, 
again indicating the presence of some galactic preference. 

The selection between the two hemispheres may be expected to affect the 
position of the preferred plane found in the previous section. Suitable 
corrections are found to move the preferred plane nearer to the galactic plane 
and also to shift the centre of the distribution closer to the Sun’s apex of motion. 

Allowances for other types of observational selection are more difficult to 
make but an attempt has been made by Bourgeoise and Cox (10) to consider 
correlations that may be contained in any cometary list. They conclude that 
certain statistical features of cometary orbits, such as the general distribution 
of numbers with g and of the observed mean value q with 7, are well explained 
on the basis of observational selection. In particular an ellipsoid of dispersion 
is obtained for the spatial distribution of perihelion points. The least axis of 
this ellipsoid is found to lie within 5° of the galactic pole. Earlier work by 
Oppenheim (4) is also corrected with similar results. It thus appears that even 
when observational selection is taken into account there still remains evidence 
of the non-random distribution of perihelion points. 

5. Groups of comets.—Local condensations of points can be observed in many 
regions of Plate 7, those at (275°, 0), (160°, —25°), (220°, — 30°), and (15°, 0) 
being especially noticeable. On closer examination it appears that in many 
cases the comets corresponding to adjacent points also exhibit a close similarity 
in their other orbital elements (with the exception of time phase), a fact first 
noted some go years ago by Hoek (7, 11) and later studied in detail by Pickering 
(2). The members of a group are thus comets that not only have come from the 
same part of the sky, but also in doing so have followed very similar orbits. In 
this way it is possible to discover 74 such comet groups containing in all some 
174 comets, although there is no absolutely clear-cut criterion for deciding 
whether or not particular comets form a group. The members of a group are 
usually considered to be the component parts of a larger parent comet and hence 
it is clearly necessary to consider whether the non-randomness in the distribution 
of perihelion points can be explained in terms of comet groups. 

Consequently the analysis of the previous section has been repeated after 
replacing each group by a single point at its ‘“‘centre”. Non-randomness is still 
present and details of the preferred plane for the group-corrected distribution 
of perihelion points is given in Table II. 

Mention has been made of the difficulty of deciding whether or not comets 
form a possible group, and the inclusion or exclusion of doubtful cases has 
considerable effect on the result. Corrections similar to those applied in the 
previous section can bring the preferred plane close to that of the Galaxy. 
Nevertheless, although there is a small indication of possible galactic concen- 


tration, the evidence derived by a consideration of groups must be regarded as 
inconclusive. 


6. The period of comet orbits.—Of the elements of the orbit of a comet, it 
appears that the period is the least likely to be affected by observational selection. 
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Consequently an investigation of the distribution of the periods of long-period 
comets may be expected to furnish immediate information about the general 
distribution of orbits in space. Perturbations by the planets will be likely to 
manifest themselves most strongly on comets whose orbits have a small angle 
of inclination to the ecliptic. It is well known, for example, that the short-period 
(captured) comets exhibit markedly small values of 7. Similarly it appears that 
comets with a period in years given by log P<6 may also have come to possess 
them through planetary perturbations. 

On the other hand, as Oort (12) has pointed out, if comets exist with large 
values of log P their orbits must stretch well out into interstellar space and hence 
be affected by nearby stars. Thus comets with a true period greater than a value 
given by log P= 16 may certainly be perturbed near aphelion by the Sun’s stellar 
neighbours. Hence there remains the possibility that comets with intermediate 
values of the period may be relatively unperturbed. It is interesting to note that 
the distribution of comet frequency with log P has a weak maximum in this 
region—at log P=8. 

Consequently the positions of the points of perihelion of the 70 comets with 
6-8 <log P<11-0o—which may be termed the P class—are considered separately. 
The preferred plane for the P class comets is calculated as in Section 3. As for 
the total distribution discussed in that section, it is noted that the plane of 
maximum frequency of points lies very near to the galactic plane but closer than 
does the same plane for the total distribution of perihelion points. In fact it 
is found that the normal to the preferred plane for the P class comets lies at an 
angular distance of only 5°-5 from the galactic pole (see Table II). 

7. Bright comets.—A further relatively unperturbed class can be expected 
to be provided by those comets that have made but few passages through the 
solar system and in this sense are to be regarded as new. It seems probable that 
high luminosity is associated in a general way with new comets. ‘The absolute 
magnitude at any time can be expected to be a function (decreasing with time) 
of the mass, density and chemical composition of a comet, and hence of its age 
and number of perihelion passages. ‘Thus on general grounds it seems that 
what are termed new comets may be expected to be brighter than the average. 

The question of how the brightness of comets depends upon distance from 
the Sun is an extremely complex one. The apparent law of variation of magnitude 
with distance from Sun and Earth varies from comet to comet and in addition 
comets are sometimes observed to undergo quite large sudden irregular fluctua- 
tions of apparent magnitude. As a result, no precise unique absolute magnitude 
(defined as the apparent magnitude if at 1 2.u. from the Earth and 1 a.u. from 
the Sun) can usually be calculated for a comet. Nevertheless, Vsessaviatsky (13) 
has assigned absolute magnitudes to several comets, and his results have been 
used here to obtain a list of some 79 bright comets, defined as those having 
an assigned absolute magnitude less than 5-00. 

In order to include in the discussion comets for which no absolute magnitudes 
are given, the list of bright comets is enlarged by including those further comets 
that were observable for relatively long periods—four months prior to 1875 and 
five months thereafter. The complete list of so-called bright comets thus totals 
127 and will be referred to as the M class. The two separate criteria of 
brightness adopted here are necessarily of a somewhat arbitrary nature and were 
chosen to provide a satisfactory number of comets for investigation. It is also 
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clear, however, that the members of class M were certainly brighter, on the whole, 
than comets not included in it. As before, details of the preferred plane for points 
corresponding to the M class are given in Table II. As with the P class it is 
found that the preferred plane for the M class lies rather closer to the galactic 
plane than the preferred plane for all observed comets. . 

8. Conclusions.—The danger of drawing too detailed conclusions concerning 
comets from a discussion of the type considered here must be immediately 
apparent. In the first place, doubt has often been cast on the reliability and 
significance of many of the observations made. Owing to the continuous effect 
of planetary perturbations, the elements of a cometary orbit will change during 
a single passage and from passage to passage. The data used in this paper refer 
only to particular osculating orbits that represent the actual paths of comets for 
limited times. Thus, for example, the given tabulated value of P= 1 385 600 
years for Comet 182611 will certainly not carry the significance the number of 
figures appears to indicate. 

Again, as has already been stressed, the available data have probably suffered 
from observational selection. Nevertheless, certain conclusions will be 
tentatively advanced. 

It seems quite clear, for instance, that departure from randomness in the 
observed distribution of perihelion points of the comets considered is strongly 
suggested, even when the effects of comet groups and observational selection are 
considered in an elementary way. When enquiry is made into the nature of the 
non-randomness it appears that there is evidence indicating a preference for the 
perihelion points to lie near a particular great plane. Further, the centre of the 
observed distribution appears to be significantly displaced towards a point on 
the celestial sphere. The possible identification of the plane and the point with 
well-known and easily recognizable astronomical features is certainly striking. 

The identification of the plane with that of the Galaxy also has been seen to 
become more apparent when two special sub-classes of the comets are considered. 
These represent comets with special values of the period (P class) and bright 
comets (M class), and were selected for special study on the ground that they 
might well contain a high proportion of relatively unperturbed comets. Some 
evidence is available in support of this. 

It is difficult, then, to resist the suggestion that the comets studied were 
introduced into the solar system with their orbital principal axes exhibiting a 
preference for lying near the galactic plane, and that superimposed on this there 
has been a further tendency for comets to approach the Sun from a direction 
(given by the position of aphelion) directly opposed to that of the Sun’s motion 
relative to the nearest stars. 

g. Accretion theory.—It is interesting to attempt to interpret these suggestions 
in the light of the accretion hypothesis. According to the hypothesis, the position 
of perihelion indicates the direction from which the comet-producing cloud has 
approached the Sun. The perihelion positions would therefore be expected to 
exhibit a tendency for lying near the apex of the Sun’s motion relative to the 
local interstellar clouds. If clouds share in the local general stellar motions 


this position may be expected to lie near the stellar apex of the stellar motion. 
As regards possible galactic influence, it may be observed that if cloud distribution 
or structure possessed any galactic preference then a star would remain longest 
immersed in a cloud when moving relatively in the galactic plane. Thus if all 
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directions of relative motion were equally probable the distribution of perihelion 
points would have a preference for lying near the galactic plane. 
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A PHOTOELECTRIC STELLAR SPECTROPHOTOMETER, 
USING A FABRY-PEROT ETALON 


4. E. Geake and W. L. Wilcock 
(Received 1957 March 7) 


Summary 


Tests are described of a photoelectric method of investigating stellar 
spectra, in which a prism monochromator and a Fabry-Perot etalon are used 
to scan the spectrum simultaneously. The addition of the etalon enables the 
monochromator to be used with a wider entrance slit for a given wave-length 
resolution. The etalon is tilted in order to scan in wave-length : this method 
is simple, but leads to a loss of wave-length resolution which is assessed 
quantitatively. 





1. Introduction —We have described elsewhere (1) an instrument for the 
direct photoelectric scanning of astronomical spectra. The light passing the 
exit slit of a prism monochromator is received by a photomultiplier. ‘This is 
followed by an amplifier and a pen recorder which draws out the spectral profile 
directly as the wave-length passed by the monochromator is varied continuously. 
In order to use our instrument at its maximum resolution of about } A a narrow 
entrance slit (about 25) is needed, and this width is usually several times less 
than that of the tremor disk of the stellar image. This leads to a considerable 
waste of light, and also to large fluctuations in the amount of light entering the 
monochromator as the image moves about under the influence of atmospheric 
turbulence. ‘The effect of these fluctuations has been removed satisfactorily 
for most purposes by a compensating system, but some imperfections remain. 

It would therefore be advantageous to widen the entrance slit until the whole 
stellar image is always passed, if this could be done without loss of resolution. 
It has been shown by Jacquinot and Dufour (2) that the addition of a Fabry-Perot 
etalon to a monochromator enables wider slits to be used for a given resolution, 
and references to the history of this method are given by Jacquinot and Chabbal 
(3) and by Dunham (4). It is the purpose of this paper to describe astronomical 
tests of such an arrangement which have been carried out with our instrument 
mounted on the 120 cm telescope at Asiago, Italy, and which we believe to be 
its first stellar application. ‘ 

2. Theory.—\f a Fabry-Perot etalon, having parallel reflecting surfaces with 
optical separation ¢, is placed in a parallel beam of white light, it passes a series 
of wave-lengths which satisfy the relationship 


2t 
di ineatebiiiepoaad p=I, 2, 3,... (1) 


where is the angle of incidence on the reflecting surfaces. It therefore acts as a 
filter with many passbands corresponding to different values of the order of inter- 
ference p. If all but one of these passbands is suppressed by a preliminary filter 
the whole arrangement has a single passband: we use our prism monochromator 
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as this preliminary filter. In order to scan the spectrum the passband 
of the system must be moved continuously over a range of wave-lengths. For 
the etalon this may be achieved by varying either ¢ or 6; it is better, for reasons 
given below, to vary t, but it is so much simpler to vary 6 that we chose to do this 
for our initial tests. ‘The wave-length passed by the etalon is changed at a constant 
rate, by using a synchronous motor and cam to tilt the etalon so that cos @ varies 
linearly with time. Simultaneously, the wave-length passed by the mono- 
chromator is changed, by using another synchronous motor to turn one of its 
prisms, so that the etalon passband remains in the centre of the monochromator 
passband. 

In our auto-collimating monochromator, the width of the exit slit is made 
equal to that of the entrance slit, which in turn is fixed by the requirement that 
the whole stellar image shall be passed. The wave-length range corresponding 
to this slit width must not be greater than the separation AA between the etalon 
passbands, which according to (1) is given by: 

AA= Ao _ Ao ; (2) 
Pp 2t 
and this fixes the maximum permissible value of ¢ for a given monochromator 
dispersion and slit width. 

The wave-length resolution of this combination of etalon and monochromator 
is determined by the half-width 5A of the etalon passband used. Neglecting the 
effect of the finite slit width, 5A is given by 

dA = . Ny (3) 
where .N is the equivalent number of interfering beams, and depends on the 
reflectivity and flatness of the surfaces. When the finite slit width is taken into 
account it is found to cause an increase in the half-width of the passband. This 
increase, 5'A, is negligible at normal incidence, but becomes increasingly important 
as the scan proceeds and @ increases: it may be found by differentiating (1), from 
which 





5A ~A, sin 0660 ~A, sin #é. ; (4) 
where s is the exit slit width and f is the focal length of the collimating lens 
immediately before the etalon. By (4) and (1), 5’A (for @ small) can be expressed 
in terms of the wave-length range (Ay — 4) of a scan starting at normal incidence 
as : 


; 2sA,)? " 
sAw~ 7 (Ay — Ag)". (5) 


This progressive increase of the instrumental half-width as the scan proceeds 
is the disadvantage of scanning by tilting the etalon. We were able to use this 
method because the slit width required to accept the image of a star was small 
and because we scanned only a small wave-length range. For a larger wave-length 
range the loss of resolution could be minimized by using a saw-tooth form of 
scan (3). 

3. Instrumental arrangement.—Fig. 1 shows the etalon mounted experimentally 
on our monochromator at Asiago, with the covers of its light-tight box removed. 
The tilting lever, cam follower and cam can be seen under the etalon box. The 
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angular range of the scan is adjusted by sliding the cam follower along the arm, 
and then moving the cam and motor assembly correspondingly. ‘This adjusts 
both the rate and the range of scan, as the cam throw and rate of rotation are fixed. 
The monochromator must be adjusted to scan in wave-length at the same rate 
as the etalon, so that their passbands remain coincident. Contacts, operated 
by the cam, stop the monochromator scan automatically at the end of the etalon 
scan, and restart it at the beginning of the next scan. In the interval the mono- 
chromator is turned back to a pre-set starting wave-length by hand. 

The etalon consists of a pair of 1} in. flats, of which the central 1 cm is used. 
The reflecting surfaces have 5-layer dielectric coatings; the present flats are of 
rather poor quality, and their irregularities limit the resolution to about 1/13th 
of an order. Nickel wire spacers go in diameter are used, giving an inter-order 
spacing of 11 A at 44471. The transmission of the etalon is about 60 per cent 
at this wave-length. 

The etalon may be turned by hand through go’, so that it can be adjusted 
to parallelism in situ by viewing circular fringes against a screen illuminated by 
a helium discharge tube. 

4. Tests and performance of the instrument. (a) Calculated performance and 
laboratory tests.—Table I gives values of the half-width in angstrom units, 
calculated from (5), after scanning 11, 22 and 33 A, starting at normal incidence. 
A collimating lens with a focal length of g cm was used before the etalon and 
a mean wave-length of 4471 A is assumed. 











TaBLe I 
Angle S Half-width in A at end of scan 
‘ ocan 
tilted deine 
Te ieee Slit width in p 
in A 
degrees 100 500 1000 1500 
4°1 II | 0°35 18 3°5 5°3 
5°7 22 | o-50 2°5 50 7°5 
6°9 33 0’60 3°0 60 gl 


























These figures are also shown graphically in Fig. 2; lines 1, 2 and 3 represent 
the calculated half-width after 11, 22 and 33A of scan respectively. Line M 
shows, for comparison, the half-width of the monochromator passband without 


fe} 


HALF-WIOTH 
IN AU. 


» w @& way @ © 
T a 





*) 








| 


Montuiy Notices of R.A.S. Vor. 117, Prate 8 








(a) (b) (c) (d) (e) 


Fic. 3.—Repeated scans of the Hy absorption line in the spectrum of scattered sunlight. 


Scans (a) and (b): Etalon, plus monochromator with 200 p slits. 


(c) and (d) : Monochromator alone, with 50 slits, giving about the same resolution as (a) and 
(b), with the amplifier gain increased by 3. 


(e): Monochromator alone with 2004 slits. 


Each scan took 4 min ; recording time-constant 1 sec ; range of scan 20A. 


J. E. Geake and W. L. Wilcock, A photoelectric stellar spectrophotometer, using a Fabry~Perot etalon. 
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(i) 
(a) (b) (c) (d) (e) (f) (g) 
(ii) 
(iii) 
Fic. 4.—Repeated scans of the Hy absorption line in the spectrum of 0 U Ma, of spectral class Go. 
Trace (ti): output of signal photocell, which scans spectrum. 
Trace (tii): output of comparison photocell, showing variations of total (white) light, after the 
entrance slit. { 
Trace (i): the ratio of (ii) and (iii), the spectral profile compensated for variations of the 
total light passing the entrance slits. 


Scans (a) to (c): Etalon, plus monochromator with 400 slits. 
(d): Monochromator alone with 400 slits. 


(e) to (g): Monochromator alone with 75 w slits, giving about the same resolution as 
scans (a) to (c), with the amplifier gain increased by 3. 
Each scan took 4 min: recording time-constant 10 sec: range of scan 20 A. 


J. E. Geake and W. L. Wilcock, A photoelectric stellar spectrophotometer, using a Fabry-Perot etalon. 
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the etalon. Experimental measurements were in good agreement with these 
calculated results except at the lower end, where the natural half-width of the 
etalon passbands (due in this case mainly to irregularities of the etalon plates), 
becomes important. The actual performance in this region is shown by the 
dotted lines. 

‘The measurements were carried out by setting the monochromator to pass 
the A 4471 line from a helium discharge tube. ‘The etalon was then tilted to pass 
this line at several successive orders, starting at normal incidence, and with the 
spacer adjusted so that the line was passed at normal incidence. The half-widths 
of the recorded profiles of successive orders were measured, for several values of 
monochromator slit width. ‘This is equivalent to measuring the half-width 
after successive intervals of 11 A, with the etalon and monochromator scanning 
together. 

(b) Astronomical tests.—(i) The Hy absorption line in the spectrum of scattered 
sunlight was scanned, and the result is shown in Fig. 3. Scans (a) and (6) were 
made by using the etalon and the monochromator with wide slits; scans (c) and (d) 
were made without the etalon, but with the monochromator slit width reduced 
to give the same wave-length resolution. The increased noise fluctuations result 
from the reduced amount of light passed by the narrower slits. ‘The amplifier 
gain was increased to keep the recorder deflection about the same. Scan (e) shows 
a scan with the monochromator alone with the same wide slits as were used for 
the etalon scans. 

(ii) Similar test scans were made of the Hy absorption line in the spectrum 
of o U Ma, a star of spectral class Go. ‘These are shown in Fig. 4, and illustrate 
both of the advantages of using the etalon. The wider slits used with the etalon 
result in a gain in light, in this case by a factor of about 3 in spite of the fact that 
the transmission of this etalon is only 60 per cent; also, the fluctuations caused 
by the stellar image moving about on the slit under the influence of atmospheric 
turbulence are greatly reduced. ‘This is seen by comparing the fluctuations on 
trace (iii), which is a record of the undispersed light after the entrance slit, for 
scans with and without the etalon. The fluctuations with the etalon are largely 
caused by changes in the atmospheric transparency due to haze, as these tests 
were carried out on a particularly bad night. It is the fluctuations due to 
obstruction by the slit jaws which are the most difficult to remove by compensation 
(x), and their absence results in the improved quality of the scans with the etalon. 

5. Acknowledgments.—The authors are grateful to Professor L. Rosino and 
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Summary 

The second-order differential equation governing the non-radial oscilla- 
tions of polytropic stellar models has been integrated for values of the poly- 
tropic index varying between 3 and 4. It has been shown that the 
fundamental mode is no longer a member of the discrete frequency spectrum 
for the polytrope n=3:25; and that, as m increases, an increasing number 
of modes on either side of the fundamental mode disappear from the spectrum. 

It has been shown analytically that the equation has no solutions com- 
patible with the boundary conditions of the problem for the limiting polytrope 
n=5. 





1. Introduction.—The significance of the problem of non-radial oscillations of 
stellar models has been adequately discussed in previous works by Rosseland (1), 
Pekeris (2), Cowling (3) and Kopal (4). Rosseland pointed out that a limitation 
to purely radial oscillations could fail to reveal the possible instability of a 
particular model for a more general type of disturbance: a consideration of the 
homogeneous gas sphere emphasizes this point. Kopal has stressed the interests 
of double-star astronomy in possible free periods of non-radial oscillations. 
It is well known that any non-vanishing orbital eccentricity (no matter how small) 
is bound to invoke, through tides, non-radial oscillations (corresponding to the 
second and higher spherical harmonic distribution) of both components. These 
are forced oscillations of period equal to that of the binary orbit; and if the 
latter is sufficiently far removed from any period of free oscillation in which 
either component can oscillate then the amplitudes of the forced non-radial 
oscillations can be predicted by the equilibrium theory of tides. If, however, 
one of the free periods can resonate with the orbital period then the amplitude of 
the corresponding mode would greatly exceed that of the equilibrium tide. This 
would give rise in eclipsing binary systems to photometric phenomena which 
could become not only observable but conspicuous. Such phenomena which 
could be ascribed to this cause have indeed been pointed out in a number of 
known systems; but their identification remains uncertain as yet. Asthe problem 
at issue is of considerable theoretical as well as practical interest, the aim of the 
present work will be to throw light on the problem of whether or not resonance 
between free and forced periods of non-radial oscillation is at all possible, or 
likely to develop, in self-gravitating gas spheres characterized by a sufficiently 
high degree of central condensation. In order to attempt to answer this question, 
our aim will be to establish the specific values of free periods of non-radial 
oscillations for polytropic gas spheres (of the many other models which could be 
considered) whose index n lies between 3 and 5, and to compare them with the 
Keplerian orbital period of binary motion. 
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In doing this, we propose to regard the amplitudes (and velocities) of oscillating 
motion to be sufficiently small for their squares and higher powers to be ignored 
{i.e. to consider the oscillatory problem in its acoustic approximation). The 
system of differential equations governing such motion is one of fourth order, 
and was first formulated by Pekeris (2). He solved his equations, however, only 
for the homogeneous gas sphere, and showed this to be unstable with respect 
to non-radial oscillations. It may be remembered that Emden (5) had previously 
considered the general oscillations of a polytropic gas sphere. He simplified 
the problem by neglecting changes in the gravitational potential produced by 
the oscillations, and thereby reduced the differential equation to one of second 
order; unfortunately, however, he used an over-simplified form of the equation 
of continuity. ‘This was pointed out by Cowling (3) who derived the correct 
second-order differential equation governing this simplified problem, and also 
investigated the general nature of its frequency spectrum. He also located the 
fundamental mode and one overtone for the second-harmonic oscillation of 
Eddington’s ‘‘standard model’’. Kopal (4) has extended these results to the 
third- and fourth-harmonic oscillations, and has located the second overtone in 
each case. These results obtained by Cowling and Kopal all apply to a model 
with a ratio of specific heats, y, of 5/3 and in this numerical investigation we will 
confine our attention to models having this value of y. 

2. Equations of the problem.—Consider the small adiabatic oscillations of a 
self-gravitating fluid globe. The motion of any individual particle is assumed to 
take place adiabatically, y being the ratio of the specific heats (assumed to be 
constant throughout the sphere). Viscous forces are neglected. Let the pressure 
P, the density p, and the gravitational potential V at any point of the initially 
spherical body be assumed to have increments 5P, dp and dV, respectively, due 
to the disturbance. Let the displacement of the matter from its equilibrium 
position be represented by h, which has a radial component 6r. The Eulerian 
equations of motion are then 


~ 


ae +V8P+8pVV +pV8V =0. (1) 


Cowling (3) has shown, by means of a perturbation method, that a neglect of V6V 
in this equation leads to values of the period which differ, in the case of the 
standard model, by less than 4 per cent from the true periods. The approximation 
should yield even better results for the higher degrees of central condensation 
which will be considered in this investigation, and it will therefore be retained 
throughout the present study. 

The equation of continuity is 


3p= —V. (ph); (2) 
and the equation of hydrostatic equilibrium 
VP=-—pVV. (3) 
The condition for adiabatic motion can be expressed by the equation 
dP d 
sP+arF = 2 (Spar), (4) 
where 
e=y-, 
”p 


¢ corresponding to the velocity of sound in an ordinary gas. 
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We now assume that the dependence of the perturbed quantities on the polar 
coordinates 6 and ¢ is given by a surface harmonic of order m; and further suppose 
that their dependence on time is a simple harmonic one. The period of the motion 
is taken to be 27/0. The former assumption provides the relationships 





app- 14 (atBP) _ mm) ' 
V8P= 45 (+ 7 = — 6P, (5) 
5p _1d/,dpdP\ | 
v.(2vP) = 55(+ ee) (6) 
and the latter gives 
oh 
ae = oh. (7) 


The foregoing developments have so far been quite general and the equations 
apply to any initial distribution of density and pressure. In order to obtain 
specific results, we propose from now on to restrict ourselves to a consideration 
only of the class of stellar models for which this distribution is polytropic. This 
means that in equilibrium the density and pressure, at a distance r from the 
centre, are connected to their central values p, and P, by the equations 


P=P"t!, p=p0" for o<n<s; 


where @ satisfies the Lane~-Emden equation 


1d { 2 dé 1 ae 

al" ix) +0" =0 (8) 
with boundary conditions 

6=1 and 4 =O at <%*#=0, 
dx 
n is the polytropic index, and x is given by r=r,x where 
rot = a . 
4nGp, 


The equations can then be simplified considerably if we introduce, following 
Cowling, the pair of non-dimensional variables z and w given by 


z= 8F pa snyiy w= SP pasmiy, ( ) 
To : A 9 
Under these conditions our fundamental equations can be shown (see Cowling (3)) 
to reduce to 


. (x*z)=gw, (10) 
and 
oe af; (11) 
where 
no [RtD 2) paeefer x] 
m yO |’ t) 





21+ 
pa 8 we K= (= -2) +m, “a= trot. 
c 


Y 





\ 
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‘The boundary conditions require a node (no displacement) at the centre, and 
a loop (no pressure variation) at the free surface of the oscillating globe, i.e. 
ér=o at r=o, 
$P=0 at r=n,, 
where p(r,)=0. ‘These conditions can only be satisfied for a certain discrete 
spectrum of values of «. here is, in general, a fundamental mode having no 
loops or nodes between those at the centre and at the surface. Cowling has 
shown, by applying Sturm-Liouville theory, that there is an infinite set of 
oscillations having indefinitely long periods and also an infinite set of indefinitely 
short period. Cowling referred to these as the gravitational (g-) and the pressure 
(p-) oscillations respectively. ‘They exhibit an increasing number of nodes and 
loops in the interior. The number of nodes or loops is denoted by the suffix. 
Thus the g,- mode has one node and one loop in the interior. 

3. Method of solution.—In the location of the eigenvalues of equations (10) 
and (11) use was made of the Manchester University Electronic Digital Computer, 
Mark I. We used a trial and error method identical to that used by Kopal (4). 
A number of trial integrations were made; the choice of « being designed to con- 
verge on the particular eigenvalue sought. We decided to use the Runge-Kutta 
method, as described by Gill (6), to integrate the equations. ‘This method can be 
applied to a set of simultaneous first-order equations in which the independent 
variable does not occur explicitly. We found it convenient to incorporate the 
Lane-Emden equation (8), in the form of two simultaneous equations 

a =0 and c= - =o, (12) 
with equations (10) and (11); and to integrate the complete system of equations 
(10), (11), (12) together with the further equation dy/dx=1. ‘This latter equation 
allows x to occur explicitly in (10), (11) and (12). The system is to be integrated 
subject to the initial conditions w=0, z=0, 06=1, 0’=0 at x=0. 

These equations cannot be integrated as they stand because both w and z 
are zero at the origin. The fact that w must be zero at the origin can be seen 
if we consider the series expansions, about the origin, of w and z. Kopal (4) 
has shown that close to the origin 

W = X"(dy + agx* + ayx4+...), (13) 

3 = x™-1(by + bax? + byx* +...), (14) 
where a, and 6, are related by ab)=ma,y. We will take a)= $a and b,=4m, for 
consistency with previous results. 

In order to overcome this difficulty we can either make use of the series 
expansions about the origin, derived and used by Kopal (4), to start the integration ; 
or we can make use of a transformation to a more well-behaved pair of dependent 
variables. ‘The latter course is clearly the more convenient for our machine. 


We put 


w= xl, (1 5) 
and 
sagt, (16) 
and equations (10) and (11) become 


dw, %1 _ 1% 
ye “a Ns (17) 
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dz, @, 2 
= =g— —-m— 18 
dx * x x (18) 
with initial conditions 
s, @ 
oe —=— at *=0. 
oa es 
Using a more systematic notation the five equations to be integrated are 
dy, at 
is 
dy, 
dx 
Os, 90 
dx x " 
dy, _ Ys ¥s 
— =f— —(m—1)— 
dx rt; ( Vs : 
Ya _ gY2 _ mY 
dx x a” 


I he initial conditions at x=0 are 
) I y of ) m 
Ye ities as a and ya 


=z 0 =JI . 
Yor M=t y Se os x 2 


Due to the general nature of the solutions, a variable step-length was used. 
The functions tended to vary most rapidly at the two ends of the range of 
integration and a shorter step-length had to be used there. This could be carried 
out automatically by the machine to keep a constant accuracy. The eigenvalues 
were located correct to four decimal places. 











| 2 


A A 4 rT 


+25 “50 +75 1-0 








Fic. 1.—Pattern of machine solutions for the polytrope n=3:25. 


The integrations were started with trial values of « Kopal found that if 
this value was too large the values of w near x, tended to + oo for the fundamental, 
the p,- or the g,- modes and to — oo for the p,- or g,- modes; while for 
small values of « the converse was true. The eigenvalue could thus be bracketed, 
to any desired degree of accuracy, by two values of «. 
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We first tried to locate the fundamental second-harmonic mode for the 
polytrope n=3:25. ‘The initial trial value was taken to be 0-25. In this case 
w, had two zeros close to the boundary and finally tended to + 00 (see Fig. 1). 
This implied that 0°25 was too small for the p,- mode and too large for the f- 
mode. We therefore made a second trial integration with «=0-24; and found 
that, although w, had no zeros near the boundary and tended to + © less rapidly, 
it had two zeros in the centre of the range. On further reducing « we found 
that the curve tended to that appropriate to the g,- mode. Thus for n=3°25 
there is no fundamental second-harmonic mode of oscillation: there is no mode 
which has no loops in the interior of the star (o<x<x,). This was also found 
to hold for the polytropes n=3-5 and n=3-75. Other small modes were also 
absent for these models. Fig. 2 shows the development of a minimum in the 
w, curve as m increases from n=2'5 to n= 3°25. 


fie 2:50 


M2795 


OS} 


N=3-00 


—=_ET —_———— 





Ne3 25 % 


— sn bn ne 


25 50 75 10 





Fic. 2.—The development of a minimim in the graph of w, as n increases from n=2°5 to n=3°25. 


4. The polytrope n=5.—We have seen in the previous section that as n 
increases an increasing number of modes on either side of the fundamental 
mode gradually disappear from the characteristic spectrum. It is natural, 
therefore, to enquire what happens when the degree of central condensation 
becomes infinite. This is true when the polytropic index reaches the value of 5. 
In this case it is well known that the solution of the Lane~Emden equation (8) is 

9=(1+ Jat) (x9) 
This polytrope clearly has an infinite radius but its mass remains finite. It is 
tantamount to a finite gas sphere whose whole mass is condensed to the centre, 
enveloped by an atmosphere whose mass is infinitesimal for gravitational purposes, 
i.e. the Roche model. For such a model the perturbation in the gravitational 
potential produced by the disturbance would be zero; so that the second-order 
system of equations (10) and (11) represents the rigorous formulation of the 
problem. 

In order to ease the formal analysis of this limiting case it is convenient to 
make a change of variables x=4/3tany; which has the effect of reducing the 


27 
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range of integration fromo <x < @too<y<n/2. The Lane-Emden solution (19) 
is then 6=cos y and equations (10) and (11) become 
dz 22 





ah 20 
dy r cos y sin y ~ ae 
and 
dw 
7 74% (21) 
where 
ss a . 
A=v/3(cosy)-? [2 + ry sin? y cos*y |, (22) 
and 
= 3 10 m(m + 1) ws 3 
B= Vy) Ee y cos*y ; (23) 


If we eliminate z between (20) and (21), the differential equation governing w 
assumes the form 





d*w dw 
dy? +007, +9q(y)w=0, (24) 
where we have abbreviated 
d,_tan*y 
P(y)= a ae ae (25) 
and 
q(y)= — AB. (26) 
Thus, near the origin 
2 m(m+1 
py) ~= and gy) ~— MEE, (27) 


Putting y,;=y—7/2, we have near the point y=z/2, 





291) ~£ and 9») ~- §,, (28) 


Equation (22) therefore has, in general, a regular singularity at the origin 
and an irregular singularity at y=7/2. This means that there can be no solution 
of (24) which satisfies the boundary condition that w must be zero at y=7/2. 
This statement follows from a much more general existence theorem (7). We 
will state here only the particular form of the theorem which applies in the above 
instance. 

If m, and n, are the orders of the poles which p(z) and g(z) have respectively 
at the point 2, which is an irregular singularity of the equation 


d*v d 
qa +P(2)5- +4(2)0=0, (29) 


then a necessary condition for the existence of a regular analytic solution of this 
equation, at z=2g, is that m, and m, should satisfy the relationship n,<n, +1. 

There are two exceptional cases for which the above proof does not hold. 
If «=o or y= 00, then g(y) would have a pole of order only two at y=z/2. The 
equation (24) would then have a regular singularity at both y=o and y=7/2. 
Neither of these cases is, however, of great physical significance as far as real 
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stars are concerned. ‘The former requires an infinite period for the oscillation, 
and this means that the star is in a state of neutral equilibrium: while the latter 
means that the fluid is an incompressible one. 

This proof applies equally well to the case m=o. The polytrope n=5 is thus 
incapable of performing either radial or non-radial oscillations of finite period. 

















TABLE I 

Polytrope Oscillation m=2 m=3 
3°25 pi 0°3129 0°3372 
f None None 

£ 0°1400 0°1704 

Le 0:0816 o0'1198 

3°5 Pe 0°3573 0°3870 
Pi None None 

¥ None None 

Li None None 

Re 00874 O-1174 

3°75 Ps 0°3795 0°4095 
: None None 

pi None None 

f None None 

21 None None 

Le None None 

23 00618 None 

£4 0°0442 0°0697 

















5. Results and conclusions.—The most important result which has transpired 
from our numerical work is the gradual disappearance of low modes of non-radial 
oscillation with increasing degree of central condensation of the respective 
configuration. For n>3 the polytropic family of models is seen no longer to be 
in a position to oscillate non-radially in the fundamental mode (i.e. in one for which 
the variations remain of the same sign throughout the interior), and with increasing 
degree of central condensation (i.e. as m5), an increasing number of adjacent 
*‘ pressure”’ and “ gravitational”’ overtones vanish from the discrete characteristic 
spectrum of free periods until, for n=5 (Roche model), the corresponding 
configuration loses completely the ability to oscillate non-radially in any free 
period at all. For n<5, the discrete frequency spectrum of free periods consists 
essentially of high modes of pressure oscillations (whose periods are very short), 
and correspondingly high modes of gravitational oscillations (whose periods are 
very long). The fact that an infinite number of indefinitely long periods of 
oscillation are consistent with the boundary conditions of our problem has been 
inferred by Cowling (3) from the Sturm-—Liouville character of our equations 
of motion, for very small values of «; the present results supply quantitative 
evidence as to the structure of characteristic periods corresponding to the low- 
frequency tail of the characteristic spectrum. ‘The density of this spectrum makes 
it, moreover, very probable that the periods of some of its members may come 
very close indeed to the period of the binary orbit and resonate with it. On the 


a7° 
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other hand, high modes of oscillation for which resonance becomes possible are 
characterized by so many nodes between the centre and the surface that their 
damping is likely to be severe, restricting the amplitude which they may 
eventually attain. 

As a result, it seems at present impossible to predict the height of the partial 
tides to which a close coincidence of the orbital period with one of the permissible 
free oscillations can give rise. The mere possibility, however, that this can 
happen should underline the need for a dynamical theory of stellar tides which 
is so far completely lacking, and whose development represents an important 
requirement of the theory of close binary systems. 
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Summary 


Transition probabilities have been computed for magnetic dipole and 
electric quadrupole radiation for all transitions between levels of the 3d‘ 
configuration of Fe v and the 3d® configuration of Fe 111. 





1. Introduction.—Many gaseous nebulae and peculiar stars have spectra 
showing emission lines due to forbidden transitions between the lower energy 
levels of atoms in the first three rows of the Periodic Table. A detailed study 
of these stars and nebulae presupposes a knowledge of the transition probabilities 
of the forbidden lines. Forbidden lines cannot be studied in the laboratory, and 
their transition probabilities must be obtained from quantum theory. Much 
work has been done on the calculation of transition probabilities of forbidden 
lines (for a general review see (g)), and for atoms in the first two rows of the 
Periodic Table the available data are relatively complete, although further 
improvement in the accuracy of the calculations may eventually be possible. 
The chief gap in our knowledge of forbidden line strengths of astrophysical 
importance is in the elements of the first long period. The atoms from scandium 
to nickel, and many of their ions, have low-lying metastable levels arising from 
configurations of the types d", d"-1s and d"-*s?. So far as the writer is aware, 
the only detailed study of the transition probabilities of the forbidden lines of 
atoms with incomplete d-shells is that of Pasternack (15) on the 3d? configuration 
in [Fe vit] and on the 3d* configurations in [Criv], [Mnv] and [Fevi]. The 
great labour involved in calculations of this type renders a study of all the atoms 
in the first long period impracticable. Iron is the most important of these atoms: 
iron lines of many different stages of ionization are present in celestial spectra. 
The lines of other atoms are usually much fainter, and for purposes of identi- 
fication and qualitative discussion transition probabilities for iron will serve for 
other atoms with the same electron configurations. With these considerations 
in mind calculations have been performed on the 3d‘ configuration in Fev and 
on the 3d® configuration in Fei. The results are reported in the present paper. 
The d‘* and complementary d* configurations contain several terms of a kind, 
and throughout this paper these are distinguished by using Racah’s seniority 
numbers as subscript prefixes, e.g. 3F, 3F. 

2. Method of calculation.—The first step in the calculation of forbidden line 
strengths is the setting up of the line strength matrices in LS-coupling. Shortley 
(18) gave explicit formulae for magnetic dipole line strengths in LS-coupling, 
valid for all configurations, and the numerical values required in the present work 
are obtained by substituting the appropriate quantum numbers in his formulae. 
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He showed that transitions only occur (in LS-coupling) between the levels of 
each term, and not between two terms. The only new point to be noticed is 
that this rule still applies even if there are several terms of a kind in a configuration. 
For electric quadrupole radiation the methods introduced into the theory of 
atomic spectra by Racah can conveniently be utilized. A detailed discussion 
has already been published (10), together with the numerical results for the 
d‘ configuration studied in the present paper. The d* results apply also to the 
d® configuration except that all the phases are reversed (18). 

The second step in the calculations is the setting up of the energy matrices 
for the electrostatic and spin-orbit interactions. Throughout the calculations 
configuration interaction was neglected. This is probably very satisfactory 
for Fe v, where the 3d*4s configuration is well separated from the 3d‘ configuration. 
The approximation is less satisfactory for Fe 11, where there is some overlap 
between the 3d°4s and 3d® configurations, but in the interests of simplicity the 
3d°4s configuration was neglected. The electrostatic energy matrix for the d* 
(and d®) configuration has been given by Racah (16) and by Trees and Harvey (32). 
The spin-orbit matrix for the d* (and d*) configuration was calculated* using 
formulae of Racah (16,17). This calculation was independently checked by 
D. E. Osterbrock. ‘Trees has shown (e.g. (31)) that the addition of a term 
«L(L +1) to the diagonal elements of the electrostatic energy matrix substantially 
improves the fit of the theory to spectroscopic energy levels. A term of this 
kind was included in our calculations. 

The next step in the calculations is the adoption of numerical values for the 
parameters occurring in the energy matrices. In accordance with normal practice 
the parameters are determined so as to obtain as good agreement as possible 
between the eigenvalues of the energy matrices and spectroscopic energy levels. 
It is not possible to treat all the centres of gravity of the terms as independent 
parameters (as has been done for simpler atoms) because of the occurrence of 
non-zero off-diagonal elements in the electrostatic energy matrix. The theoretical 
formulae were therefore retained. The electrostatic matrices involve three 
parameters, A, B, C, together with a fourth parameter « from the «L(L +1) 
correction. ‘The spin-orbit matrix involves one parameter, (,. The ion Fei 
has been studied in detail by Trees (30, 31). There appears little prospect of 
improving on his results, so that his parameters can be adopted without alteration. 
For Fev only the lowest five terms (seventeen energy levels) have been observed. 
The four parameters A, B, C and « were determined to give the best fit to the 
five observed terms. This method was checked by several trials on Fei; it 
was found that while the use of four parameters and five terms was not as good 
as the use of four parameters and all fourteen observed terms in Fett, it was 
much better than using three parameters A, B, C and five terms. The spin-orbit 
parameter (¢,,) was estimated by trial and error to give the best possible fit of the 
theoretical and experimental term intervals. The final values of the parameters 
ate listed in Table I. 

Using the adopted parameters the combined electrostatic and spin-orbit 
matrices were set up numerically, the sign of £, being changed to get the spin-orbit 
matrix for d® from that for d*. These matrices were diagonalized numerically 
by successive partial diagonalization, dealing with the large off-diagonal terms 
individually, and using perturbation theory when all the off-diagonal elements 


* The writer has deposited in the Library of the Society copies of the spin-orbit matrix for 
the d‘ configuration. 
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TABLE I 
Adopted Parameters 
Parameter Fe 11 Fe v 
A 19969 24303 
B 953 1144 
Cc 3652 4459 
a 67 85 
ba 429 510 


had become small. This method, which has recently come into prominence 
and has been used by many authors, is in fact due to Jacobi, and has been used 
in the classical problem of the secular variations of the elements of the planetary 
orbits (4). The final energy levels are given in Table II, together with the 
observed values (14). The transformation matrices were then used in the 
standard manner to transform the line strengths for both magnetic dipole and 
electric quadrupole radiation from LS-coupling to intermediate coupling. 

The electric quadrupole line strengths obtained by the above procedure are 
multiples of the square of the radial integral 

sy= | r® P*(3d) dr. 
0 
The evaluation of this integral requires a knowledge of the 3d radial wave functions 
for Fett and Fev. Professor D. R. Hartree has recently developed methods of 
estimating wave functions (with exchange) for atoms in the first long period. 
Using his wave functions, Professor Hartree estimates that for Fe1, s,= 1-28, 
and for Fev, s,=0-94. ‘These estimates have been used in the present paper. 

The final stage in the calculations is the conversion of the line strengths, 
obtained in atomic units, into transition probabilities. 

3. Results.—The final values of the transition probabilities for all possible 
transitions downwards from all levels up to 'F; are listed in Table III. For the 
nine highest energy levels, which are of little astrophysical importance, it has 
not been thought justifiable to publish in full the transition probabilities of the 
170 possible transitions from and between these levels.* For these levels, 
multiplet transition probabilities have been computed by summing the 
individual values over the lower levels and averaging over the upper levels: 

I 
nm (2L, +1)(28, +1) y 


where E, >E,,. The multiplet transition probabilities are listed in ‘Table IV. 

All quintet-singlet transitions (except *D-{G) have been omitted from 
Tables III and IV. For each atom there are 22 such transitions predicted as 
having finite strengths. The calculated values are very small, hardly greater 
than the rounding-off errors in the calculations, and no real significance can be 
attached to the values. 

It is virtually impossible to assess the accuracy of calculations of this nature. 
One criterion, unfortunately not a sensitive one, is the agreement of observed 
and calculated energy levels. Inspection of ‘Table II shows that this agreement 
is quite good for the terms as a whole, and excellent for the term intervals. An 
approach to the problem of assessing the accuracy of calculations of this type is 


A > (2J, + 1)A(L,, S.J» i LinSinJ m)s 


a on 





* The writer has deposited in the Library of the Society copies of the detailed results for all 
transitions for which the probabilities have not been published in full in this paper. 
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TasLe II 
Observed and calculated energy levels 
Fe v Fe 111 
Term J Observed Calculated Observed Calculated 
'-> eo ° —38 1027 952 
i. —- a ee. 84 
74 ed 274 193 202 
2} 419 8 377 7390 sa 
8 3 5 68 391 6 393 312 
3 04 a 7 ai © te 340 
4 | 1285 1261 ° — 103 aes 
aP oo | 24054 24034 21209 219i2 
4 a 
I | 24973 htt 24870 Be. 20688 ‘atin 21376 hs 
2 | 26466 493 26328 45 19405 3 20070 —_— 
3H = 4. |: 24037 ons 24811 7 20482 _—s 20520 sg 
5 | 25229 208 25120 309 20301 <i 20360 ee 
6 | 25527 25463 343 20051 5 20144 
3F 2 | - 26765 26822 21857 8 
4 8 21057 en 22081 - 
3 | 26846 Pe 26909 bei 21700 4 21910 Bir 
4 | 26973 27058 49 21462 3 21678 — 
5G = 3 |_-29820 29831 25142 248 
° 62 i -° se 4971 — 29) 
4} suse es Else 7 lasts | 7 
5 | 30429 30489 24559 3 24293 = 
DD 1 37070 30726 
nei 7 r 30119 
2 36932 m8 30716 a 30144 = 
3 36822 30858 4 30268 — 
ae 37500 30356 30427 
iG 4 37712 30886 30623 
iS 0 808 8 
4180 34812 34090 
iD 2 44950 35804 36732 
5 ee 52050 42897 42655 
-F 62 6 2 
2 1242 8 
3 6 4 2g | 5019S 110 i 80 
1271 __ 12, | $0295 he 50311 
4 61150 50276 19 50266 ——_ 
3P 6198 
2 1981 49148 7 49298 
I 61464 507 | 49577 — won a 
2 60484 ~~ 98° 835 856 
404 50412 50571 
3G 4 69916 57222 57390 
1 
2D 2 92241 76019 
iS o 118975 98065 
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Tasce III 


Transition probabilities of [Fe 111] and [Fe v] 


All entries in this Table are probabilities of spontaneous emission from the upper to 
the lower of the levels mentioned in column 1. The relative positions of the two levels 
involved in any particular line should be ascertained by reference to Table II. Units : 
reciprocal seconds. The letters ‘ v.s.’ denote a line with a theoretical finite strength, 
which is however too weak for transition probabilities to be calculated by the methods of 
this paper. 


Fev Fe 111 
Transition 
r A, An A, 
5pD—D o-I 1°6 x 10-4 oes 14x 107-* ae 
o-2 = Jaxto™ : roxX10"" 
I-2 raxw 6°5 x 107" 6-7 x 10-* 3°5 x ro-* 
1-3 se 6-0 x 107° _ ss X10-" 
2-3 2%6 x 10 74x 10°" 1°8 x 107% 6:0 x 107"! 
2-4 ~ ra xXs0"° ; Es x10" 
3-4 3°0 X 107% 2°6 x 107" 2°8 x 107° 4°7x 10" 
*D-{P I-o 1°3 ae 0°67 
2-0 : 28 x 1074 a 2°4 X 107 
o-I 0°13 0°09! eb 
I-1 9 x10 1a X20-* 8 x 107% 7 xrT07* 
2-1 I'l ar X10~* 0°53 73x t0° 
3-1 4°1 x 1075 3°3x to-* 
o-2 ns 4°15 x10" us m5 X 10° 
I-2 0°036 3°1 X 107° 0:038 7°3.X 107° 
2-2 79x 10~* rexXs ry x20" 3°5 <s0°* 
3-2 o’7I1 raxs0"* 0°40 6°5 x 107° 
4-2 - 1-5 X10°° 7-1 x 107° 
$P-3P o-I O°O14 a 75x 107° 
o-2 a 6-1 X 107° shi 8-4 x 107° 
I-2 0°045 1°3 X 107° 0°047 16 x 107° 
5D-*H 2-4 ee V.S. ae 1°6 x 1078 
3-4 4 xX1074 v.S. 8-6 x 1074 2°7x 107° 
4-4 rr X 107° 5°8 x 107” 4°8 x 107% v.S. 
3-5 a 10 X 107° ss 2°8 x 107* 
4-5 : x<ro™ v.s v.s 3°3 x 10° 
4-6 8-5 x 10-* 47x18" 
3P-H 2-4 gas 2°5 X 107° aioe 1°3 X 10710 
*H-*H 4-5 6-5 x 10-4 6°9 X 10714 1'9 X 1074 5°3 X 10718 
4-6 se 1:0 X 1078 33x z0"" 


5-6 5°8x 10-4 4:1x10"* 4°1 X 1074 2°5 X 1071# 
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Tas_e III (continued) 


Fe v Fe 111 
Transition 
A, A, Ay A, 
°‘D-j]F 0-2 wn 2°2x 1075 ins 9°3 x 107% 
1-2 o'10 2°6 x 1075 0049 1°3 xX 1075 
2-2 0°20 v.S. o'10 V.s. 
3-2 0°047 g°0 X 10° 0026 8-5 x 107% 
4-2 zs X 107° a°IX10° 
I-3 vit 8-6 x 107-8 ae I*r xX 1075 
2-3 0°16 1°8 x 1075 0°087 a4 x 10~° 
3-3 0°40 20x 1075 0°27 6°5 x 10-8 
4-3 0:066 2°6 x 1075 0°038 2°7 x 1078 
2-4 eae 8:6 x 107? rae 7°6 x 107% 
3-4 0°16 rsx 107° 0'081 2°4X 107° 
4-4 0°74 1°3 x 107° 0°44 6°8 x 10-5 
jP-jF o-2 = 1°9 X 107? , 2°5 X 107° 
I-2 9°4 X 10-* 27x10°° 3°70 X 107° 3°S x t0~* 
2-2 3°6 x 1077 7°3 x 10-"* 20x 10°* rs X10°° 
I-3 ee 3°8 x 10° es 3°38 x 10°* 
2-3 15 X1077 9°0 X 107” 2°9 X 107-8 7°9 X 1078 
2-4 7°97 X 1071 2°0 X 107? 
5H-} 4-2 - 6°7 x 1078 a 2°8 x 1078 
4-3 rex 30" 8:2 x 107° 74x<eo"* I°o x 19078 
5-3 a 4°5 X 107 ‘ 1°8 x 1078 
4-4 4°I X 1073 2°2 X 1077° 1°3 xX 107° 1'o x 107° 
5-4 5°8 x 1074 §°3 x 10-* $°9X10-* 7°0 xX 107* 
6-4 a°§ X 07" 2°0 x 1078 
}F-}F 2-3 1-3 X 107% 17 xX 10—% E*3 X 1078 2°9X 1075 
2-4 2°4 X 10715 16x 10°" 
3-4 4°1 X 1075 1-1 X 107** 3°4x 1074 rrxIz0°" 
‘DG I-3 =n 2°9 X 1075 “i 2-2 X 107% 
2-3 7°0 X 107% 5°8x 10-5 2°7 X 107% 2°3 x 1075 
3-3 0'O17 g*I X 1078 8-0 x 1073 2°1 xX 1075 
4-3 26x 10° v.s. 7°9 X 1074 5°9 x 1078 
2-4 tee 4°7 X 1075 ion 2°4 x 1075 
3-4 7°38 x 107-3 1°1 X 1074 3°7 X 107% 4°9 Xx 1075 
4-4 0'032 3°6 x 10-8 ©'O19 4°7 X 1075 
3-5 os 5°0X 1075 aia 2°8 x 1075 
4-5 3°5 X 107 13 X 1074 2°4X 10~* 9°3 x 1078 
$P-°G 1-3 si 3°8 x 10-8 ids r9x10-* 
2-3 4°6 x 107% 5°9 xX 107° V.S. #°-3 X 107° 
2-4 4°7 X 107° I°7 x 107? 
3H-*G 473 0°036 5°9 x 1075 0°044 7°4X 107% 
5-3 te 2°3 X 107% Sa 1°5 X 107° 
4-4 0°033 7°6 x 10-8 0°025 mrxX10~ 
5-4 3°9 X 10-4 56x 1075 1°5 X 1074 6-3 x 1075 
6-4 ke 1°7 X 1078 wis 5 xX 207° 
4-5 1°2X 107% 8-5 x 1078 2°7 xX 1074 1°77 X 107? 
5-5 O'O4I 5°2 X 10~* 0029 8-0 x 10-8 
6-5 0'041 6-1 x 1075 0042 6-9 x 1075 
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TABLE III (continued) 


Fe v Fe 111 
Transition 

Ay A, a A, 
j7F-°G 2-3 0°030 4°4X 107? 0°025 1°r X 107% 
3-3 0'037 8-1 x 107® 0°039 4°6x 107” 
4-3 22x10°* 5-5 X<s07° 1°§ X 107° 978x107" 
2-4 . 4°4x10°° - 2°4x 107" 
3-4 73x 30"° 5°6x 107" 3°2X 107° 7:0 X 1077 
4-4 027 ¢5.x10"° 0°037 1'2X 1077 
3-5 ‘ 3°0 X 107° i 1°6 X 107’ 
4-5 0°037 9°2 x 107? 0°023 8:5 x 1077 
3G-G 3-4 9°2 x 107-4 2°9 X 10738 2°6 x 1074 5°5 x 10" 
3-5 = rsxi0°™* _ 1°8 x 10713 
4-5 4°4X 1074 2°4 x 10738 1°4 X 107% 1°5 X 1071 

5pD-*D o-I 0°22 er o'l ; 
I-I o'19 2°1 x 1074 O'ls 1°5 X 10-* 
2-1 1°9 X 107% 4°2X 1075 1°3 x 107% 3°3 x 107° 
3-1 4°4 x 1075 4°1 x10 
o-2 5°7 x 10°° aes 3°8 x 107% 

1-2 0°20 v.S. 0°095 v.S. 
2-2 018 2°4 xX 1074 Ovll 2°2 xX 10-* 
3-2 o'r! rx 1075 0'027 17 x10-° 
4-2 I°9 X 1074 z7 xX 107° 
1-3 = 2°6 x 107 ane 2°1 x 107% 

2-3 0'097 V.S. 0°044 v.S. 
3-3 0089 3°4x 1074 0°047 3°0X 10-4 
4-3 0°37 9°: X 1074 0°23 5°8x 10-4 

;P—"D o-I 0°049 ans 0'023 ve 
I-I o'12 2-2x10°°* 0085 1°5 X 107 
2-1 0'036 2:2 x 107° 0'019 2X 207° 
0-2 oe 1°6 X 10-4 a 4°5 X10 
1-2 7°8 x 1074 3°1 X 1075 26 x 1074 1°4 Xx 1075 
2-2 052 70x 10° 0066 27x 10~* 
I-3 9°4 Xx 1075 in 6°7 x 1075 
2-3 0°056 8-3 x 1075 0'063 2°8x 1074 
3H-°D 4-2 6-4 x 1075 = 2° xX 107% 
4-3 5°8 x 1075 7X io" 1°0 X 107 5°§ X 107° 
5-3 rx t0"* 16x 107 
jF—-D 2-1 O°O14 x*7< 207° 49X16" 1°6 xX 107° 
3-1 as 3°4 X 107-4 ios 1°4 xX 107% 
2-2 O17 6:8 x 1074 20x 107% 7°2 X 1074 
3-2 1°6 x 107% 1°4X 107% ‘2X 107% 1'°0X 107% 
4-2 as 2°0xX 1074 “= 1°0 X 107% 
2-3 5°6x 10-5 4°3 x 107 1°39 Xx 107° 56x 1075 
3-3 6°4 x 1073 4°9 X 10~* 3°8 x 107% 57 X10 
4-3 6-9 X 1073 1°9 x 107% 6°6 x 1073 1°9 X 107° 
GD 3-1 in 6°5 x 1074 on 22x 107* 
3-2 9°4 Xx 1075 4°0X 107° 4°9 X 10-* 6°5 x 107% 
4-2 os 3°7x10—* a 2'2 x 107 
3-3 1°8 x 10-* 3x10" V.S. 71x10 
4-3 8-3 x 107% 1°4 X 1075 gsxX10"° 8-0 x 1075 
5-3 2°5 X 107* hen 4°3 X 1074 
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Transition 

3pD—*D 2-1 
3-1 

3-2 

3H-"T 4-6 
5-6 

6-6 

*F-"I 4-6 
3G-1T 4-6 
5-6 

5D-j}G 2-4 
3-4 

4-4 

3P-1G 2-4 
3H-jG 4-4 
5-4 

6-4 

jF-jG 2-4 
3-4 

4-4 

3G-1G 3-4 
4-4 

5-4 

3D-1G 2-4 
3-4 

1[-1G 6-4 
¢P-jS I-o 
2-0 

3F-!S 2-0 
°D-1S 1-0 
2-0 

3P-1D o-2 
I-2 

2-2 

3H-1D 4-2 
3F-1D 2-2 
3-2 

4-2 
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m 
1°0X 1074 


3°3 X 107° 


ovr! 
O°14 


2°8 x 1074 


7 x107* 
K zo" 


o'18 
0°25 


O's 
0°32 


0°042 
6°4 x 1073 
0'026 


3°6 x 1078 


3°3 X 1075 


0-062 
o'18 


O'2I 
0°42 


Fe v 


2°3 X 10714 
9°6 x 10715 
1°6 x 190715 


I°I X 107° 
2°1 X 1078 
2°5 X 107° 


“4X 1077 
2x10" 


a4 s67* 


y 1078 
< 1075 
1°6 x 1075 


— = me 
AdL 
_— 

° 
=x 


7°6 Xx 1078 
1°5 X 10714 


1°7 X 10718 


7-7 X 1075 


8-8 x 107% 


78 7 107? 
8-1 x 1075 
9°6 x 107° 


3°r X r0~* 


3°3 X107* 


R. H. Garstang, Transition probabilities for 


A 
3°9 x 1078 


s°0X107° 
0°058 
0089 
x 107* 


< 10-4 
K 107° 


N 
N 


x 1074 
¥ 1074 
KX 107* 


w 
Ww A) 


1°3 x 10712 


x 1075 


*033 
‘096 


oo 


° 


‘12 


Fe 111 
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A, 
5°3 X 107 
2°5 x 10715 
2°6 x 1071 


1'2 xX 107% 
5°9 \ 10~’ 
1°8 x 107% 


2°3 x 10° 


x zo~* 
x 107" 


1°5 X 107° 


3 x 107% 
x 1078 
1°7 xX 1075 


<10°° 
< 1077 
<10 7 


oo 
west 


x 107" 
xX 1077 
x 1078 


N 
wnt a 


3°4 . 10716 
1°r X 1078 


7 10718 


5°3 x 10°* 


6°7 xX 1077 
4°9X 10> 
3°7 x 107~° 
3°1 X 10° 


8-2 x 1075 


2°0 X 1078 
9°8 x 1075 
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TasLe III (continued) 


Fe v Fe 111 
Transitions 
An A, An A, 
°G-iD 3-2 9°5 X 107% 2°6 x 1075 2°5 xX 107% 20x 1075 
4-2 1°4 X 1075 aes 3°0x 107° 
*‘D-}D I-2 0080 i: Xz0"* o'021 27x10" 
2-2 0'017 2°5 X 107° 3°7 X10 5°6x 107° 
3-2 0°090 6-5 x 107° 0°020 9°9 X 1077 
4G-4,D 4-2 pick 2°0 X 107 ios $°3Xz0~ 
i5-1D o-2 3 1°9 X 107° a 2 x26"° 
jP—"F I-3 ; 1:4 X 10-4 zs X z0-* 
2-3 1'o xX 107% 4°7 X 107° 5°4 x 10-4 20 x 1078 
3H-'F 4-3 4°0 X 107% rs x10" 4°¢x<10°° 8-9 x 1074 
5-3 1°8 x 107% 2°8 x 107° 
3F-1F 2-3 2°09 x 107* 1°2 X 107° 4°4x 1074 1°o X 107° 
3-3 5°1 X 107° 7°8 x 107 I°g X 107° 1°3 X 107° 
473 o°o12 1°2xX10~° 5°6 x 107% 4°8 x 10-* 
3G-'F 3-3 orl 72x 10~* 0'050 20 X 1074 
4-3 O'l 10 X 107 0°098 73 xX 107° 
5-3 6-3x 107° 2°axs0—* 
®D-1F 1-3 ae 1°3 X 1075 _ 8:2 x 107% 
2-3 0°070 9°9 X 1075 0°034 7X ro 
3-3 O'15 4°0 X 1075 0'061 1°7 xX 107% 
iG-'F 4-3 v.s. 4°8 x 107% v.s. 3°7 X 1073 
,D-'F 2-3 v.S. 3°2x107-4 v.S. 5°9X 107 
TABLE IV 
Transition probabilities for multiplets 
Fev Fe 111 
Multiplet 
Am A, bs A, 
‘D-3F 0°16 5°9 X 107% 0'077 4°4 x 107% 
<P-3F v.S. 0°22 V.S. o'18 
3H-3F 27x 107% 1°8 1:2 x 107% 1°3 
3F-3F O'15 0°89 0:096 0°66 
*G-3F 0'27 0°83 O'15 0°57 
8P-3F 0°26 O°OI5 0°097 go X 1073 
1]-3F - 6-0 x 1078 wi 12x10" 
{G-3F 0022 10x 107% 6:0 x 107% 7°0X 10-* 
19-3F one 1'1 X 1078 _ 1:2 x 107° 
,D-3F 0°035 74x10 * 0'016 “sxs0— 


iF-3F 0081 8-7 x 10-* 0048 5°9X10-* 
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Multiplet 


sp_is 
2F-1S 
°D-1S 
1D-1S 
SFIS 
sP-1S 
1D-Is 


* Transition probabilities cannot be quoted for the 2F-3P multiplet because the two 
terms overlap. The transition probabilities of the individual lines are all small, not greater 
than 10~* sec~, and are unlikely to be of any importance. 
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V.S. 
v.S. 
o'19 


0°043 
O'14 


3X s0°° 


V.S. 
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70 X 1073 
3°3 X 107% 
0°20 
oOrll 


0°44 


3°3 
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Fe v 


0°O14 

I'o 

7°2 X 1075 
0°13 

3°1 Xz0°° 
o's! 


1°6 x 107% 
6°6 x 107% 
2°0 X 107° 
4°7 x 107% 
2°3 X10~° 


0'020 
or 
0°43 

33°0 
1°4X 107% 
0'032 
3°2 


Fe 111 


0°36 

0°022 
V.s. 
v.S. 

0°097 


0'026 
o'll 
7°2 X 107° 


O'14 
o'19 
0°39 
7°5 x 10-8 
4°6 X 1074 


I'i X 1073 
0°040 


0°055 
2°7 X 107% 

V.S. 
o'61 


1°7 x 10~* 
O'10 
0'061 


o'I2 


I°9 
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0°010 
082 
7°8x 1074 
o'r! 
17x 107° 
0°30 
3°3 x 10-* 
6-1 x 1075 
2°4 xX 1074 
4°4x10°° 
* 


1°6 x 1078 
3°8 x 1073 
9°9 X 1074 
4°5 X 107 
ms x<z0"° 
I°4 

0°13 
0'021 
0°044 
16x 107% 
1°7X 107% 


0°027 
Ovll 

O's 
1°4X 107% 
0°036 
I1°4 

0°96 

ce 

0°50 

4° x 107° 
4x 207° 
0°29 


0°021 
0-090 
0°29 


25°9 


6-6 x 107* 
0023 
2°3 
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to repeat the calculations for various reasonable values of the parameters, and 
examine the changes in the transition probabilities. This approach is not 
practicable, however, as the labour involved would be prohibitive. Previous 
experience with simpler atoms suggests that the strengths of the stronger lines 
should be accurate to better than + 50 per cent, but the weaker lines may be 
uncertain by larger amounts. 

4. Relative intensities in multiplets.—It is interesting to examine Table III 
and the data on which Table IV is based to see whether any general rules are 
suggested regarding relative intensities in multiplets. Among the triplet-triplet 
transitions are several in which electric quadrupole radiation is stronger than 
magnetic dipole radiation, the former being permitted in LS-coupling for these 
transitions, the latter arising only through departures from LS-coupling. In 
these cases transitions with AL=AJ are the strongest. This is as expected by 
the usual formulae for relative intensities in electric quadrupole multiplets. For 
the remaining triplet-triplet transitions and for triplet—quintet transitions 
magnetic dipole radiation predominates, and in these cases there appears to be 
no rule by which the strongest transitions can be predicted. For the singlet-triplet 
transitions (AS= +1) the combinations (AL= —1, AJ =o), (AL=o0, AJ = +1) 
and (AL= +1, A/ = +1) give the strongest lines, but there are a few exceptions 
to these rules. ‘There is no obvious theoretical reason why any particular rules 
should hold; the intensities depend on the precise nature of the departures from 
LS-coupling in the configuration concerned. 

In his paper on 7 Carinae (28) Thackeray pointed out that for his observed 
[Fe11] transitions between terms of the same multiplicity AJ =AL gives the 
strongest lines. ‘This suggests that most of these transitions are due predominately 
to electric quadrupole radiation. 

5. Comparison with observations.—[Fe 111] lines were first identified by Edlén 
and Swings (7) following their analysis of the Fe 111 spectrum (8). ‘They showed 
that five strong lines of previously unknown origin in RY Scuti (11) could be 
ascribed to the °D-F multiplet of [Fe11]. The strongest line was present in 
a number of other objects. [Fe] and [Fev] were tentatively identified by 
Dufay and Bloch (5) in Nova DQ Herculis 1934. The first star found to have 
strong [Fev] lines was AX Persei, studied by Swings and Struve (22, 23). 
Many objects are now known with either or both [Fe 111] and [Fev]; among them 
are RY Scuti (8, 11, 19, 26), BF Cygni (1, 12, 23, 24, 27), AX Persei (22, 23), 
RR Telescopii (28, 29), 7 Carinae (28,29), FR Scuti (3), MWC17 (20), 
MWC 349 (21), Z Andromedae (13), R Aquarii (20), the Orion Nebula (25, 33) 
and NGC 7027 (2). [Fe 111] and [Fev] have also been observed, usually blended, 
in various novae, and other peculiar stars and nebulae show the one or two 
strongest lines of these ions. 

The most extensive data available for comparison with the theory are for the 
5‘D-3F multiplet of [Fer]. Much’ of the available data for this multiplet are 
collected in Table V. The theoretical intensity is on a scale on which 4658 
is taken as 100, and it is further supposed that a Boltzmann distribution is valid 
between the three levels ?F,, ?F;, 3F,;. The observed intensities are quoted 
directly, no attempt being made to place them on a common scale. The general 
agreement between theoretical and observed intensities is remarkably good, and 
may be interpreted as confirming the theoretical intensities. It also suggests 
that there are no very substantial departures from the Boltzmann distribution 
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TABLE V 


Intensities of [Fe 111] multiplet 7D-{F in celestial spectra 


























J 
51) 2F A Theory |1 2 3 gh >? & -o 6 
° 2 | 4800 | 0-0012 I 
I 2 | 4778 6 ° 0°83 «1-2 o'7 
2 2 | 4734 13 1 4 Oi + ¥ Os § OO 2 F 
3 2 | 4667 3 ° bl : 6 o% 
4 2 | 4574 | 070003 | Oo? 03 
I 3 | 4814 | o-0019 bl bl 
2 3 | 4769 15 : 4 £2. @ os 2 08 I°5 
3 3 | 4702 48 314 68 5 221008 7 P 
4 3 | 4607 ” ° 9 o? 8 o°5 o's 
2 4 | 4824 | 0°0017 6 
3 4 | 4755 18 16 25 3 I o9 82 04 
4 4 | 4658 100 “20 39% 8 gs @ 15% 3h %.20 
1. RY Scuti (19) 6. RR Telescopii (28) 
2. RY Scuti (26) 7. Nova Herculis 1934 (6) 
3. BF Cygni, 1950 July 15 (1) 8. Orion Nebula (33) 
4. BF Cygni (23) g. Orion Nebula (25) 
5. FR Scuti (3) 10. NGC 7027 (2) 


in the °F term. Lines of the *D-$P multiplet of [Fe11] and of the *D-3F and 
>D-3P multiplets of [Fev] are observed in various objects, and the theoretical 
and observed relative intensities are in general agreement. Special mention 
must be made of observations of the spectrum of RRTelescopii in 1956. 
Dr A. D. Thackeray kindly sent to the writer prints of an excellent plate showing 
the [Fe v] spectrum in great strength: all the strong predicted lines of the °D-{P 
and ®*D-{F multiplets are present with approximately correct relative intensities. 

Several other points should be noted: 

(a) The low transition probabilities of the °D-°H and *D-°G multiplets 
provide a reasonable explanation of the weakness or absence of these multiplets 
in the various celestial spectra (e.g. (22)). 

(6) The absence of the *D-°D multiplet of [Ferm] in BF Cygni (24), the 
Orion Nebula (8) and other objects cannot be explained by low transition 
probabilities, the values of the latter being comparable with those for observed 
multiplets. 

(c) The identifications of 44800 in BF Cygni (23) as the *D,-3F, transition 
of [Fer] and of 43736 in AX Persei (23) as the *D,-3F, transition of [Fev] 
are not supported by the theoretical probabilities of these transitions. 

(d) The identification of 45057 in various objects as the ®D,-%P, transition 
of [Fe f11] is not confirmed. 

(e) The identification of 44986 and 5033 in RY Scuti (19) as lines of the 
*D-*H multiplet of [Fe 111] seems unlikely in view of the low probabilities of 
these transitions. 
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we both were at Yerkes Observatory. Thanks are due to Professor D. R. Hartree 
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PHOTOELECTRIC PHOTOMETRY OF GALACTIC CEPHEIDS 


Olin }. Eggen*, S. C. B. Gascoigne and E. 7. Burr 
(Received 1957 January 28) 


Summary 


Photoelectric light curves have been observed in two colours for 26 galactic 
cepheids, and these, together with another four variables observed only near 
maximum and minimum and the 33 reported on in a previous paper, bring the 
total to 63 for which photoelectrically determined light curves are available. 
This number includes all the known, nearer cepheids. The relations 
between shape of light curve, amplitude and period length are re-examined 
and it is shown that the magnitudes used by Blaauw and Morgan in their 
discussion of the zero-point of the period-luminosity relation are too faint by 
o™-17, relative to the International Py-system. 





1. Introduction.—A previous application of photoelectric photometry to 
cepheid variables (Eggen 1951) showed that little reliance could be placed on 
much of the photographic photometry of these stars except for the determination 
of periods. The photoelectrically-determined light curves indicated, among 
other things, a division of the classical cepheids into three sub-groups, but the 
need for more observations was clearly indicated. A second programme was 
accordingly commenced in 1951 by Eggen with the g-inch Mt Stromlo refractor 
and continued by him with the 12-inch Lick refractor, and by Gascoigne and 
Burr with the g-inch refractor and the 30-inch reflector at Mt Stromlo. In this 
programme emphasis was placed (i) on the brighter cepheids, since these would 
presumably be the least reddened by galactic absorption; (ii) on those of the 
longest periods, for comparison with similar cepheids in other galaxies; and 
(iii) on cepheids that might be expected to show population type II characteristics. 
The new programme has resulted in reasonably complete light curves for 26 
cepheids which, together with an additional four variables observed only near 
maximum and minimum light, bring the total number for which photoelectrically 
determined light curves are available to 63. This number includes most of the 
known nearer cepheids. 

The details of the observational procedure are discussed in Section 2 and the 
individual observations are set out in Table II and plotted in Figs. 1 to 4. The 
maxima and minima colours and magnitudes and other data relating to all 63 
cepheids on this and the preceding programme are collected in Table IV. In 
Section 3 we discuss the correlations between shape of light curve, period length, 
and the light and colour amplitudes, and it is shown that the variables with 
asymmetric light curves have greater amplitudes, in general, than those with 
symmetric light curves. We also conclude that, because the magnitudes used 
were in the mean too faint, the zero-point of the period—luminosity relation 
derived by Blaauw and Morgan (1954) needs a correction of —o™17. 


* Formerly at Lick Observatory, now at the Royal Greenwich Observatory. 
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2. Observations.—All of the colours and magnitudes in this paper are on the 
(P, V)z-system which has previously been fully described (Eggen 1955); for 
convenience, (P-V), has been abbreviated to Cy in the figures for this paper. 
The observations by Gascoigne and Burr were made with a 1P21 photomultiplier 
and photometer, on loan from the Lick Observatory, that are almost identical 


TABLE | 
Comparison stars for cepheid variables observed by Gascoigne and Burr 

Comp. Name Ve A.D. (P-V)z A.D. n 
HR 1917 8 Dor 5°33 +o-o1 0°93 + 0:005 4 
HD 69066 RS Pup 6:88 03 0°67 ‘Ol 7 
HR 3445 RZ Vel, T Vel 3°78 "02 0°565 "O15 5 
HD 93306 VY Car 747 03 0°95 ‘Ol 7 
HR 3914 1 Car 5°51 02 1°21 ‘Ol 6 
HD 108500 R Cru, T Cru 6°82 "02 0°54 ‘OI 7 
HD 109761 R Mus, S Mus 7°39 02 0°81 ‘Ol 7 

HD 161852 W Sgr, X Sgr, 
RY Sco 6-60 "02 0°25 "02 5 
HR 7130 x Pav 5'97 + 0°03 0°90 +or-or 6 


with those used by Eggen; their filters were 2mm of Schott GGu1 for the 
yellow region and 1 mm of BG12 with 1 mm of Chance OY ro for the blue, 
the latter glass being used to exclude the ultra-violet region. The observational 
procedure used by Eggen has been previously described (1955). Gascoigne 
and Burr compared the cepheids with the stars listed in Table 1; the magnitudes 
and colours of these stars were determined through observations of some bright 
stars and stars in Harvard E-Regions Nos. 3, 6 and g, which had been observed 
at the Cape Observatory (1953), and the reduction to the (?, V),-system was 
made with the relations previously published (Eggen 1955). 


TABLE II 


Observations of cepheid variables; those designated by an asterisk were made by Gascoigne and 
Burr, the remainder by Eggen 


SW Tauri 42 19™-3, +3° 54’ SW Tauri (continued) 
Max =JD 2419730°3536+ 19-5836468 E JD Phase V P—V 
Pp m m 
JD Phase V P-—V 2887'809 0327, 9°35 0365 


2860-899 0°335 9°36 0°365 


Pp m m 

2434418°743 0041 9°73 (0°45 
4334'944 0°126 9°38 0°36 0°331 9°36 
3178°947 07167) «9°36 = 0°35 3180967 0443 9°67 = 0545 

4425°736 0°457 9 9°71 = -0°545 





0°365 








0146 §=99°37, "355 
0450 899°69 = 0545 














2882-840 0-189 9:28 80°30 
4369°889 O'192 9°36 = -0°3.45 4341°934 0540 9°78 0-60 
2882-869 0:208 9‘22 0°34 3190°933 0736 9°89 0-705 
2883-759 0°779 =. 993. :0"705 

0'196 9°29 0.33 
2882-986 o'281 g‘21 o'3I °°753 hs re 
2887-736 0281 9°38 0°32 4424°710 0809 =«0"rO 0°705 
2887°751 0290 9°36 0°32 2883838 o'819 10°00 0°725 
0'284 9°32 0°315 0814 10°05 o’71I5 


28* 





Olin }. Eggen, S. C. B. Gascoigne and E. }. Burr 























408 
SW Tauri (continued) 
JD Phase V P—V 
p m m 
2883-865 0'900 10°06 0°69 
3206°853 0996 «9°89 0°57 
RX Aurigae 4" 54™-5, + 39°. 49’ 
Max=JD 2420651°543+119-623731 E 
JD Phase V P—V 
Pp m m 
2433216'917 o'o10 7°34 0°68 
4425 °840 O'Ol5 7°40 o'7I 
4460°705 O'O!l5 7°39 0°68 
0013. 7°38 ~—(0 69 
4428-809 o°185 7°41 o°75 
4335°049 0204 7°45 0°825 
4369°941 0°206 7°44 0805 
0198 7°43 «0°79 
2882-976 o0'281 7°50 0°83 
2883 °952 0°365 7°55 =: O'915 
4337°016 0374 «97°58 0895 
0°370 = 7°56 0"905 
2861005 0°391 7°64 0'935 
4418-844 0°413 7°70 =—-0"95 
0402 «7°67 80°94 
2885-995 O54 7°94 = =I1°O1 
2944°903 0609 8-06 1‘00 
4341986 o'801 7°81 0-905 
3181°031 0°923 7°63. :0"795 
4424 "901 0934 «= 7°63. (0785 
0928 7°63 0'79 
SV Persei 4" 42™-8, + 42° 07’ 
Max=JD 2419055°345-+119-128877 E 
JD Phase Vg (P-V)z 
Pp m m 
2434424°783 0°04! 8°65 080 
2944°896 0°064 865 o-81 
07052 865 ~— so -Bo5 
3178-962 0°096 8-70 0°85 
4425°769 0135 «8°74 O'S 
4337°012 0154 8-76 o-gI5 
0'144 8-75 o-gI5 
3181-023 0-281 8-95 1°00 
2947°791 0°224 9°00 1°00 
4418°816 0°460 89°31 39099 
2860°865 O°513s_- 9°33 1°02 





SV Persei (continued) 


JD 


2883-938 
4341°958 


2885 °979 
2887 °947 
4335 °042 


B Doradus 
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Phase Vz P—Vz 
Pp m ™m 

0°587 9°30) = 0995 
0°599 9°30 1°02 
0593 9°30~—C—i=éd''021 
0°770 8-99 0°865 
0°947. 8:54 = 0°67 
0977. 8°55 = 0°72 


g® 3a™-7, —6a° 33' 


Max=JD 2426013:93 +. 9%:842603 


JD 


2434705 °147* 
4420°080* 
4449°958* 
4834°016* 
4451°009* 
4844°952* 
4422°044* 
4451°957* 
4826°083* 


4855 °936* 


4433°031* 
4442°963* 


4462°974* 
4434°050* 
4700°158* 
4463 °994* 
4847°942* 


4828-917* 
4820°003* 
4829°985* 
4840°981* 
4438°008* 
4419°030* 
4438-960* 











Phase V P-V 
p m m 
0°020 3°48 = 0058 
0058 3°54 0°60 
07093355 (0°63 
O°113 3°58 0°64 
0°'200 3°69 o’7!I 
0°224 =©3°79 8°75 
0°257 3°79 = 9°75 
0'296 3°83 o's 
0°307 3°86 080 
0°302 3°84 0805 
0°340 3°91 0°82 
0°373. 3°93 0°89 
0°383. 3°97. 089 
0378 §=63°95 += 089 
0416 4°00 0°84 
0°477 4°05 0°86 
0°513 4°04 0°84 
o'519 4°02 0°89 
0°528 4°00 0°88 
0°520 4°02 0°87 
0°596 §=— 3°93. 0°76 
0'689 3°78 0°67 
0°704 3°75 0°66 
o'821 3°58 0°65 
0°879 3°56 054 
O'95I 3°39 «86046 
0°976 §=93°40 39052 


RT Aurigae 65 22™-1, + 30° 33’ 
Max=JD 2420957°466 + 39-728261 


JD 


2434342°049 
4427°819 





Phase Vg (P-V)z 
Pp m m 
0°034 «55°08 )—(0'3 
0°039 «= 5§'°07-—s«OO'B1 
0°036 = 508 )—o'531 
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RT Aurigae (continued) RS Puppis (continued) 
JD Phase Ve (P-V)z JD Phase V P—V 
p m m Pp m m 
3249892 0°094 518 0°38 4834°026* O°145 6°82 1°36 
2825-012 0°132 5°21 O41 4840°051* 0°290 7°02 1°48 
3238-8909 O'145 5°22 0°43 4840°991* 07313. 7°05 ~=—«I"50 
4335062 o'160 8=6§°27 = 0465 4469'100" 0°330-7°07”—sI'st 
4424865 0247 5°39 (0505 3764979 0-331 7°05 ~—sa'52 
2833°028 0282 5°41 O°515 
2944°931 0207 5°44 (0°55 0-330 7°06 8515 
2822007 0°326 5°47 0°57 ‘ a 
3285°906 = 0465, 5°59 (064 — — = we 
4425°856 0°513, 5°70 «= -0°675 4244°035° 0386716159 
2848-962 0556 85°74 0°65 4844'978*  o409 87°16 1°57 
4459°728 0598 §8=95°72 = 0665 4472°083* 0405 719 1°53 
3020°741 0630 §©=— 5°75 = (0675 
4418-885 0°643 5°75 0°68 0°407. 7°18 1°55 
2882-903 0659 $°97 0°655 4847-985* 0-482 7°32 1°57 
2841°975 0682 5°77 0°675 3775°915 0595 7°53 1°61 
3437°052 0694 «5°79 — (0685 3776-906 0619 7°57 1°60 
3777°943 0°644 7°59 = 1°63 
0688 ‘78 0°68 3778°917 0°668 = 7°53 1°S7 


ec: * . . . 
77 -o-6ss 4855°979"  0°675 7°58 161 


‘41 oar 0672 7°56 1°59 
‘26 = -0°365 4859°003*  0°748) 7°34 1°45 
4820°'014* 0806 6°84 I‘I9 
4451°988* 0°920 6°45 1°05 


2860-984 0°780 — 
3233°920 
4460°760 0875 
2883896 0926 





° 
io.) 

° 
UuUuUuUun wm 
~I 
rs 








AH Velorum 8! 08™-9, —46° 21’ 4826°087* 0°953 6°49 1°08 
Max =JD 2426145°708 + 44-227172 E 4867 -964* 0'964 6°55 ae ce) 
JD Phase Ve (P-V)e 0958 652 1-09 


p m m 
2433775892 0°033 5°49 0°385 RZ Velorum 8" 33™-6, —43° 46’ 

LL laa dl pe od 0°49 Max=JD 2421057°710+ 2049-39645 E 

4454°110% 0°47 5°85 0°50 , > 
Ph P- 

3777°909 "510 583. O"50 ‘i = ; m 

° x . . . ) 

5841°212 ° one 5°76 pa 2433765135 0'021 6°49 «=: 066 

3778901 0°745 5°75 0°495 4826°1 I 2“ 0°039 6°48 0°64 

3783 °993 0950 5°47 (O41 , 





0°030 )=-« 6485065 
RS Puppis 8" o9™-2, —34° 17’ 


. 867°977* 0092 6° 078 
Max = JD 2427083°613+419-414 E ee ° 57 7 


4848 -002* O'112 6°65 0°80 








JD Phase V P-V 4808-151* o'159) — 674—Ss«O0'838 
p os _ 4830°022* 0°231 6°84 1°03 

2435160-025* O'OI7 6°58 1°22 4871 -996* 0289 6-91 1°12 
4828-931* 0°022 6°63 1'22 4873°008* 0°338 7°07 1°18 
4833°062* 0380 7°10 1°24 

0:020.~ = 6 ‘60 1°22 4834°056* 0°429 7°21 1°25 

4829°996% 0°047 6°64 «1°27 we pont 0"554 7°48 I 38 
4871°987* o-061 6°64 = 1°28 toe ph sc 7st nh 
872972" 0085 668 1°34 ae oe ee 
re 5 3778°040 0°654 7°61 1°40 
3796°896 o'102 6°74 1°32 3779010 0702 «= 762 1°28 
4833°046* o-121 6°75 1°33 4820°036* 0-741 7°51 1°25 
4844°081* 0°920 7°25 1°02 

Orrr = 6°74 -1°325 4844°907* 0965 645 0°59 








410 


T Velorum 8) 34™-4, —47° 00’ 
Max =JD 2426274°432 + 41-639769 


JD 


2434872-001* 
3797°005 


4826-107* 
4840°'058* 


4867-984 * 
4844°985* 
4859°015* 
4872°995* 


3778°070 
4808°158* 


3796°915 
3779°015 
3765 "140 


4855°991* 
4833 °069* 
3776010 

4820°022* 


4834°035* 
4847°992* 


4830°003* 
4844°069* 


SW Velorum 08 40™-4, — 47° 
Max=JD 2434872-95 +239-473 


JD 


2434872 °999* 
4830°007* 
4808 -170* 
4833°072* 
4834°040* 
4840°063* 
4844°988* 
5481°238* 
4847°9905* 
4872:004* 


Olin J. Eggen, S. C. B. Gascoigne and E. }. Burr 


7 Carinae 




















Phase V P-V JD 
Pp m m 
0'O17 7°71 0°62 2434872-‘010* 
0018 =©7°67 = 0605 4873°030* 
3775986 
018 ‘6 6 
rhe ee 4557°960* 
0°125 7°81 0°73 
0132 «67°88 073 
‘ ' 3776-988 
— I~ we 4808 -116* 
O'IS5I 7°84 0°75 
0194 7°93 «0°78 
0218 7:97. 0°77 wns 
0231 7°97. 0°83 3778°035 
; 3779°003 
“—-— 7 = 4845004* 
0°245 7°99 0°785 
0°257 8:02 o'8o 
4527°064* 
o'251 8:00 = -00'79 
3852°859 
0°305 8-10) 085 4848-016* 
0448 815 0°88 
0°455 8:16 0875 
3784°010 
0°452 8-16 0°88 3855°877 
0°566 8:28 o'"9I 4818-067* 
0°626 = 828 —o'95 4820°043* 
0800 8:25 3087 4469°187* 
0814 8-09 0°75 4826-201* 
0°834 8:07 0°74 3797°O15 
0°842 8-o1 0°75 4472°166* 
0838 =68'04—0'745 
0965 7:69 0°60 4830°028* 
0996 «7°68 ~=— 063 3765 097 
4867 -989* 
4833 °084* 
02’ 4087°153* 
E 
Phase V P-V 
Pp m m 
0°002 7°55 0°66 
O'171 §=67°80— «005 VY Carinae 
0240 8= 788 1°02 
0°301 7°95 1'14 
0342 -7'99—Ss«16 JD 
0°5909 846-147 
ms 8 27 2434833 °099* 
0914 8-70 = 1°22 4928-072* 
0°937. 849—Ssar12 
0960) = 80g—ssoO'TE09Q1 
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gh 42-5, ~62 03, 
Max =JD 2427089°82 + 354-5323 E 
Phase V P-V 
p m m 
0017. 3°34—s«0'8Q 
0°046 3:28 0°89 
o'r7I 3°44 1°04 
0°179 3°40 1°04 
O°175 3°42 1°04 
0°200 3°44 1°09 
o'219 3°47 1°09 
0'210 3°46 1°09 
0°229 3°49 I*I2 
0'256 3°52 1°175 
0257 3°52 118 
0°256 8= 352 1°18 
0°309 3°60 1°20 
0°335 3°63 1°29 
0°342 3°63 1°29 
0°338 «= 33°63,—s«1'29 
0°397 3°69 1°32 
0°420 3°75 1°335 
0499 «3°34 sa°35 
O'555 3°91 1°38 
0°681 3°99 1°35 
0°722 »= 401 r3a 
0°764 4°02 1°305 
0°764 3°99 1°36 
0764 4:00 1°33 
0835 4:00 1-29 
0865 3°07 = 1°25 
0904 3°84 I*I4 
0922 3°76 1°08 
0°929 3°75 1°I25 
0°926 «=. 3.76 1‘10 
10) 40M.6, — 57° 02’ 


Max = JD 2434529'90 + 184-9370 E 





Phase V P-V 
Pp m m 
o’oll 6:90 90°74 
0'026 6°89 o'71 
0018) §=66'90)— 0°72 
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VY Carinae (continued) 


JD 


4834°136* 
4873°018* 
3851°900 
3777°O15 
3796-920 
3779°025 
4557°981* 
4845°018* 
4560°989* 
4920°971* 


4543°988* 
4848-023* 
4527°076* 
4472178" 


S Muscae 





Phase V P-/Y 
p m m 
0:066 7°03 0°86 
O'lIg 7:06 0'90 
0°197 7°27 1°06 
0°243 7°31 1°07 
0°2904 7°44 I‘ll 
0°349 7°49 1185 
0483. -7°79°—s«d1"19 
0°640 7°90 1°25 
0°642 7°94 1°21 
0°651 7°88 5°23 
0°644 7°QI 1°22 
0°744 7°76 I‘10 
°0°799 «7°65 — 1°03 
o'8s1 7°64 0°96 
0'952 6°87 0°67 

12h o7M-4, 69° 36’ 


Max=JD 2419541°800+94-658688 E 


JD 


2434833°142* 
4544°030* 
4921°003* 
4950°000* 
4834°163* 


4873°066* 
4950°979* 


4845 °093* 
4903 °060* 


4527°110* 
4952°949* 
4953°943* 
3852°943 

4558°031* 


4848 -086* 
3776°118 


4926°043* 
3777095 
4927 °009* 
4917°985* 
3855 °888 
4561°019* 

















Phase I P-V 
p m m 
o0°169 5°91 0°63 
0°237 5°93 0°63 
0'266 5°93 0°70 
0°268 5°97 0°66 
0°275 6:00 = 070 
0°270 5°97 0°69 
0°303 6°03 0°70 
0°370 6°14 0°72 
0°407 6°15 o°75 
0°408 6°14 0°75 
0408 6:14 0°75 
0°485 6°21 0°79 
0°574 635 0°84 
0676 «66°41 = 0°83 
0686 6°40 8=6_ 082 
0686 6°45 0°83 
0°683 6°42 0°83 
0°717 6°45 0°79 
0°723 6°40 80°78 
0°720 6°42 0°785 
0°788 = 6°37, 0°75 
0°833. 6:20 80°67 
0888 6°13 0°68 
0954 «5°98 0°59 
o'991 5°96 0°60 
0996 = 55°95 058 
07994 5°96 = 0°59 


T Crucis 
JD 


2434925 008* 
4845 °069* 
4925°949* 
4952°969* 


3855 °891 

4873°045* 
4953°934* 
4927°017* 


4833°119* 
3776°183 

4920°992* 
4544°015* 


4558°020* 
4834158" 


3777 186 

4848°054* 
4949°992* 
3852°950 

4950°967* 
4527°096* 
4561 °003* 
4917°999* 


R Crucis 


12h 15-9, 


411 


-61° 44’ 
Max=JD 2427578-96 + 69-733178 E 

















Phase V P-V 
Pp m m 
0'022 6°36 0°67 
O'150 6°41 0°74 
o*162 6°43 0°76 
O'175 6°43 0°75 
0°162 6°43 o'75 
0239 «645,079 
0°305 6°52 0°83 
0°318 6°54 0°85 
0°32! 6°52 0°82 
O°315 6°53 0°83 
0°375 6°57 0°89 
0°40! 6°56 0°875 
0°426 6°63 0°88 
0°438 6°62 o"90 
0°432 6°62 0°89 
0518 6°70 = 0905 
0°529 6°69 0°93 
0°524 6°70 0°94 
0°550 6°69 0°935 
0°593 6°72 0'92 
0°733 6:80 0-92 
0°803 6°66 0°84 
0°878 6°56 0°73 
0°925 6°48 0°70 
0-961 6:40 80°66 
0-981 6°36 = 0 68 
0'971 6°38 0°67 

12) 78M-1, —61° 04’ 


Max JD 2419432'010+ 59-825589 E 


JD 


2434544°007* 
4527°091* 
4848-058* 
4917°994* 
3776°185 
4952°963* 


3776-986 

4953°928* 
4925°015* 
4558°015* 








Phase V P-V 
p m m 
0'072 6°42 0°47 
0°168 6°57 0°55 
0°265 6°64 0°67 
0'270 6°67 0°67 
0'270 6 «46'66—s(0. 665 
0'272 6:69 0°68 
0'269 6°66 0°67 
0408 6:80 0:'765 
0°438 6°81 0°78 
0°475 6°89 0°80 
0°476 6°88 0°78 
0°476 6:88 079 








410 


JD 


2434872‘001* 
3797 °905 


4826:107* 
4840°058* 


4867 -984* 
4844°985* 
4859°015* 
4872°995* 


3778°070 
4808 °158* 


3796°915 


3779°015 
3765140 


4855°991* 
4833°069* 
3776-010 

4820'022* 


4834°035* 
4847 °992* 


4830°003* 
4844'069* 


SW Velorum 08" 40™-4, — 47° 
Max=JD 2434872-95 +231-473 


JD 


2434872°999* 
4830°007* 
4808 -170* 
4833°072* 
4834°040* 
4840°063* 
4844°988* 
5481°238* 
4847°995* 
4872:°004* 
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T Velorum 8" 34™-4, —47° 00’ 
Max = JD 2426274432 + 419°639769 























Phase V P-V 
Pp m m 
0°'O17 9°71 0°62 
0018 7°67 0°605 
o'018 §=6—67°69~—s«O'61 
O°125 7°81 0°73 
o-132 738 0-793 
0128 7°84 0°73 
Ovrsi 7°84 0°75 
0194 7°93 0°78 
0°218 7°07 0°77 
0°231 7°97 0°83 
0°224 7°97 080 
0°245 7°99 0°785 
0°257 8-02 o'8o 
0'251 8:00 §=-00'79 
0°305 8:10) 085 
0°448 8-15 o'88 
0°455 816 0°875 
0°452 8-16 0°88 
0566 §=8-28 8o-g: 
0°626 = 828 — 0705 
0800 «= 825s (087 
0814 8:09 0°75 
0°334 807, 0°74 
0°842 8:01 0"75 
0°838 §=68:04 0-745 
0965 7°69 0°60 
0'996 7°68 0°63 

02’ 

E 
Phase V P-V 

p m m 

0°002 7°55 0°66 
O’171 7°80 = 0705 
0°240 7°88 1°02 
0°301 7°95 1°14 
0°342 -7'99—Ss«a1"16 
0°599 8:46 1°47 
080g. 870-137 
0'°914 8-70 1°22 
0937, 84912 
0-960 = 809 O'01 


? Carinae 
JD 


2434872:'010* 
4873°030* 
3775 °986 
4557°960* 


3776°988 
4808-116* 


3778°035 


3779003 
4845 °004* 


4527°064* 


3852°859 
4848-016* 


3784°010 
3855-877 
4818-067* 
4820°043* 
4469°187* 
4826-201* 


3797 O15 
4472°166* 


4830°028* 
3765 097 

4867 -989* 
4833 °084* 
4087°153* 
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9g" 42™-5, —62° 03’ 
Max =JD 2427089°82 + 35-5323 E 
Phase 4 P-V 
p m m 
O'OI7) = 33°34 «08 
0°046 3°28 0°89 
o'171 3°44 1°04 
0°179 3°40 1°04 
O°175 3°42 1°04 
0*200 3°44 1°09 
o°219 3°47 1°09 
o°210 3°46 1°O9 
0229 3°49 28S 
0256 ge 1°75 
0°257 3°52 1°18 
0°256 3°52 1°18 
0°309 3°60 1°20 
0°335 3°63 1°29 
0°342 3°63 1°29 
0°338 3°63 1°29 
0°397 3°69 1°32 
0°420 (3°75 1°335 
0°499 3°34 1°35 
O'555 3°91 = 1°38 
0°681 3°99 1°35 
0°722 4°01 1°32 
0°764 4°02 1°305 
0764 3°99 1°36 
0-764 400 1°33 
0°835 4°00 1°29 
0865 3°07 1°25 
0°904 3°84 I'I4 
0°922 3°76 1°08 
0°929 a 1125 
0926 3°76 I‘10 
105 40M-6, —57° 02’ 


VY Carinae 


Max=JD 2434529°90+ 184-9370 E 


JD 


2434833 099* 
4928:072* 





Phase V P-V 
Pp m m 
ool! 6:90 0°74 
0°026 «§=666'89-—so0'71 
0018 )§=6 690 0°75 
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VY Carinae (continued) T Crucis 12% 15™-9, —61° 44’ 
JD Phase V P-V Max=JD 2427578-96 + 64-733178 E 
p m m . . 
4834°136* 0066 7:03 0°86 JD Phase hig 
a din : Pp m in 
4873018" o119 7°06 0°90 2434925'008*  0°022 6°36 0°67 
3851 °g00 0°197 7°29 1°06 
3777°O1S 0'243 7°31 1:07 4845 069* O'150)3=—s-« 6 “41 0°74 
3796°920 0'294 7°44 I'll 4925°949* o°162 6°43 0°76 
3779°025 03497491185 4952°969" 0175 643 0°75 
4557°981* 0483, 7°79 sdT"19 — 


0162 6: "75 
4845°018* 0640 7°90 1°25 I 1°43 0°75 


4560-989* 0°642 7°04 3855891 0°239 6°45 0°79 


) 











4920°971* 0651 7°88 1°21 4873°045* 0°305 6°52 0°83 
— ’ 4953°934" 07318 6°54 0°85 
0°644 7-9! 5°22 4927°017* 0°321 6°52 0°82 
4543°988* 0°744 7°76 1°10 
4848-023* 0°799 7°65 1°03 O°315 6°53 0°83 
4527°076* o851 7:64 0°96 4833°119* 0375 6:57 0°89 
4472'178* 0'952 6°87 0°67 3776°183 0°401 6°56 0875 


a * , , 2 
S Muscae 12! 07™-4, —69° 36’ 4920°992 0°426 «66°63 088 


























- ° * ° 2 “A> . 
Max=JD 2419541-800+91%-658688 E 4544 °O15 0°43 tal —— 
JD Phase V ital 0°432 6°62 o'89 
Pp m m > " 4 ‘ac 
2434833°142* 0169 5:91 0°63 4558:020* o'518 6°70 0°95 
4544°030* 0°237 5°93 0°63 4834°158* 0°529 6°69 0°93 
4921°003* 0'266 5°93 0°70 0'524 6°70 0°04 
4950°000* 0'268 5°97 0°66 
4834°163* 0°275 6-00 0°70 3777186 0°550 6-69 0°935 
2 4848-054* 0°593 6-72 092 
0270 5°97 0°69 4949°992" 0°733, 680 0-92 
- a ; » 3852-950 0°803 6:66 0°84 
4873"066* 0°303 6°03 0°70 4950°967* 0:878 656 0°73 
4959°979 0°370 «614 0°72 4527°096* 0-925 6:48 0-70 
4845 "093° 0°407 O15 vin 4561°003* o-961 6:40 0°66 
4903060 0-408 86614 0°75 4917°999* 0-981 6:36 0°68 
0-408 6°14 wre 0'971 6°38 0°67 
4527°110* 0°485 6°21 0'79 
4952°949* 0°574 6°35 0°84 RCrucis 12" 18™-1, —61° 04’ 
4953°943* 0'676 6°41 0°83 Max JD 2419432-010+ 5 9-825589 E 
3852°943 0°686 640 80082 JD Phase V P-V 
4558°031* 0686 6°45 0°83 p m m 
ssiiittiaaian, ~ “leiapanibis 2434544°007* 0°'072 6°42 0°47 
0°683 6°42 0°83 4527°091* 0168 6°57 0°55 
4848-086* 0-717 6°45 0°79 4848-058* 0°265 6°64 0°67 
3776°118 0°723 6°40 0°78 4917°994* 0°270 6°67 0°67 
— 3776°185 0°270 6°66 0°665 
0°720 0 6©66"42.—s«0"785 4952°963* 0'272 6:69 # 0°68 





4926°043* 0788 6:37 09°75 


3777°095 0°833 6:20 = 067 0'269 ©=— 6°66 0°67 


4927°009* 0888 6:13 0°68 3776°986 0408 6:80 0°765 
4917°985* 0954 5°98 0°59 4953°928* 0438 6°81 0°78 
3855 °888 0'991 5°96 0-60 4925°015*  0°475 6°89 0°80 
4561°019* 0:°996 5°95 0°58 4558°015* 0476 6:88 0°78 








0994 5°96 9°59 0°47 6°88 0°79 








412 


R Crucis (continued) 


JD 


4873°053* 
4825°947* 
4833°129* 
4845 °062* 
4949°984* 


4920°987* 


3825 °955 
4927°020* 


4834°148* 
4950°963* 
3855885 


4561 -008* 
4928°059* 


R Muscae 


JD 


2434952 °952 
4848 -093* 
4833°147* 


4953°947* 
4834°166* 
4527°115* 
3776°198 


4917°977* 
4903 °042* 
4873 °089* 


4558°035* 
4926°040* 


3777°238 


3852°963 
4927°006* 
4544°024* 
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Phase V (P-I *) JD Phase V P-V 
p m m " p m m 
0°555 7°00 0°83 4950°008 0°630 «=6.6'70—s 086 
0635 7:07 0°87 4845°097* 0661 6°74 0°87 
0°702 7°13 0°89 4950°983* 0°760 6-71 0-76 
0'750 7:22 088 4921 008* 0'769 6°66 0°79 
0-761 7:21 0°89 
0°764 6°68 0°775 
0°756 =97°22 (0885 4561°024* 0°836 6:47 0°64 
0°783. 7°16 = 083 §200°154* 0938 6:01 0°43 
0809 7°09)—_(0'76 3855890 0°946 6:06 0°42 
o819 «67°08 -~——io 8 
0°942 6°04 0°425 
o0°814 7708 = =—-00"'785 
0877. 6-88 067 W Virginis 13" 20-9, — 02° 52’ 
0929 6°68 056 Max=JD 2420472'543+ 17%-271671 E 
0'°952. 6°58 os JD Phase V P-V 
0'9900 6: oO: p m m™m 
uae of an 2433063°771 O’OIO 10°30 0°50 
3082-767 O'1IIO =. "9403.85 
2944°976 0°132 9°79 = O40 
—-— Se oe 3083 °762 0168 9°73 0°345 
12 36™-0, —68° 52’ es O°310 = 9°47 = 
Max= JD 2426496-288 + 74-51o211 E oe — =, oe 
— 7" 3087°743 0°398 89 9°42 0°55 
ius V P-V 3089°726 0°513 9°60 38 o'61 
M ua at 3090°729 "57 9°66 0-655 
0°022 5°95 0°43 3091°722 0°629 9°76 0°725 
0060 ©«600-— «049 3040°819 0°682 10°05 0805 
0-070 6:02 0°47 3092°740 0688 I0°1! o0'8o0 
0065 6:01 0°48 0685 10°08 0-80 
O°155 6°13 0°57 3059°740 O-777 10°51 = 
0206 46:19 0-61 as 0°838 10°60 vir 
a te whe 3062°741 O'951 10°49 0°635 
0335 622 0-625 AL Virginis 14% 05™-8, —12° 50’ 
Max=JD 243020260 + 1049-30438 E 
0328 §=66:20) (0-63 JD ‘., "a P-V 
0°365 6°31 o°7I Pp m m 
0°377. 633. 072 2433095 °740 0063 «-9'22—s«O'44 
0388 = 635 0°74 3096-736 0160 9°34 0°54 
a , 3027°855  0°476 99°69 = 0°705 
0°377, 633, -0"725 3089°733.  0°480 9°82 0°69 
0438 650 «= 0-79 
0439 «49648 0°77 0478 9:76 0-70 
: ; 3059°785 0574 9°82 0°565 
— oR 3090°735 0578 9°76 0545 
0°473 648 0-78 
0556 659 3=— 086 0576 =899°79 809555 
0°568 660 0:87 3060°795 0°672 9°76 39048 
0°573, 6620-87 3091°729 = -§ "674 9'69— (0°46 
0°566 6:60 = 0-87 0673. 9°72, (0°47 





sie 
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AL Virginis (continued) X Ser (continued) 
JD Phase V P-V JD Phase V P-V 
p m m p m m 
3081 °743 0705 46965 = 046 4918°045* 0630 4°78 0°79 
2948-063 0°732 9°57 0°42 4925°065* 0°631 4°76 083 


3040°826 °°734 «9°59 = -0"405 











0°630 4°77 o'81 

0°733 «9°58 Og 4869:102* o650 4°75 o'8 
3061°753 0°765 9°55 0°40 4848-115* 0657 4°77 0°81 
3092°750 °°773 «9°36 0°35 

0769 = 9°46 0°38 
3062°750 0°862 = g'21 0°295 
3093 °736 0869 Qg'I2 0°255 











0°654 4°76 o'81 
4603 -051* 0°709 4°82 0°80 
4561°124* 0°730 4°81 o'81 














0866 «916 = 0275 BL Herculis 17" 56™-8, 19° 15’ 
3063°781 0962 = g'I2.— 030 Max =JD 2415196°499 ~ 19-307472 E 
JD Phase V P-V 
RU Camelopardalis 17" 10™-9, +-69° 51’ 
F i 9 p m m 
Max=JD 2423153°88-+22%-1725 E 2433093 °764 ©°450 10°42 0°42 
JD Phase V P-V 4195965 0°452 10°34 0°45 
p m m 
2432887883 o'o12 8-28 0°90 sada -28 Siete 
2822-012 0°042 8-26 0°89 ae orass 
2825-017 0°177 8:26 0-92 3059°951 0°589 10°47 0°46 
2848-969 0°257 8:28 oro 3092809 0°720 «10°56 = 0°37 
seeraes 0°33! 4 die, 3062°944 0878 10°18 0°245 
3009°972 O°519 ° 1°34 4193°934 0°898 10°16 0°24 
2833°031 0°539 «= 8 '83.—s«1°355 a eee 
2944°965 0°587 8-88 1°36 0888 10°17 0°24 
2860°922 0°796 9°05 1°235 
2884-008 0838 876 1-21 4567-908 0926 «6991 = 0165 
2885-899 0'923 8-53 1‘00 
~ x : og1 ‘80 0950 “82 O°105 
2821°017 0°997 8-38 0°925 — . = : eg 


3095 °764 0'980 9°75 o'll 





X Sgr 17% 41™-3, —27° 48 


° 6> “78 . 
Max — JD 2403169°58+-7%-01216 E =e, fe ha 








JD Phase V P-V 3116°757 0036 «6987 orl 
: p m m 4570°828 0°160 ~=10°00 0°285 
irene 0°072 4°20 0°47 4251°812 0165 10°08 O°215 
5012°997 O17! 4°24 0°55 4213°910 0176 10°13 0°24 
4873°109* 0°222 4°32 0°60 aes = 
5308°039* 0°247 4°30 0°61 0167 10°07 0°25 
e * . . ‘6 
4544°149° 0°309 4°43 0°62 3004767 «0'217—«dT002_—S«"25 
4558176 0°309_ 4°43.—ss«O'6 
sors" ow 4°38 ot; 4549°971 0'208 10°03 0°265 
: 7 4574°809 0'205 10°02 0°25 








0°309 4°41 062 


4951°038* 0335 4°43 0°63 
5014°948* 0449 4°51 0°74 
4903°083* 0496 4°54 0°74 3098-792 
4868°254* 05290 4°67 9°75 

4602°038* 0564 467 0°77 

5002°987* 0600 84°73 0°82 3066 °947 0°350 10°30 0°345 
4959°951* 0606 4°74 0°79 3123°728 0°367 10°38 =60-34 


O°210 10°02 0°26 


3073°938 0°286 =10°12 0°275 
0296 10°20 0°28 








o'291 10°16 0°28 








0°603. 4°74 = (0805 0°358 10°34 0°34 
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414 
W Ser 17"58™-6, —29° 35’ 
Max = JD 2403198-30 + 7%-59466 E 
JD Phase V P-V 
p m m 
2434868:258*  0°030 4°25 0°45 
5012'992* 0087 4:23 0°48 
4959°944* O102 4°30 0°49 
4603°055* oO'110 4°31 o'51 
0*106 4°30 0°50 
4869-111* 0°142 4°35 0°55 
4558°190* 0°203 4°47 0°60 
5013°974* 0216 4°48 0°58 
o'210 4°48 0°59 
5014°956* 07346 445 0°64 
4544°144* 07353 447 0°61 
07350 «6. 4446S (00625 
4848-118* 0°378 4°50 0°67 
4925 :060* 0*509 4°78 082 
4918-040* 0°585 4°86 0°84 
4561°138* o'591 4°87 0°85 
0588 4°86 0845 
4903°102* 0618 4°87 0°86 
4873°102* 0668 4°98 °'90 
5002 °980* 0-769 4°92 0°88 
4951°043* 0930 4°62 0°59 
5308'055* 0938 4°59 0°60 
0°934 4°60 0°595 
4602048" 0°977 439 0°53 
4921'150* 0°9904 4°35 0°46 
BB Herculis 18" 41™-3, + 12° 14’ 
Max = JD 2426244:630+74-50720 E 
JD Phase V P-V 
Pp m m 
2433098813 O'Ols 9°74 0°69 
3091 °840 0086 9°76 = 0765 
3122°757 0204 990 0°875 
3093816 = 0349 999 0825 
3116°767 = 0406) 9°95 (0°95 
3094868 0°489 10°11 1°075 
3087889 0°560 =10°28 1'07 
3095°785 O°611 10°30 1°135 
4574885 0°636 =—:10°38 1°09 
0°624 10°34 I‘il 
3096-764 0°742 10°39 I‘1I5 
3119°771 0806 =10°24 1°07 
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«x Pavonis 18 46™-6, —67° 22’ 
Max=JD 2426436°5 +. 94-0738 E 
JD Phase V P-V 
p m m 
2433777295 o7o10 = 406 0375 
4848°128* 0:024 3°98 0°42 
O°01I7 4°02 0°40 
4921°119* 0077 4°07 0°46 
3923 °903 0167 = 4°33. (0'59 
5013°013* 0°195 4°34. 0°68 
5004°002* 0°204 4°37 0°64 
0°200 4°36 0°66 
3852°177 0°262 4°44 -0°735 
5013°984* 0*302 4°52 0:76 
3898-042 0317 4°57) 9°75 
O'310 4°54 90°75 
4959°978* 0°350 4°62 0°78 
4951°0o1* 0°361 4°60 0°79 
3853°167 0°372 «4°69 (0'79 
0361 4°64 °o'79 
4561°166* 0°398 4°71 0°77 
3917°O17 0°408 4°71 080 
5014°962* 0°410 4°71 o'81 
0405 4°71 0°79 
3980°920 0°451 4°71 080 
4925°108* 9 0°507, 4°74 (0°75 
3856-170 0°703 4°57. _-0°485 
3946°913 0°703 4°51 = 050 
5308-021* 0°707 4°49 o'si 
0°704 4°52 0°50 
3928-906 0'718 4°53 0°485 
4918'143* 0°740 4°45 0°47 
4982°004* 0778 4°33 0°40 
3983°917 0-781 = 4'29 (0425 
0'780 )=—s 4431 o'4I 
4910-090* 0852 4°00 0°32 
3948°913 0°924 3°97 0°335 
4920°244* o'971 3°94 0°38 
5002°031* 0985 3°96 0°37 
0978 = 3°95: 9°375 
V496 Aquilae 19% 03™-0, —07° 36’ 
Max =JD 2428719°90+ 6%-8069 E 
JD Phase V P-V 
Pp m m 
2434567844 o7120 7°59 (0965 
74°816 0144 7°60 1*O2 
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V496 Aquilae (continued) AU Pegasi (continued) 
D Ph F P-V 
JD Phase Vo P-V . oe 
p m m 


Pp m m 58-762 O'rrl 9°04 0°69 
75°868 0°299 7°68 1°10 O°113 9°08 0°675 
70°837 0560 §=6—7'86—sar'19 





65-878 0°831 7°99 1°065 O'112 9:06 0°68 

52°962 0°934 7°66 0975 

66-847 0974 7°59 0°99 22°861 0138 =6.9'06—S (068 
51°778 0198 8:96 0°64 
80°701 0260 9°04 0°65 


25°854 0°385 8-99 0°68 


’ 858 0°48 . “765 
AU Pegasi 21" 19™-4, 17° 51’ 54°95 453 g*00 0°705 


em 23°854 0°552 9°06 o'815 
Max=JD 2428729°76 + 29-3978 E 26-844 0-799 «9°23 0845 
JD Phase V P-V 79°705 0°845 9°29 0°825 
p m m 19°872 0892 9°30 0°84 
2433153°893 0-080 9°05 0°73 48-809 0'960 g'22 o'815 


Individual observations of magnitudes and colours of the cepheids are collected 
in Table II where those obtained by Gascoigne and Burr are designated by an 
asterisk. The light and colour curves are plotted in Figs. 1 to 4, where open 
circles indicate observations by Gascoigne and Burr, filled circles those by Eggen, 
and half-filled circles those means in Table II that are formed from observations 
by all three observers; the crosses in the curves for RU Camelopardalis in Fig. 3 
indicate the photoelectric observations made by Lenouvel and Jehoulet (1953). 
The light elements have been taken from the General Catalogue of Variable Stars 
(Kukarkin and Parenago 1948) or, when the results of the present discussion 
indicated that revision was necessary, from Table III. For « Pavonis we have 
used the period of 94-0738 that is indicated by a combination of our 1951 and 
1953-54 observations. 

Tase III 


Revised elements for seven cepheid variables 


d 
W Sgr: max =JD 2403198-42 + 7°594760E 


S Mus: rising median=JD 2434926°89 + 9°6598 E 
SV Per : max =JD 2419055:040+11'12880 E 
VY Car: max JD 2434529°90 +19°9370 E 
RZ Vel: max JD 2434844°96 +-20°39739 E 
SW Vel: max =JD 2434872:95 +23°473 E—3°75.10-5 E* 
RS Pup : max -JD 2434823°23 +41°38535 E 


There is a substantial overlap between the cepheids observed by Eggen and 
those observed by Gascoigne and Burr, from which it can be seen that the 
agreement is satisfactory. Eggen found S Muscae to be o™-o5 brighter at 
minimum than did Gascoigne and Burr and he also observed the post-maximum 
branch of « Pavonis consistently fainter. In the second case at least it seems 
probable that this difference reflects a real variation in the light curve, but we 
have drawn only the mean curve in the figure. 

Improved light. elements are suggested for seven stars in Table III and the 
data relative to four of these are shown graphically in Fig. 5. In the figure, the 
earlier results originated almost entirely from the Harvard College Observatory 


(Gaposchkin and Payne-Gaposchkin 1946); the points for RZ Velorum are 
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. 1.—Light and colour curves for individual cepheids. Filled circles represent observations 
by Eggen, open circles by Gascoigne and Burr, and half-filled circles indicate those means in 
Table II that are formed from observations by all three observers. 
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Fic. 2.—Light and colour curves for individual cepheids. Filled circles represent observations 
by Eggen, open circles by Gascoigne and Burr, and half-filled circles indicate those means in Table II 
that are formed from observations by all three observers. 
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Fic. 3.—Light and colour curves for individual cepheids. Filled circles represent observations 
by Eggen, open circles by Gascoigne and Burr, and half-filled circles indicate those means in Table II 
that are formed from observations by all three observers. 
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Fic. 4.—Light and colour curves for individual cepheids. 
by Eogen, open circles by Gascoigne and Burr, and half-filled circles indicate those means in Table IT 


that are formed from observations by all three observers. 
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from a discussion by Oosterhoff (1936). Our observations of this latter star 
do not support the term in E? in the light elements suggested by Oosterhoff, 
and, indeed, there is little reason to suppose that its period has changed in the 
last 40 years. The change in the period of SW Velorum, however, continues 
to be large and regular as found by Oosterhoff (1936), but we find a somewhat 
smaller secular term than his. The period of RS Puppis has apparently not 
altered since the curious phase change noted by Gaposchkin and Payne-Gaposchkin 
(1946) and shown schematically in Fig. 5. 
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Fic. 5.—Period variations for RZ Velorum, VY Carinae, SW Velorum and RS Puppis. 


For most of the stars in Fig. 5 the period variations—when they exist—are 
small and regular. There are four cepheids, however, for which these variations 
are large and erratic—x Pavonis(Parenago 1946), AL Virginis (Payne Gaposchkin 
1938), WVirginis (Nielsen 1954) and RUCamelopardalis (Lenouvel and 
Jehoulet 1953); all these have been so fully discussed as to merit no more than 
brief mention here. For « Pavonis we can add the following two points to 
Parenago’s (1946) table: 


E=1576, O—-C= +3%-86 (Eggen), 


E=1692, O—C= +1%-98 (Gascoigne and Burr), 


Both W Virginis and AL Virginis are recognized members of population type II 
and evidence will be presented later to indicate that x Pavonis and RU Camelo- 
pardalis may also be members of that population. 
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The 63 observed variables, including those discussed previously (Eggen 1951), 
are listed in Table IV together with the following information: Column (1) name: 
(2) logarithm of the period in days; (3) maximum colour; (4) minimum colour; 
(5) colour amplitude; (6) maximum visual magnitude; (7) minimum visual 
magnitude; (8) light amplitude; (g) a classification according to the shape of 
the light curve, where II indicates a probable member of population type II; 
(10) radial velocity amplitude in kilometres per second; (11) source of the value 
in (10). The magnitudes and colours for the 33 variables given in the previous 
F publication (Eggen 1951) have been re-reduced from the original observations 
to the (P,V),-system; additional observations made since 1951 have been 
secured for most of these stars. Since the magnitudes and colours for 6 Cephei 
obtained with the 12-inch refractor may have been affected by the light of the 
visual companion, the values of V,; and (P-V), in Table IV were obtained by 
conversion from the six-colour observations by Stebbins (1945). 











4 TaB_e IV 
i Photometric data for 63 galactic cepheids 
; (P-V zg Ve 
(1) (2) 3) © (5) (6) (7 8) (9) (10) (11) 
Name log P max min amp max min amp Cl. Vel. Ref. 


BL Her o°116 O10 60°46 «0°36 9°78 10°56 0-78 II 27 3 
SW Tau 0200 0°32 0°72 0°40 9°29 10°06 0°77 II 38 3 


A eb NN a Ph ETS ie? DS 


SU Cas* 0-290 0°49 40°64 O'15 5°75 6°13 0°38 ? 22 3 
AU Peg 0°379 0°64 0°84 0°20 8:98 9°30 0°32 Cc ae 
DT Cyg  0°398 0°36 O51 OTIS 5°62 5°92 0°30 . 17 2 
SZ Tau 0498 0°62 080 018 6°34 6:68 0°34 Cc 22 2 
SS Sct 0°565 0°66 o'92 0°26 7°82 8:30 0°48 B 26 3 
RT Aur 0°572 0°30 600°68 )=—0 38 5°05 5°78 0°73 A 33 3 
SU Cyg = o585 0°30 0°55 0°25 6°45 7:06 o-61 A 47 3 
~UMi_  0°s599 0°45 0°50 0°05 1:96 2°06 or10 c 6 2 
ST Tau* 0-605 0°52: 0°92 0°40: 7°75: 8°50 0°75 ? 38 3 
AH Vel 0°626 0°39 O'50 «~O'll 5°47. 5°84 0°37 Cc 17 6 
Y Lac 0°635 0°44 0°75 O'31 8-77. 9°45 0°68 A 39 
V Vel* 0640 O’°5I 0°77 0°26 7:26 87°80 0°54 ? - 
T Vul 0°647 0°36 60°67 0°31 5°43. 6:09 0°66 A 35 2 
FF Aql 0°650 0°57 o'69 oO712 5°17 5°53 0°36 c 16 2 
T Vel 0°667 061 0°94 0°33 7°68 8:28 0-60 B 39 6 
V350 Ser o'712 0°60 0°93 0°33 6:90 «=. 7°66 0°76 J 39 3 
BG Lac 0°727 0°63 «(0°94 «0°31 8:46 9g:02 0°56 B 34 3 
5Cep  0°730 0°34 0°72 «20°38 3°46 4:28 082 A 39 2 


Y Ser 0°761 0°53 O'gr 0°38 5°37. 6:10 


0°73 A 40 2 

R Cru 0°765 0°44 o'89 0°45 6°40 7°21 O81 A 34 6 
FM Aql 0°786 0°96 1°32 0°36 7°82 8°51 0°69 B 35 3 
X Vul o'8o1 1°03. 1°47 0°44 8:36 gtI5 0°79 A 47 3 
RR Lac 0°808 0°55 096 0°41 8:38 9°24 0°86 A 26 3 
XX Sgr 0808 0-79 «6°17 «00°38 8:30 9:02 0°72 B 57 3 
T Cru 0828 0°67 0°94 0°27 6:36 ©6679. 0°43 A 34 6 

U Ser 0829 0°75 I'I5 0°40 6:25 6:92 0°67 B 38 3 
V496 Aql 0833 0:98 I'IQ O°21 7°59 «©=6.7'86_-«00'27 C 14 6 
U Aql 0846 0°68 1°08 0-40 6:05 6:80 0°75 B 42 2 
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(1) 
Name 
X Sgr 
7 Aql 

BB Her 
R Mus 
W Ser 


RX Cam 
S Sge 
x Pav 

FN Aql 

YZ Ser 


S Mus 
B Dor 
¢ Gem 
AL Vir 
Z Lac 


SV Per 
RX Aur 
TT Aql 
RW Cas 

X Cyg 


Y Oph 
SZ Aql 
W Vir 

VY Car 
RY Sco* 


RZ Vel 
WZ Ser 
RU Cam 
SW Vel 
T Mon 


1 Car 
RS Pup 
SV Vul 
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TABLE IV (continued) 


(2) 


log P 


0846 
0°856 
0876 
0876 
0880 


0898 
0°920 
0°958 
0°977 
0'980 


0985 
0°993 
1:006 
1°O13 
1°037 


046 
065 
138 
‘170 
‘215 


~~ — 


‘234 
‘234 
°237 
278 
308 


— ot 


*310 
°339 
*345 
371 
*431 


~— oo 


1°550 
1°617 
1°654 


(3) 


max 


0°45: 


0°48 
0°69 
0°42 
0°45 


0'90 
o'sI 
0°32 
o'91 
0°74 


0°59 
0°46 
0°59 
0:28 
0°74 


0°67: 


0°69 
o'81 
0°77 
070 


1‘07 
0°88 
0°34 
0°67 
ra: 


0°58 
0°86 
0°90 
0°66 
0°79 


0°89 


1°03: 


0°94 


(P-V )xz 
(4) 
min 


o'81 
o'9I 
I‘12 
0°87 
090 


1°32 
0°93 
080 
1'23 
1°09 


0°84 
087 
0°90 
0°70 
5°27 


‘02 
‘OI 
“46 
‘40 
-36 


~—_ om 


1°36 


~~ ot 
Sab Sat ae 
a | 


1°38 
1°62 
1°66 


(5) 


amp 


0°36: 


0°43 
0°43 
0°45 
0°45 


0°42 
0°42 
0°48 
0°32 
0°35 


0°25 
0°41 
0°31 
0°42 
0°53 


0°35: 


0°32 
0°65 
0°63 
0°66 


0°29 
0°68 
o'sI 
0°55 


0°3: 


082 
0°68 
0°47 
080 
0°57 


0°49 


0°59: 


0°72 


Ve 


(6) 
max 
4°16: 
3°50 
9°72 
5°94 
4°22 


7°20 
5°25 
3°96 
7°95 
6°83 


5°90 
3°39 
3°68 
g°12 
7°87 


8°54: 
7°38 
6°38 
8-72 
5°76 


5°91 
7°84 
9°44 
6°87 
74: 


6°44 
7°26 
8-20 
73S 
5°60 


3°28 
6°42: 
6°64 


(7) 
min 
4°84 
4°27 

10°39 
6°74 
5°00 


8-03 
6:03 
4°73 
8-66 
7°58 


6°42 
4°05 
4°15 
9°79 
8°77 


9°33 
8-06 
7°50 
9°76 
6°82 


6-38 
9°03 
10°60 
7°91 
8-15 


7°62 
8°34 
9°10 
8°75 
6°45 


4°00 
7°59 
7°77 


(8) 


amp 


0°68: 


0°77 
0°67 
0:80 
0°78 


0°83 
0°78 
0°77 
o'7I 
0°75 


0°52 
0°66 
0°47 
0°67 
0-90 


0°79: 


0°68 
I‘I2 
1°04 
1°06 


0°47 
I‘I9 
1°16 
1°04 


0°75: 


1°18 
1°08 
0'90 
1°20 
0°85 


"72 


3°37 : 


1°13 


* Observed only at maximum and minimum. 


(9) 
Cl. 
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(10) (11) 


Vel. 


32 
40 
24 
43 
44 


48 
36 
32 
40 
34 


29 
36 
28 
35 
44 


36 
27 
52 
53 
56 


18 
64 
67 
37 


57 
62 


20 
54 
40 


58 
53 


Ref. 
6 


NW WwW WwW WwWwn AD ww Ar WwW Daw N 


wre 


6 


“SW 


mn a 


The sources given in Column (11) for the radial velocity data in Column (10) are as follows :— 


In addition to RU Camelopardalis mentioned above, observations comparable 
with ours have been published for only one other cepheid in Table IV. Harris 
(1953) has observed the bright variable {Geminorum and the results, when 


SN QDuUt wn & 


. Abt, H. A., 1954, P.A.S.P., 66, 65. 
- Quoted by Eggen, O. J., 1951, Ap. 7., 113, 367. 
. Joy, A. H., 1937, Ap. F., 86, 363. 
. Herbig, G. H. and Moore, J. H., 1952, Ap. #., 116, 348. 
. Sanford, R. F., 1956, Ap. F., 123, 201. 

. Stibbs, D. W. N., 1955, M.N., 115, 363. 
- Jehoulet, D., 1953, C.R., 236, 663. 
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reduced to the (P, V),-system, show the following agreement with ours: 
(P-V) max (P-V) min V max V min 
Harris o™-59 o™-89 3™-67 4™-16 
Present © *59 © “go 3 °68 4°15 


3- Discussion.—The study of the relations between the various observationally 
determined parameters of cepheid variables has received much attention, but 
the subject is complicated and clear-cut results have been difficult to establish. 
The two most important results are those of Hertzsprung (1926), who showed 
that, for cepheid light-curves with humps, the phase of the hump decreases as 
the length of the period increases, and of Parenago and Kukarkin (1936), who 
found that most of the characteristics of the classical cepheids undergo dis- 
continuous changes at a period of about 10 days. The relationships between 
shape of light curve and length of period have been well summarized by 
Payne-Gaposchkin (1954). 

Ina previous paper (Eggen 1951) the cepheids were divided into three groups, 
A, Band C, on the basis of the relationship between light and colour amplitudes 
and periods. From the material then available it was found that this classification, 
in general, divided the stars also into (i) those with 5 Cephei-like or smooth, 
asymmetric light curves—type A; (ii) those with 7 Aquilae-like or asymmetric 
light curves showing a secondary hump—type B; and (iii) those with a nearly 
sinusoidal light curve, similar to ¢Geminorum—type C. In the light of the 
augmented material considered here, we have retained the classification system 
but altered the basic criterion to one of symmetry: if the maximum follows the 
minimum by an interval of more than o”-37, we assign the star to type C; if 
less, to types A or B according to whether humpsare absent or present, respectively. 
The somewhat arbitrary value of o’-37 was suggested by the previous work. 
These criteria are not always free from ambiguity, and in particular there are a 
few stars, of which 8 Doradus and VY Carinae are examples, with humps super- 
imposed on symmetrical light curves. In such cases the property of symmetry 
is taken as the fundamental one. 

The colour and light amplitudes of the variables in Table IV have been plotted 
against the periods in Figs. 6a and 6), respectively. ‘The stars of type A, 
including three doubtful cases labelled type AB in Table IV, are shown by open 
circles, type B by filled circles, type C by crosses and population type II cepheids 
by plus signs; the four stars for which the classification is unknown or uncertain 
because of incomplete data, as indicated in column (9) by a question mark, are 
plotted as half-filled circles. A comparison of Fig. 6a with the corresponding 
Figs. 25 and 26 in the previous paper (Eggen 1951) shows why it was necessary 
to modify the original classification of the types A, B and C on the basis of colour 
amplitude alone. The increased scatter in Fig. 66 is due partly to revision of 
some of the previous data but mostly to newly observed cepheids such as 
RT Aurigae, ST Tauri and X Sagittarii. However, although the criteria have 
been changed, the classifications, except for seven stars, have not. These seven 
stars are: DT Cygni and SZ Tauri, which become type C instead of type B; 
SZ Aquilae, B instead of A; WZ Sagittarii, A instead of B; T Monocerotis 
and SV Vulpeculae, A instead of C; and Z Lacertae, C instead of B. 

Preliminary to further discussion of Figs. 6a and 66 we should point out the 
existence of certain selection effects, The cepheids studied in this paper were 


29* 
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chosen primarily on the basis of apparent magnitude and, considered per unit 
volume of space, the intrinsically brighter cepheids, that is, those of longer 
period, are therefore disproportionally numerous, a bias that has been accen- 
tuated by the deliberate addition of long period cepheids to the programme. 
Also, as discussed below, the five faintest variables in Table IV may belong to 
population type IT. 
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Fic. 6.—Colour amplitudes (a) and light amplitudes (b) plotted against the logarithm of the periods. 
Open circles indicate cepheids of type A, including three doubtful cases indicated in Column 9 
of Table IV as type AB, closed circles type B, crosses type C and plus signs population type II 
cepheids; the four stars indicated by question marks in Column 9 of Table IV are plotted as 
half-filled circles. 

It is apparent that the same broad features are common to both Figs. 6a 
and 6b. For example, omitting for the present the five faintest stars mentioned 
above, (i) the cepheids of type C have generally smaller amplitudes than those 
of types A and B; (ii) type A variables greatly predominate in the period range 
3°5 to 6°5 days, and type C in the range from g to 12 days; (iii) all the type A stars 
in the period range from 3°5 to 6-5 days have greater light amplitudes than those 
of types B and C. 

A selection of light curves illustrating the three types of cepheid is shown in 
Figs. 7, 8 and g. Fig. 8 clearly shows the phenomenon noted by Hertzsprung 
(1926), that is, the relatively smooth progression of the hump from minimum 
light at periods near 3°5 days, up the descending branch of the curve to maximum 
light at 11 days, then down the ascending branch to minimum light again for 
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Fig. 7.—Light curves of typical type A cepheids. 
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Fic. 8.—Light curves of typical type B cepheids. 
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Fic. 9 
The light curves for 8 Doradus and Z Lacertae, 


10 


which show type C light curves when the classification by symmetry is made but also have the hump 
characteristic of type B, are shown for comparison. 


Fic. 9.—Light curves of typical type C cepheids. 


Two type A variables, T Crucis and | Carinae, 


and the population type II cepheid, W Virginis, are included in the figure for comparison. The 


peculiar variable VY Carinae is intermediate between types A and C. 
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periods near 20 days. The only apparent break in this progression is the slightly 
misplaced hump in the light curve of BB Herculis which, since it has the high 
radial velocity of 100 km/sec relative to the stars near the Sun, may be a population 
type II object. The later star is included in Fig. 8 because such marked similarity 
between the light curves of population type I and type II variables is unusual ; 
the light curves of such well-known population type II variables as W Virginis, 
RU Camelopardalis and, to a lesser extent, of « Pavonis, are all quite distinctive. 
Two stars with periods near 10 days that are classified as type C on the basis of 
symmetry but which show the humps characteristic of type B, 8B Doradus and 
Z Lacertae, are included in Fig. 8. 
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Fic. 10.—(a) Light and colour curves of RT Aurigae, ST Tauri and Y Sagittarii. (6) Light 


curves of the population type II cepheids BL Herculis, SW Tauri and variable No. 60 in w 
Centauri. 
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Four additional variables that stand out from the general run in Figs. 6a 
and 66 need mention. The two of shortest period, near 14-5, are BL Herculis 
and SW Tauri, for which both the light and colour amplitudes seem much too 
large. In Fig. 10b we compare their light curves with that found for variable 
No. 60 in w Centauri by Martin (1938); the similarity is striking. Spectroscopic 
evidence has led Joy (1949) to identify SW Tauri (6= — 28°) with stars of similar 
periods in globular clusters and apparently BL Herculis (b= + 18°) is also a 
galactic counterpart of the cluster variables. The systematic radial velocities 
of both variables are small but there is little doubt that the stars belong to 
population type II. The longer period variables RT Aurigae and ST Tauri also 
appear to have abnormally large light and colour amplitudes and Fig. 10a shows 
that their nearly identical light and colour curves match closely those for 
Y Sagittarii which has a fifty per cent longer period. Both RT Aurigae and 
ST Tauri lie within 10° of the galactic plane and except for their anomalously 
large amplitudes, are indistinguishable from the run of classical cepheids. 
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There appears to be little or no correlation between period length and light 
or colour amplitude for the six type C cepheids of shortest period although 
Fig. 9 shows their light curves to be quite similar; two type A variables, T Crucis 
and / Carinae, and the population type II cepheid, W Virginis, are included in 
Fig. 9 for comparison. The peculiar variable VY Carinae is intermediate between 
types A and C. The well-known break in the correlations between various 
observed characteristics of cepheids at periods near 10 days is reflected in Fig. 6.4 
where only « Pavonis and Z Lacertae have light amplitudes large enough to fit 
in with the general run of the A and B types. The radial velocity of « Pavonis 
is not particularly high but the variable has a peculiar light curve (Fig. 2), 
variations in the shape of which are well established (Payne-Gaposchkin 1945), 
and an erratically varying period discussed above; it is probably a population 
type II cepheid. Although the symmetrical light curve of Z Lacertae classifies 
that variable as of type C, the presence of the hump at maximum light may mask 
the true maximum and therefore lead to a wrong classification. 
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Fic. 11.—The correlation between the velocity and light amplitudes. The symbols are those 
used in Fig. 6. 


The observed radial velocity amplitudes of the variables in Table IV are 
listed in column (10) and have been drawn from the sources listed in column (11) 
and plotted against the light amplitudes in Fig. 11, where the symbols are the 
same as those used in Fig. 6. The scatter in Fig. 11 is considerably greater than 
has been found in previous discussions of this correlation (Hellerich 1937 ; Eggen 
1951) and is due mainly to the addition of new material including variables of 
population type II. A particularly clear-cut example of dispersion in this corre- 
lation is afforded by R and T Crucis, both of which have well-established velocity 
amplitudes of 34 km/sec but with visual magnitude ranges of o™-81 and o™-43, 
respectively. Also, T Monocerotis and SV Vulpeculae have very accurately 
determined velocity amplitudes of 54 and 53 km/sec, respectively, but the 
corresponding visual magnitude ranges are o™-85 and 1™13. ‘Two stars that 
depart even more from the general run in Fig. 11 are XX Sagittariiand RR Lacertae, 
but for these the velocity data are less complete. 
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The present photometric data for 17 of the 18 cepheids used by Blaauw and 
Morgan (1954) to establish the zero-point of the period-luminosity relation are 
assembled in Table V. The magnitudes that are used by Blaauw and Morgan 
were obtained from a variety of sources but, presumably, are on the International 


TaBLe V 
Magnitudes of zero-point cepheids (Blaauw and Morgan, 1954) 


Name Myg Px Mpg-PE 
a UMi 2°6 2°48 +o°12 
SU Cas 6°5 6°50 0°00 
SZ Tau 7°3 7°22 +0:08 
B Dor 4°9 4°38 +0°52 
RT Aur 6°0 5°90 +0°10 
¢ Gem 4°7 4°66 +004 

1 Car 5'0 4°78 +0°22 

X Ser 5°4 5°13 +0°27 
W Ser 5°5 5°29 +o°21 
Y Ser 6°5 6:46 +0°04 
U Sgr 7°38 7°54 +0°26 
FF Aql 6°0 5°98 +0'02 
» Aql 5‘0 4°58 +0°42 
S Sge 6°3 6°36 —o'06 
T Vul 63 6°27 +0°03 
DT Cyg 6°5 6°20 +0°30 
5 Cep 4°7 4°40 +0°30 





Mean +0°17 


System; their adopted values are in the second column of the table while those 
determined in the present study are in the third. Our median photographic 
magnitudes, }(Pmax + Pmin), average 0™-17 brighter than those used by Blaauw 
and Morgan and thus their determination of the zero-point requires a correction 
of —o™-17 to reduce it to the International P,-system. 

Acknowledgments.—The participation of the Mt Stromlo observers (S. C. B. G. 
and E. J. B.) in this programme was made possible by the generous loan of a 
multiplier and photometer from the Lick Observatory and the construction and 
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CEPHEID VARIABLES AND GALACTIC ABSORPTION 


S. C. B. Gascoigne and Olin }. Eggen 
(Received 1957 February 9) 


Summary 


From a discussion of the colours and magnitudes of galactic and Magellanic 
Cloud cepheid variables, we have : (1) derived a reddening of +o™-11 for 
8 Cephei and +o™-o8 for a Ursae Minoris from the colours and spectral 
types of their visual companions ; (2) found similar values for the reddening 
of these bright cepheids by spectrophotometric comparisons with F- and 
G-type supergiants; (3) determined, from the proper motions of the brighter 
galactic cepheids, a correction of —1™-7 for the zero point of Shapley’s 
period-luminosity relation; (4) derived the value 17°5+ 1°9 km/sec/kpe for 
the rotation constant A and (5) discussed the large values of the mean absorp- 
tion in the galactic plane that are derived from the assumption that the galactic 
and Magellanic Cloud cepheids have the same intrinsic colour. 





1. Introduction.—Several years ago one of us (Eggen 1951) published 
photoelectric magnitudes and colours for 32 Population I galactic cepheid 
variables. Later, Gascoigne and Kron (1952) and Gascoigne (1955) made 
similar observations for a number of cepheids in the Magellanic Clouds and 
found these stars to be so much bluer than the variables observed by Eggen 
that it was difficult to attribute the same intrinsic colours to both groups without 
invoking an improbably high value for the interstellar absorption. These results 
raised two questions: what are the intrinsic colours of galactic cepheids, and to 
what extent can they be used as distance indicators or, more generally, what 
degree of similarity exists among galactic and extra-galactic cepheids? These 
problems have been briefly considered by Gascoigne and Eggen (1957) and 
the intrinsic colours of the galactic cepheids recently discussed by Stibbs (1956) 
and by Code (1957). 

In this paper we have attempted to (a) determine the reddening of 5 Cephei 
and « Ursae Minoris from the colours and spectral types of their visual 
companions (Section 2), (6) make spectrophotometric comparison between 
the cepheids and some F- and G-type supergiants (Section 3), (c) redetermine 
the zero-point of the period—luminosity relation from the proper motions (Blaauw 
and Morgan, 1954) of the brighter galactic cepheids (Section 4), (d) derive a 
value for the Oort rotation constant, A (Section 5), and (e) study the mean 
absorption in the galactic plane (Section 6). 

The properties important for distance determination are absolute magnitude 
and colour. Our main conclusion can perhaps be best expressed by saying 
that as far as these two properties are concerned, there does not at present seem 
to be any convincing reason for supposing that the cepheids in our galaxy differ 
from those in other galaxies. 

It has been necessary, in this discussion, to ignore possible dispersion in the 
period-luminosity and period-intrinsic colour relations; additional observations 
of cepheids in the Magellanic Clouds are needed before this dispersion can 
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be taken into account. Unless otherwise noted, the magnitudes and colours 
are on the (P, V’),-system (Eggen 1955 a). 

2. Cepheids with common proper motion companions.—The bright variable 
5 Cephei forms ADS 15987 AC with the 6™-5 star BD +57°2547, 40” distant. 
There has been no relative motion in this visual pair since its discovery by Struve 
in 1835. An unpublished velocity of the companion, — 19 + 4 km/sec, has been 
determined from five Mills plates by G. H. Herbig and, since the systematic 
velocity of the variable is — 17 km/sec, the two stars apparently form a physical 
system. Stebbins and Kron (1956) have observed the companion in six-colours, 
and we use their results to determine its spectral type and reddening, following 
the O method described by Johnson and Morgan (1953). Rather than make a 
somewhat uncertain transformation from the six-colours to the (U, B, V)-system, 
we will work with the six-colours directly; six-colour observations are here 
denoted by square brackets. 








(G-1) ’ ' ' 
+0°6 a 
+02 a 
-O2 _ 
-0'6) - 
-1d “ 
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1 j | ] 
fu-p)-'4¢ -1lO -O6 -O2 #02 
Fic. 1.—The colours [G—I] and [U-B] for Og-5 to B1+5, Ia and Ib supergiants. 


The following stars of spectral type Og-5 to B1-5, Ia-Ib, have also been 
observed on the six-colour programme (Stebbins and Whitford 1945; Stebbins 
and Kron 1956): 


Star [U-B] [G-I] Spectrum O 
mn m 

HD 195592 —0°32 0°66 Oog:'s la —o-'76 
¢ Ori —1°66 -1°38 Oog:s5 Ib 0°75 
e Ori 1°64 1°33 Bo Ia 0°76 
69 Cyg —1°59 —1°21 Bo Ib 0°79 
HD 194839 —o0-21 +o'81 Bo-s la 0°75 
26 Cep 1°04 0°55 Bo-s Ib 0°68 
«x Cas 1°34 0°77 Br Ia 0°83 
HD 216411 0°69 o'12 Br la 0'77 
HD 169454 0°27 +o-78 Bi la 0°79 
HD 190919 1°13 0°62 Bi Ib 0°72 
¢ Per -1°28 0:84 Bi Ib 0'72 
HD 194297 O14 0°95 Bi-s Ia 0°77 
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In Fig. 1 we plot [G-J] for these stars against [U-B]. The line on which the 
points lie—the reddening line—is given by the equation: 


[ U-B] =0°662[G-J] —o™-765. 
This enables us to define a quantity 
Q=[U-B] — 0°662[G-J], 


which is independent of the degree of reddening and depends only on the spectral 
type of the star. 

To find the intrinsic colour of a star with spectral type given by Q, we use 
the following list of stars considered by Johnson and Morgan (1953) and by 
Johnson (1955) to be almost or entirely unreddened: 





Star (U-B} [G-l Sp. Q Egr 
7 U Ma —1°27 —1°34 B3 V —o-38 —o'oI 
7 Aur —1°23 —1°'28 B3 V —o'38 +0°05 
7 Her —I1'I4 —1'29 Bs IV —o'29 —0'02 
wu Eri —I1'Io —1'23 Bs IV —o'28 +0:03 
19 Tau —0'95 —1'16 Bo IV —o'18 +0°04 
17 Tau —o'89 —1°13 B6 III —Oo'l4 +0°05 
B Tau —1°00 —1'17 B7 Ill —o'22 +0°06 
a Leo —o'84 —I°Is§ B7 V —o-08 —o'ol 
18 Tau —o'82 —1°l4 B8 V —o'06 —o'ol 
B Cyg C —0o'74 —1°04 B8 V —0'05 +008 
5 Cep B —o'80 —o'98 (B6) —or'ls +0°20 
e | ' ' 
[G-1] e 











! L i 
Q -0O3 -02 -O1 





Fic. 2.—The Q-values for B3 to B8 type stars. 


In Fig. 2 we plot [G—J] for these stars against Q, and, giving greater weight to 
the bluer stars, find for the intrinsic colours, 


[G-I] = 0°620 — 1™-09. 
The colour excess E,,, on the [G—I]-system is then 
Eg;=[G-I] —0°620 + 1™-09. 
The measured colours of § Cep B give Q= —o™-15, corresponding to a spectral 
type of B6 and colour excess of +0™-20. Confirmatory evidence is provided by 
W. P. Bidelman who informs us that the spectrum of the star is similar to that 
of 16 Tauri, which is B7 IV on the MK system. It remains to put the colour 


excess on the (P, V), system. From the data given by Johnson and Morgan 
(1953) we can find the relation between [G-J] and B-V; that between B-V and 
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(P-V), has been given by Eggen (1955 a). Combining the two, we have, for 
reddened B-type stars only, 

(P-V )p=0°444 +.0°55[G—] 
whence the colour excess for 5 Cep B, and hence presumably for 5 Cep itself, 
is +O™rT, 

The important assumption in this procedure is that the companion is 
a normal star and not, for example, a binary with dissimilar components. 
Unfortunately the large absolute magnitude dispersion of the B-type stars, and 
the resulting uncertainty of the spectroscopic luminosity classification of 5 Cep B 
preclude, at present, the use of this star for calibrating the zero-point of the 
period-luminosity relation. 

Polaris also is a wide visual binary, and with an eighth magnitude star 18” 
distant forms ADS1477; the companion shares the proper motion, 0°-046, 
of the variable. The radial velocity of the companion from four Lick plates 
is — 13 km/sec, and from seven Mt Wilson plates it is — 5 km/sec; these figures 
are in only fair agreement with the systematic velocity of the variable, — 16 km/sec 
(Roemer 1955). ‘The magnitude and colour of the companion, determined by 
Kron and Eggen from one night of good seeing with the Crossley reflector, are 
V,=8™45 and (P-V),= +0™-31. Since the median brightness of the variable 
is V,=2™-00, the photoelectric magnitude difference is 6™-45, which agrees 
well with the visual magnitude difference found by the Harvard observers 
(Pickering, 1879) 6™-62 and by Stebbins (1907), 6™-63. 

The most recent determinations of the spectral type of the companion are as 
follows : 

Mt Wilson dFr1 (Adams et al 1935) 
W. W. Morgan F2V_ (1938) 
W. P. Bidelman F3 V_ (unpublished). 


For comparison we have the following F1 V to F3 V stars for which accurate 
colours and absolute magnitudes are known: 


Star (P-V)z My Sp 
37 UMa~ +o-21 +3°1I Fr V (Roman 1949) 
78 UMa~ = +0°235 +-3°14 F2 V (Roman 1949) 
HD 106946 -++-0°23 +3°O1 F2 V (Weaver 1952) 
¢ Ser +0°255 +2°79 F3 V (Johnson and Morgan 1953) 
Means +023 |-3°O1 F2 V 


The absolute magnitudes of the first three stars are found from cluster parallaxes 
and are taken from a previous paper (Eggen 1955 b) the luminosity of the last 
star was determined from the trigonometric parallax of 0-043 +0°005. If we 
adopt a spectral type of F2 V for the companion, the reddening is +o™-08. Also, 
if we adopt + 3:0 as the visual absolute magnitude of the companion, we find 
M,, = —3™5 for the variable at maximum; the corrected distance of Polaris 
is then about 100 parsecs. Since Polaris is a single-lined spectroscopic binary, 
the absolute magnitude of the cepheid component will, of course, be somewhat 
fainter than this. 

3. Spectral types and multi-colour photometry.—One means of estimating the 
intrinsic colour of a cepheid is to compare the observed colour at some phase, 
say maximum light, with the colour of a non-variable, and supposedly unreddened, 
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star of the same spectral type. Code (1947), who has classified the spectra of 
18 cepheids on the Yerkes system, finds that at maximum light they all have 
nearly the same spectral type, F61b, and that the general features of cepheid 
spectra match those of standard supergiants well, although the hydrogen lines 
are abnormally strong at maximum and the CN band at A 4215 is abnormally 
weak at minimum light. 

The only reliable estimate of the intrinsic colour of an F-type non-variable 
supergiant is that for « Persei, F5 1b, which is a member of the moving cluster 
in Perseus. It was previously shown (Eggen 1955 b) that the B-type stars in 
that cluster are reddened by +0™-10 on the (P— V),-system and this result has 
recently been confirmed by Harris (1956). Since the observed colour of « Persei 
is (P—V),= +0™-36, an intrinsic colour of +o™-26 will be assumed in the 
following discussion. 

The following colour observations of non-variable supergiants have been 
taken from Stebbins and Whitford (1945) and Stebbins and Kron (1956); those 
of » Aquilae are from Stebbins, Kron and Smith (1952) and for 5 Cephei from 
Stebbins (1945). 


Star Sp [B-G] [R-I] 
m m 

a Per Fs Ib —0'0g 0°09 
35 Cyg Fs Ib o'o! +-0°06 
44 Cyg Fs lab +-0°27 +-O°42 
HR 690 F7 Ib o'12 +O'2I 
45 Dra F7 Ib ~0'03 0:06 
7 Aql (Max) F6 Ib -0°02 +O°OI 
5 Cep (oP-o05) F6 Ib —0:08 -0°09 
B Aqr G2 Ib LFOrrs +008 
a’ Cap G3 Ib +-O°24 +0:20 
¢ Cap G4 Ib + O17 + 0°04 
9 Peg Gs Ib 0°23 +023 
7 Aql (Min) G4 Ib +o°22 +or1g 
5 Cep (Min) G2 Ib Lorr2 + 0°09 


The ultra-violet, which serves as such a useful discriminant between tem- 
perature and space reddening among the early type stars, loses its effectiveness 
for the later spectral types (Morgan, Harris and Johnson 1953) and _ has, 
therefore, not been included in this discussion. The values of [B-G] and 
[R-I] are plotted in Fig. 3 where the filled and open circles represent stars 
of types F5—F7 Ib and G2-Gs Ib, respectively. The continuous line is assumed 
to be the reddening path for the F5-F7Ib stars; the broken line that of the 
G2-Gs5 Ib stars. The colour of 5 Cephei near maximum light, indicated in the 
figure by a filled triangle, is almost identical with that of « Persei and the assump- 
tion that the + o™-10 reddening of the latter also applies to the variable is consistent 
with the reddening of +o™-11 found above for 5 Cephei (B). Also, the colour 
of » Aquilae at maximum light, indicated in the figure by a filled square, is 
apparently reddened by + 0™-06 in [B-G] relative to5 Cepheiand« Persei. From 
acomparison of the colours of the F5—F7 Ib stars we find (P-V),, = 1-5 [B-G] + o™-5 
and, therefore, a reddening of + 0™-10 +0™-09 = +0™-19 for 7 Aquilae. 

The argument in the case of a Ursae Minoris, indicated in the figure by a 
cross, is not as straightforward as that given above for the other two variables 
because (a) the spectral type, F8Ib, is intermediate between the two groups 
of non-variable supergiants plotted in the figure as open and filled circles and 
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(6) the known spectroscopic companion may influence the observed colours. 
In any case, the colour of « Ursae Minoris is slightly bluer than that of « Persei, 
which is consistent with the reddening value of + o™-08 for the former, obtained 
above from the visual companion. 

4. Intrinsic colours and absolute magnitudes.—T he available colours for galactic 
cepheids have been given in a previous paper (Eggen, Gascoigne and Burr, 1957). 
The colours of the observed variables in the Magellanic Clouds (Gascoigne and 
Kron 1952, and Gascoigne 1955), which have now been put on the (P, V), 
system, are as follows: 





Small Cloud Large Cloud 
Var. Log P (P-V) max Var. Log P (P-V) max 
m 

HV 1342 1°254 +o°28 HV 886 1°379 0°35 

HV 817 1°276 0:26 HV 1003 1°387 +0°35 

HV 823 1°504 0°40 HV go2 1°421 +O'25 

HV 824 1818 1+-0°39 HV 1002 1°483 0°26 
HV 2294 1°563 +031 
HV 953 1‘680 +O'44 

{B-G) 1 ; on T | T 1 


+02 


fone) 











A. J | 4 ! 1 
{R-1] 0-0 +02 +04 





Fic. 3.—The colours [B-G] and [R-I] for stars of types Fs5-F7 Ib (filled circles), 
G2-Gs Ib (open circles), 8 Cephei at maximum (filled triangle) and minimum (open triangle), 
7 Aquilae at maximum (filled square) and minimum (open square), and « Ursae Minoris (cross). 
The continuous line is the reddening path for the F5-F7 Ib stars and the broken line that of the 
G2-Gs Ib stars. 


The maximum colours of both the galactic and Cloud cepheids have been 
plotted in Fig. 4, where the crosses and open squares indicate respectively those 
of the galactic variables of population types I and II, and the open and filled 
circles those of the Small and Large Magellanic Clouds. The arrows attached 
to the colours of three galactic cepheids « Ursae Minoris, 5 Cephei, and 7 Aquilae, 
indicate the amount of reddening derived for these variables in Section 2 and 
Section 3. 

Fig. 4 brings us to the central question of the paper: do the galactic cepheids 
have the same intrinsic colours as those in the Clouds? The following points 
are relevant: 
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(i) With the exception of the population II variable W Virginis, all the Cloud 
cepheids are bluer than all the galactic cepheids of comparable period. In 
this connection, however, it should be remembered that (a) Cloud cepheids with 
periods less than 17 days have not been observed and (6) the longer period galactic 
variables, being more luminous, will in general be more distant and therefore more 
heavily reddened than those of shorter period (cf. Stibbs 1955 b). ‘The marked 
difference between the colours of the observed galactic and Cloud cepheids is 
emphasized by the colour of the galactic variable 8 Doradus which, although 
only 6° from the centre of the Large Cloud and in the same galactic latitude, is 
decidedly redder than the Cloud cepheids. However, the absorbing material in 
the galaxy is known to be very patchy, and perhaps too much emphasis should 
not be placed on the colour of this one star. 

(ii) The maximum intrinsic colours derived in Section 2 and Section 3 for 
three galactic cepheids are similar to these observed for the Cloud cepheids. 
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Fic. 4.—Observed maximum colours of populations I and II galactic cepheids, crosses and 


open squares, respectively, and those of the Small and Large Magellanic Cloud, open and filled circles, 
respectively. 


(iii) The Cloud cepheids may themselves be reddened. Appreciable 
reddening seems unlikely in the direction of the Small Cloud, but unpublished 
colours, determined at Mt Stromlo by Gascoigne and Rodgers, indicate that some 
of the Large Cloud B-type stars are reddened +0™-1 or more by patchy 
obscuration either between us and the Cloud or within the Cloud itself. On the 
other hand, the existence of stars in both Clouds bluer than any so far observed 
in the galaxy (Code and Houck 1956) indicates that the Cloud cepheids are bluer 
than their galactic counterparts only because they are less reddened. 

(iv) Although the observed sample is smaller, it may be significant that the 
dispersion in maximum colour of the Cloud cepheids is much smaller than that 
for the galactic variables of population I. 

(v) The population I galactic cepheids do not themselves form a completely 
homogeneous group, but differ appreciably in almost all their observed character- 
istics, including the shape of light curve and the light, colour and velocity 
amplitudes. There is, moreover, the possibility that the intrinsic maximum 
colours may not be precisely the same for the three types A, B, and C, 
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(vi) There are a number of results, some well established, which tend to 
differentiate between the cepheids in the Small Cloud on one hand and those in 
the Large Cloud and in the galaxy on the other. Notable examples are the 
dissimilar correlations between period and light amplitude and the very marked 
differences in period distribution (cf. Gascoigne 1955). 

(vii) ‘The distinction between variables of populations I and II is neither 
clear-cut nor, in practice, easy to apply to individual cases. Of particular interest 
in this connection is the relation between such population II galactic variables 
as BL Herculis (14-3) and SW Tauri (14-6) and the variables with periods between 
one and two days that are so numerous in the Small Cloud (cf. Eggen, 1951). 

(viii) ‘The blueness of the galactic population II and of the Cloud variables, 
which lie, respectively, far above and below the galactic plane, and the redness 
of the population I galactic cepheids, which are closely clustered near the plane, 
may indicate that the stratum of absorbing material near the Sun is very patchy 
or very thin. 

This evidence, though not conclusive as it stands, encourages us to investigate 
in detail the consequences of assuming an intrinsic colour at maximum of 0™-25 
for both galactic and Cloud cepheids. There are some spectroscopic (Code 1947) 
and photometric indications of a dependence of maximum colour on period, but 
the effect seems small, and in what follows will be ignored. It has no effect on 
our main argument (cf. Section 6). 

Support for these intrinsic colours comes from measures of polarization. 
Hiltner (1951) has published the following figures for seven cepheids observed 
by us (those for BG Lac were communicated privately): 


Variable l b © (per cent) E, 
BG Lac 62°9 9'7 1°8 O'19 
Y Lac 66:5 4°3 1*2 0°09 
© Cep 72°9 O°4 o'6 0°04 
RR Lac 73°4 2°1 2°6 O'15 
Z Lac 73°5 1°7 2°3 0°24 
RW Cas 96'S 4°1 2°8 0°26 
SU Cas IOI‘O g'! 1°9 O°12 


The listed values of E, were derived from the values of Apg in ‘Table II with the 
assumption that A, =0-75 Apg and Ay=6-1 E, (Morgan, Harris and Johnson, 
1953). The percentage polarizations 6 (per cent) and the colour excesses EF, 
are plotted in Fig. 5, together with those for the B-type stars in the same galactic 
longitudes; the values of E, for the B-type stars were computed from the colours 
observed by Stebbins, Huffer and Whitford (1940) and the intrinsic colours 
derived by Morgan, Harris and Johnson (1953). ‘This figure hardly suggests 
that the reddening of the cepheids has been overestimated. Although the 
material is too meagre for definite conclusions, the apparently better correlation 
between reddening and percentage polarization for the cepheids than for the 
B-type stars certainly deserves further investigation. 

Blaauw and Morgan (1954) have recently redetermined the zero-point of the 
period-luminosity relation using accurate proper motions of 18 bright galactic 
cepheids. Accurate photometry is now also available for all these variables 
except W Gem and, after allowing for the absorption derived from the assumed 
intrinsic colour of +o™-25, we have repeated the computations; the details are 
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Fic. 5.—Percentage polarizations, 8 (per cent), and colour excesses, E\ of cepheids (open circles) 


and B-type stars (crosses). 


listed in Table I. Assuming that the slope of the period—luminosity relation is 
as indicated by the cepheids in the Magellanic Clouds (Shapley 1940), the median 
absolute photographic magnitude, M, is given by the following expression: 


é 


M=C-— 1-74 log P. 


For convenience we will replace C by (C’—2™-0) and use the proper motions 


to determine C’. 


Star 


x UMi 
SU Cas 
SZ Tau 

B Dor 
RT Aur 


¢ Gem 

1 Car 
X Sgr 
W Ser 
Y Ser 


U Ser 
FF Aql 
n Aql 
S Sge 
T Vul 


DT Cyg 
5 Cep 





v-components of proper motion for galactic cepheids ; 


m 
48 
"50 
*22 


4°38 
S 


SIGN 


TABLE I 


Proper motions of bright cepheids 


Moorr Tkpe v 
1°68 0088 557 
5°54 0°405 99 
5°74 0°526 2 
3°54 0°284 17 
5°70 0°550 104 
3°30 0°257 23 
2°22 0°242 55 
4°33 0°363 105 
4°48 0°400 17 
5°34 0°540 r 155 
5°54 0°625 116 
4°70 0°368 26 
3°66 0269 +- 333 
5°32 0°608 53 
5°84 0°622 38 
5°76 0°490 9 
4°04 0°290 - 167 


‘ 


v 
+ 163 
134 
40 

16 
190 


The second column of Table I contains the observed median 


unit is 0”-ooo! p.a. 
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magnitudes, which are corrected in the third column for absorption as discussed 
above. ‘The fourth column contains the computed distances in kiloparsecs with 
C’=o0. The v-components of the cepheid motions (Blaauw and Morgan 1954) 
are listed in column 5 and, in column 6, are reduced to a mean distance of 
300 parsecs. ‘These reduced motions yield a mean parallax of 0”-0328 as compared 
with 0”-00333 derived from the distances in the sixth column of the table. The 
resulting value of C’ is then 5log(0’:00328/0"-00333) or —o™-o4, and the 
period—luminosity relation then takes the following form: 


M= —2™o- 1°74 log P, 


which represents a correction of —1™-7 to the zero-point adopted by Shapley 
(1940). ‘The distances of the galactic cepheids in Table I have been computed 
from this relation. Since Shapley also found that, for the observed median 
photographic magnitudes in the Magellanic Clouds, 


M= 17™-04 — 1°74 log P, 


the distance modulus of the Clouds is then 19™-0, if the magnitude systems are 
comparable (Gascoigne 1955). 

5. Distance and galactic rotation—We can now compute distances for the 
population | cepheids for which magnitudes and colours are known, and this 
is done in Table II. The photometric data are taken from Table IV of Eggen, 
Gascoigne and Burr (1957), and the colour excesses EF converted to photographic 
absorptions Apg by the result Apg=4E (Blanco 1956). ‘These are listed in 
column 4. In the other columns are: the longitude measured from the 
galactic centre, /+32°; the galactic latitude, 6; the heliocentric, rectangular 
co-ordinates X, Y, and Z; the distance, R, of the variable from the galactic 
centre; and the radial velocity, v,, corrected for a solar motion of 20-4 km/sec 
towards /=17°:1, b=19°°6 (Weaver 1955). ‘The X-axis is directed towards the 
galactic centre, /= 328° and b=o", the Y-axis towards /= 58° and b=o° and the 
Z-axis towards the galactic pole. ‘The Sun is assumed to be 8-4 kpc from the 
centre of the galaxy. 

We have used the distances and radial velocities in Table II to discuss the 
differential galactic rotation in the manner first described by Camm (1938). 
If the circular velocity at a distance R from the galactic centre is V, then 


y,/sin (1+ 32°) cosb=RQ(Vo/Ro—V/R) 


_ Rk adv 
A=-2. 3(z)- 


No account is taken here of a possible K-term or of peculiar motions, that is, 
all the orbits are assumed to be circular. We also ignore the higher derivatives 
of V with respect to R. The 37 cepheids for which both radial velocity and 
accurate photometric data are available and for which |sin(/+32°)cosb|>0-3 
are plotted in Fig. 6. V Velorum, the radial velocity of which departs from the 
general run of values in the figure by nearly 50 km/sec, was omitted from the 
discussion and is plotted as a cross. Whilst most of the scatter in the figure 
comes from dispersion in the radial velocities, errors in distance determination 
cannot be disregarded, and give rise to a regression effect which we take 
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and the Oort constant is 
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Distance, absorption, galactic coordinate 


Star 


X Sgr 
W Ser 
WZ Ser 
Y Ser 
U Ser 
V350 Ser 
XX Sgr 
YZ Ser 
Y Oph 
SS Scu 
V496 Aql 
U Aql 
SZ Aql 
TT Aql 
FN Aql 
7) Aql 
FM Aql 
FF Aql 
S Sge 
X Vul 
SV Vul 
SU Cyg 
T Vul 
DT Cyg 
X Cyg 
BG Lac 
Y Lac 
5 Cep 
RR Lac 
Z Lac 
x U Mi 
RW Cas 
SU Cas 
RX Cam 
SV Per 
RX Aur 
SZ Tau 
RT Aur 
ST Tau 
¢ Gem 
T’ Mon 
RS Pup 
AH Vel 
RZ Vel 
T Vel 
SW Vel 
B Dor 
V Vel 
/ Car 
V7 Cr 
S Mus 
T Cru 
R Cru 
R Mus 
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Taste II 

s and radial velocity for 55 population 1 galacti 
cepheids 

Tipe X 4 Z R 
0°34 + 0°34 0°04 Oo! 8-06 
0°40 0°40 O'o! 0°03 8-00 
I‘51 1°47 0°31 0°08 6°94 
0°54 0°53 O'12 0°03 7°87 
0°63 o'61 O'14 0°07 7°79 
0°95 o'9g2 0°21 0°16 7°48 
1°51 1°46 0°35 0°09 6°95 
0°94 089 0°28 O'15 7°51 
0°44 + O'4! O'ls 0:06 7°99 
1°09 o'99 +0°46 0-06 7°42 
0°74 0°65 0°36 O'll 7°76 
0°65 0°54 0°32 O'ls 7°37 
1°81 1°48 1°04 O'l2 7°00 
0'92 rO°74 +0°54 0°07 7°68 
1*22 0°95 +0°76 —0°09 7°49 
0'27 o'1Ig +0'16 0°06 8-21 
o'92 + 0°66 0°63 O'ol 7°76 
0°37 0°24 +0°28 0°03 8:16 
0°61 0°35 0°49 0°07 8:06 
I‘13 0°50 +1°OI 0°04 7°96 
1°33 0°59 I‘Ig 0°02 7°90 
1°00 0°43 0°90 0°02 8-02 
0°62 +Or1g 0°59 —O'lI2 8:23 
0°49 o'll 0°47 o'10 8°30 
0°84 0°20 +o0'81 0°08 8-24 
1°79 0°09 +1°77 —0°30 3:68 
2°56 0°38 232. -@19 913 
0°29 0°08 --0'28 0-00 8:48 
2°26 0°61 2°17 —0°09 9°27 
ryt 0°46 1°65 0°05 9°01 
o'09g 0°04 0'07 0°04 3°47 
2°86 1°78 2°22 —O'IQ 10°42 
0°40 0°27 0°29 0°06 8°67 
0°86 0°70 0°48 0°08 Ql! 
2°21 2°10 0°67 O'OI 10°52 
1°36 1°32 +0°34 0°00 9°73 
0°53 O'si 0-00 Or15 S91 
0°55 0°55 —0-02 Ol! 8°95 
1°45 1°40 —0°32 O°17 9°80 
0°26 0°25 0°02 +0°'06 8°65 
0°77 0°70 —0O°3I1 —o0'02 Q'Il 
I‘OI 0°33 1°02 +0°01 8°79 
0°53 0°08 —0'54 —0'06 8-52 
1°57 0°20 1°56 —0'03 8°75 
1°25 O'll 1°34 —O°O4 8-58 
2°46 -0o'18 —2°45 —o-or 8:93 
0°28 0°00 —0'24 —o'l4 8-40 
1°12 -0'13 she 0°07 8-60 
0°24 + 0°06 0°24 0°03 8°33 
1°50 +O°4I —1°44 0°04 8-12 
0°70 +0°34 —0'61 —o-0g 8:08 
0°67 r0°32 —o'58 o700 811 
1°02 +0"§2 —o'89 -+o-01 7°93 
o'92 +9045 —0O°73 —O'11 7°98 
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into account by Hertzprung’s method (1922). The regression of R on 
v,/sin (/+32°)cosb gives A=22-1, and that of the v, term on R gives 16:5. 
Assuming the errors in v, to be about twice those in R (this is a minimum figure) 
we have the final result A = 17-5 + 1-9(s.e.) km/sec/kpc. 

In addition to the internal error of this result, there may bea further, systematic, 
error arising from the uncertainty of the zero-point of the period—luminosity 
relation; the dispersion in the proper motions alone (Section 4) introduces an 
uncertainty in A of 10 per cent. An independent determination of the zero-point 
would of course eliminate this uncertainty. 

Stibbs (1956) has discussed in some detail the dependence of A, determined 
in this way, on the accuracy of the assumptions and of the observations involved. 
We may add that if the calculation is repeated with a value of 4°75 for the 
Apg/E ratio instead of the 4-0 used here, the constant C in the period—luminosity 
relation becomes —2"™-2 if the agreement with the mean parallaxes derived 
from the proper motions is to be preserved, and the resulting value of A is 
increased to 19-7 km/sec/kpce. 








‘dh sina canbe tdibisalaiamen 
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Fic. 6.—Values of vy/sin (1+ 32°), and the distance, R, of 37 galactic cephetds. 


It is worth emphasizing that the value of A derived in this way is relatively 
insensitive to reasonable changes in the assumed intrinsic colours of the cepheids. 
For example, a change in the assumed intrinsic colour from +0™25 to +0™35 
reduces the total absorption by +0™-4 for each variable and alters the constant, 
C, in the period-luminosity relation from —2™o to —1™-6 but leaves the 
distances, and hence the derived value of A, unchanged. Ina previous computa- 
tion (Gascoigne and Eggen, 1956) in which we assumed that the intrinsic colour 
was given by the expression +0™-10+0"™-13 log P, the derived value of A was 
18-4 km/sec/kpc, or only five per cent greater than the value derived above. This 
follows from the fact that the proper motions determine distances, not luminosities, 
and that any assumptions on luminosities and absorptions are, for the present 
purpose, only a means of extrapolating the distance scale from the nearer to the 
more distant cepheids. 

Values of the rotation constant A have also been recently derived by Stibbs 
(1956), Weaver (1955) and Petrie, Cuttle and Andrews (1956). Stibbs found 
A=19:5 or 19:9km/sec/kpc, depending on particular assumptions for the 
K-term and for the longitude of the galactic centre. Our discussion differs from 
his in that it is confined to cepheids for which only photoelectric light-curves and 
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hence only measured colour excesses are available, and in the explicit derivation 
of the zero-point of the period—luminosity relation. It is the last point which 
enables us to establish the directness of the connection between the proper 
motions and Oort’s A. Petrie, Cuttle and Andrews, using B-type stars with 
measured colour excesses and new spectroscopic absolute magnitudes, found 
17°7 km/sec/kpc. Weaver, also using the cepheids, found A = 10°8 km/sec/kpe, 
but a comparison between his distances and ours suggests that he may have 
underestimated the effect of absorption in the direction of the galactic centre 
and overestimated it in the direction of the anticentre. 

6. Galactic absorption. —In a medium with large-scale variations of absorption, 
values of the “mean absorption” may have little significance, especially over 
large distances, but they do offer a convenient way of expressing relative amounts 
of photographic absorption derived from various assumptions and for that purpose 
will be used here. ‘The mean absorption for the cepheids in Table I is 
+1°58 mag/kpc. The variables with periods respectively less than and greater 
than eight days yield absorptions of + 1-44 and + 1-73 mag/kpc; for the 18 stars 
in the quadrant 13°=>/>289° the absorption is + 2°32 mag/kpc. ‘These figures 
are appreciably smaller than the 3 mag/kpc previously derived (Gascoigne and 
Kron 1952), for three reasons: (a) the cepheids are here assumed to be o''-4 
brighter than previously, (b) we have used the ratio Apg/E=4-0 instead of 4:75 
and (c) the mean has been computed in a different way. If the absorption and 
distance of the ith cepheid are A; and r; respectively, the mean used here is 
2A,/=r; instead of, as previously, X(A,/r;), which gives too much weight to the 
nearer stars. 

The fact that the longer period cepheids show the larger mean absorption 
suggests that their colour excesses may have been slightly overestimated relative 
to the variables with shorter period. That the effect is small is shown by a 
previous computation (Gascoigne and Eggen 1956), in which the assumption 
that the intrinsic maximum colours are given by +0™-10+0-13log P led to 
mean absorptions of +1™-89 and +1™-66 for the cepheids with periods 
respectively less than and greater than eight days. 

These mean absorptions are a good deal higher than those found by Joy 
and earlier workers, but larger values have come to be accepted in recent years, 
especially near the galactic plane, and our 1-6mag/kpc now seems reasonable 
enough. For the nearer cepheids however the situation is different. ‘The 
relevant data for the ten within 400 parsecs are listed in Table III]. The mean 
absorption for the group is 1-5 mag for 400 parsecs, which is a rate of 3:5 mag/kpce. 
Such heavy absorption for the nearer cepheids implies that either they are redder 
than we have assumed, and especially that they are redder than the cepheids in 
the Magellanic Clouds, or that heavily obscuring clouds occur in the immediate 
vicinity of the Sun. Similar results were obtained from an analysis of the colours 
of the B-type stars (Oort 1938) which indicate an average photographic absorption 
of + 2-4 mag/kpc for distances up to 800 parsecs and + 2-8 mag/kpc for the region 
within 300 parsecs if Apg/E,=8-1 (Morgan, Harris and Johnson 1953). 

At first sight absorptions of this order seem quite at variance with those 
determined by the statistical astronomers. To examine the point more closely 
we compare some of our results with those obtained at the Warner and Swasey 
Observatory and recently summarized by McCuskey (1956); these investigations 
were chosen because of their extent and uniformity of treatment. The following 
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Tasie III 


Distances and absorptions for the nearer cepheids 
Star l b 4 


aie "kpe 

° ° m 
» Aql 8-8 —14°5 0°92 0°27 
FF Aql 16°8 5'1 1°28 0°37 
6 Cep 72°9 o'4 0°36 0°29 
a UMi 90°3 26°8 080 0°09 
SU Cas 101°O g'I 080 0°40 
¢ Gem 163°4 13°4 1°36 0-26 
B Dor 238°5 32°3 0°84 0°28 
/ Car 250°7 — 6:8 2°56 0°24 
X Sgr 328°9 1°2 o'8o 0°36 
W Ser 329°2 — 5°4 0'80 0°40 


oe 


five cepheids in Table II lie within the ‘“‘ LF 
and his co-workers: 


fields” discussed by McCuskey 


Star l 5 Anes "kpe Reference 
° bi 
SU Cyg 32°4 1°5 0°20 1*0O LF2 Ap. F., 106, 1, 1947 
RW Cas 96°8 4'1 2°08 2°86 LFs Ap. F., 121, 643, 1955 
RX Cam 113°5 +56 2°60 086 LF6 Ap. F., 115, 479, 1954 
SV Per 130°3 03 1°68 2°21 LF7 Ap. F., 109, 414, 1949 
RX Aur 133°5 0'0 1°76 1°36 LF7 Ibid 


McCuskey’s conclusions are that in LF2 “the absorption...appears to be 
negligible to a distance of 700 pc”. For LF5 the absorption “‘sets in very near 
the Sun, within 100 parsecs at least, and increases rapidly to a value of 2 magnitudes 
at 2000 parsecs.... Thereafter the rate of increase is less rapid”. For LF6 
‘the region is covered by an obscuring cloud which sets in at a distance of 
200 parsecs and extends to 500 parsecs from the Sun.... The total absorption 
of this cloud is 2:0 mag at 500 parsecs (and) thereafter increases uniformly to 
3 magnitudes at 2500 parsecs”’. Heeschen (1951) finds for this later region 
absorptions of about +1™-o at 200 to 300 parsecs and +2™-3 at 800 parsecs. 
In LF7 the data “‘...are consistent with an absorption of 1-2 mag/kpc beginning 
200 parsecs from the Sun”, which would indicate absorptions of +2™-4 and 
+1™-4 for SV Persei and RX Aurigae respectively. Considering the nature of 
the problemand the necessity for averaging over large volumes of space, the general 
agreement is surprisingly good. 


Commonwealth Observatory, Royal Greenwich Observatory, 
Mount Stromlo, Herstmonceux Castle, 
Canberra : Hailsham, Sussex. 


1957 February. 
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THE PHOTOELECTRIC LIGHT CURVE OF V 499 SCO 
G. G. Cillié and E. M. Lindsay 
(Received 1957 March 8) 
Summary 
With the 60-inch Rockefeller reflector of the Boyden Observatory 
photoelectric light curves have been obtained in blue and yellow light for 
the star V 499 Sco. The star is a Beta Lyrae type with two fairly flat-topped 
equal maxima. The light elements are 
Min= JD 2428340%-405 ++ 24-3332077E. 

The eclipsing variable V 499 Sco, which is the same as HD 158155 
(spectrum Bs) and CPD — 32°4551, was studied by Hertzsprung (1), who gave 
the following light elements for the system 

Min = 24283399229 + 24-33328E. 
The star was stated to be either of the B-Lyrae or W Ursae Majoris type with 
equal minima, its magnitude varying between 8™-o and 8™-4. 

V 499 Sco was observed by us on a number of nights during 1952, 1953 and 
1954 with the Linell-King photoelectric photometer with 1P21 tube attached 
to the Rockefeller reflector of the Boyden Observatory near Bloemfontein, 
South Africa. Altogether 1596 comparisons—862 in blue and 734 in yellow 
light—were made with a nearby star HD 157750 (spectrum G5) = CPD — 32°4535 
(Ptm 8™-4, Ptg 8™-8). We used a third nearby star to check the constancy of our 
comparison star. On 3 nights we made 22 comparisons (11 in blue and 11 in 
yellow light) between the two comparison stars. The mean deviation from 
constancy for a single observation we found to be +0™-o21 in blue light and 
+0™-923 in yellow light. We have taken this to indicate at least a very fair 
degree of constancy of our comparison star. 

When plotting our observations according to Hertzsprung’s light elements 
we found that his secondary minimum is 0™-088 deeper in blue light and 0-045 
deeper in yellow light than his primary minimum. We have therefore taken his 
secondary minimum as the true primary minimum. According to our observa- 
tions the phase difference from true secondary to true primary minimum is 
0504. Hence we added 0-504 x 24-3333 to Hertzsprung’s fundamental epoch, 
which gives the revised value JD 24283409405. During 1953 we observed 
3 primary minima of V 499 Sco whose Julian dates are given in Table I. 


TABLE | 
Date JD of observed Wt JD of calculated O-C 
(1953) primary minimum ; primary minimum (in days) 
July 4/5 2434563 °3100 I 2434563 °3100 01-0000 
July 18/19 2434577 3047 I 24345773098 of-0051 
August 8/9 2434598 °3105 5 2434598°3004 +ot-oorl 


Using the above revised fundamental epoch we find that the period which gives 
the best mean square fit is 24-3332977. We therefore propose the following light 
elements for the system V 499 Sco 

Min = 2428340": 405 + 29-3332977E. 
These elements refer to what we consider as the true primary minimum. All 
the phases given in our subsequent tables and light curves were computed with 
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these elements. We should like to point out that the observation of an eclipsing 
variable whose period is so close to 24 days as that of V 499 Sco, presents a peculiar 
problem in view of the fact that the phases are practically exactly reproduced 
after 7 days. Even after 52 weeks there is very little shift. Hence, if one is not 
prepared to spread the observations of the star over a long period, say 5-10 years, 
it is very difficult to fill in all the gaps in the light curve. Even with our 1596 
observations of the star, all told, we found that when all the observations were 
plotted, small stretches of the curve remained which were not very well observed, 
whereas other portions were well determined. This difficulty is somewhat 
enhanced by the fact that cloudy weather conditions at the Boyden Observatory 
sometimes show a tendency to be repeated with a period of a week. 


TaBLe II 
Normal points of the blue light curve of V 499 Sco 

Phase Am, M.e. n Phase Am, M.e. n 

0°002 +oMm-273 +o™-o902 26 0569 o™-1902 + 0o™-904 22 
‘020 228 ‘003 21 “589 “185 "003 27 
"040 +O'lol "003 31 616 ‘208 ‘006 12 
063 —0'072 "005 12 ‘640 "224 ‘006 14 
‘089 "194 "004 10 667 253 "003 14 
-120 "225 "003 28 685 "264 "003 21 
138 241 "003 26 ‘700 ‘257 "004 23 
‘163 "254 "003 8 *716 ‘271 "003 20 
218 273 "002 8 -762 263 "003 7 
"244 *269 "002 20 -778 "265 "003 
‘266 285 003 44 ‘818 ‘260 "002 28 
286 278 ‘002 21 “841 "254 "002 36 
“319 263 003 8 $62 287 003 29 
384 ‘231 "004 19 "882 ‘218 "004 14 
“415 ‘190 +004 19 ‘QOI "204 "005 20 
438 —0'086 "003 18 "924 138 003 32 
“463 +0'076 "004 17 "942 0°046 "002 29 
"503 "185 "004 18 ‘960 +0060 ‘002 26 
"529 O'113 "004 25 ‘976 “188 "003 44 

0°550 —Oo'oll 0002 26 0988 +0°244 0°003 32 

Tas_e III 
Normal points of the yellow light curve of V 499 Sco 

Phase Am, M.e. n Phase Am, M.e. n 

0002 +om-644 + 0™-904 14 o'581 +ol-232 + 0-903 18 
"020 “604 "003 12 “602 *196 ‘006 12 
"039 ‘480 "004 16 ‘641 ‘175 ‘004 17 
"067 *307 "004 II 677 "149 003 2 
“101 “198 "004 19 697 153 "003 28 
‘121 "183 "004 16 ‘959 "129 "004 18 
-136 "159 "005 18 “748 138 003 10 
“155 "153 "003 15 ‘799 "144 "003 16 
*227 "144 "003 15 828 *147 "003 29 
"253 "138 "003 27 "847 *164 "003 26 
‘270 "132 "003 28 866 "172 "003 2 
-286 "139 ‘002 23 "894 *184 003 22 
-367 "165 “002 29 "922 -246 "003 26 
"425 ‘249 "004 32 ‘941 "348 ‘003 21 
‘457 "423 "004 19 “960 "459 "003 19 
"$03 "599 "004 20 ‘975 "564 "003 28 
536 "470 "004 40 0988 +0°629 0003 21 


0°563 +0O°319 0°003 22 
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In Tables II and III the normal points of the blue and the yellow light curves 
of V 499 Sco are given, computed in the usual way. In the first and second 
columns are given the phases and magnitude differences of the normal points. 
The third column gives the mean error for each normal point and the fourth 
column the number of observations used in forming it. The mean errors are 
no larger than we expected from previous investigations of the same kind (2). 
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Fic. 1.—Blue light curve of V 499 Scorpit. 
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Fic. 2.— Yellow light curve of V 499 Scorpii. 


Fig. 1 gives the blue light curve and Fig. 2 the yellow light curve of V 499 Sco. 
These curves clearly show that the primary minimum is somewhat deeper than 
the secondary minimum. ‘There seems to be no significant difference between 
the two maxima, a result which is confirmed by our preliminary light curves drawn 
through the means of all plotted observations. These preliminary curves differ 
from our published ones in that the maxima are more flat-topped and that they 
show evidence of shoulders in both cases. 

According to the Gaposchkins (3) V 499 Sco is a B-Lyrae star on account 
of its period being more than 12 hours. Our investigation reveals it as a B-Lyrae 
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star with two fairly flat topped equal maxima and evidence of shoulders in the 


light curve. Its components are presumably less elliptical than those of V 525 
Sag (4) for instance. 


In conclusion it is a pleasure to record our thanks to the staff of the Boyden 
Observatory for assisting with the observations of this star on a number of nights. 
Department of Mathematics, Armagh Observatory, 
University of Stellenbosch, Northern Ireland. 
Union of South Africa: 
1957 February 21. 
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Summary 
Three hundred and sixty-six southern stars of B, A and F types brighter 
than 6°5 mag. were classified on the MK system. Microphotometer tracings 
were used to estimate the line ratios. A mean modulus for the members of 
the Scorpio-Centaurus cluster is determined. The calibration of absolute 
magnitudes of B stars is briefly discussed. 


The spectrograms were obtained with a three-prism spectrograph at the 
Cassegrain focus of the Reynolds 30-inch reflector (dispersion 35 A/mm at Hy; 
range H{ to HB: AA3g00-4g00A). As the spectra were taken primarily for 
radial velocity measurements in another programme, they were too narrow for 
classification by visual inspection. ‘They were therefore recorded with a Hilger 
microphotometer and relative intensities of the pairs of lines used as criteria 
were estimated on tracings of small enlargement (x7). The extensions of the 
wings of the hydrogen lines were examined by a binocular microscope. The 
procedure was found most valuable for B and A stars and even for F stars which 
show a greater number of sensitive line ratios. ‘The main difficulties encountered 
were due to the plate granulation which affected the depth of faint lines and 
the differences in density. Dense spectra were recorded with two sensitivity 
ratios in order to obtain more details in the darker region. Another difficulty, 
common to all procedures, was the shortage of standards. 

The distribution of the 63 MK standards (1) recorded is given in ‘Table I. 
The main gap occurs among A type stars and for luminosities II and 1V among 
all types. 

The 366 southern stars brighter than 6-5 mag., 263 of which are south of 
— 30°, classified on the MK system (1), are listed in Table I]. ‘The columns 
give: 

(1) the name of the star ; 

(2) the H.D. number; 

(3) and (4) the galactic longitude and latitude to the nearest degree ; 

(5) the spectral type and luminosity class ; 

(6) the Harvard photovisual magnitude; an asterisk indicates that the 
magnitude has been corrected by +0-3 for a single-spectrum binary, or +0°6 
for a two-spectra binary ; 

(7) the colour excess derived from various colour indices all converted to 
the B-V system: the value of the intrinsic colours were taken from (1). ‘The 
main sources were the Cape Mimeograms (2) and the list published by J. Schilt 
and C. Jackson (3). The latter were corrected for lack of linearity by comparison 
with the converted B—-V Cape colours. A few other colours were found in the 
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TasBLe I 
MK standards recorded 

pe I II III IV V (n) 
O9 I 

O9°5 I I 

Bo I 2 

B o's I I 

Bi I I I 2 

B2 I I 2 2 

B3 I I I 

B 5 2 I I 

B6 I I 

B7 2 I 

B8 I I 2 

Bo I I (37) 
Ao I I 

Al 2 

A2 I I 

A3 I I I 

As I I 

A7 I (12) 
Fo 1 I I I 

F 2 I I I I 

F 6 I I I 

F 7 I I 

F 8 I (14) 


(n) (11) (7) (14) (7) (24) (63) 


work of O. J. Eggen (4) and some B-V indices were communicated by A. R. Hogg 
in advance of publication. When two or more different values were obtained 
for a star, the mean B—V was used; the agreement was good except in two or 
three cases (see notes at the end of the table). Negative colour excesses were 
disregarded ; 

(8) the modulus myp—-M=m-—3E,—M, where M is the visual absolute 
magnitude as given by Morgan and Keenan (5) and corresponding to the spectral 
type and luminosity. When no colour-excess is available the value of the modulus 
is given in parentheses; and 

(9) remarks and references to notes at the end of the Table. Recent classi- 
fications on the MK system taken from a variety of sources are given in parentheses. 

B stars.—The catalogue includes 218 B stars, 75 of which are members of the 
Scorpio—Centaurus cluster. Fifteen stars which seemed intermediate between 
B3 and Bs were classified B4, even though this subdivision does not exist in the 
MKK and MK systems. Type-differences in the B stars are shown mainly by 
the progressive increase of the hydrogen lines and the decrease of He1 lines: 
M4471, 4387, 4144 and 4026. The luminosity classification depends on the 
ratios 4009/4089, 4119/4144, 3995/4009, 4416/4387, 4481/4471 and the sharpness 
of the hydrogen lines. The classifications were made twice independently and 
the average difference between the two determinations is: 

0-4 of a subclass both for types and luminosities. A comparison with 
M. L. Woods’ catalogue (6), which gives the most extensive list of stars in common, 
indicates for about one hundred B stars the following average differences: 

o°6 of a subclass for types, o-5 for luminosities. 
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Name 
49 G Cet 
\’ Phe 
B} Tuc 
B? Tuc 
70G Phe 
x Sel 
fea Scl 
v Phe 
¢ Phe 
102G Scl 
12G Hyi 
x Eri 
129G Scl 
x Hyi 
B For 
@ Eri 
A Hor 
« Eri 
e Hvi 
a Cet 
93 Cet 
7 Eri 
1G Eri 
55G Hvi 
¢ Eni 
ge ‘Tau 
59G Hyi 
44 G Hor 
7G Ret 
T For 
17G Ret 
210G Eri 
4G Dor 
1G Cae 
x Dor 
253G Eri 
268 G Eri 
19G Cae 
e~=—s- Pic 
¢ Dor 
12G Col 
1i8G Col 
82G Lep 
74G Col 
94 G Lep 
90 G Col 
36 G Mon 
o Car 
74G CMa 


and luminosities of B, A and F southern stars 


HD 
2696 
2834 
2884 
2885 
4293 
5737 
617 
6767 
6882 
7312 
9896 
10144 
10538 
12311 
13709 


14225 
15008 
15233 
15371 
16978 
17081 
18853 
18978 
20121 
20313 


20320 
21364 
21722 
21981 


22252 


739 
346 
563 
604 
$246 


2 
5 
/ 
- 
/ 


ew NN NN 


29305 
29435 
30422 
31093 
31203 


33262 
34868 
36060 
40972 
41843 


42301 
44506 
45320 
45348 
49131 


49 
277 
273 


-59 


a 
7O 
60 
46 


36 


—17 
-20 
—O4 
—24 
—13 


Tase II 
Sp 
A3 V 
Ao V 
B8 V 
A2 IV 
Az II 


B8 III 
A2V 
A3V 
B6 V 
A7 III 


F2V 
Bs V: 
A2V 
Ag III 
Ar V 


B8 V 
A2V 
F2 IIIp 
Bs III 
Bg III 


Bo Ill 
B7 Ill 
A3 IV 
F6 Ill 
Fo II 


Am 

B8 V 
Fs IV 
A2V 
B7 V 


Ao V 

F2 Ill 
Bs III 
Fs IV 
F6 V 

Ao V: 
Bg IV-V 
A3 III-IV 
Ai: V 
Fo IV 
F7 V 

Ao IV 
Am 

Al V 
Ai: V 


Ao IV 
B2 Vn 
A7 IlIp 
Fo Iab 
B3 V 


45! 
m(HP,) E, my-M Remarks 
5°21 (3°41) 
4°82 +002 4°46 

4°8 
4°42 70°04 4° 
5°95 (5°65) 
4°34 —0O'05 7°34 
5°48 (4°28) 
5'10 0°13 2°91 
4°62* 0°07 5°‘51 
5°86 (5-56) 
6°05 +. (2°85) 
0°6 -O°OI 1°9 
5°60 “0°02 4°40 
2°95 +004 2°43 HD: Fo 
5°18 —0'02 4°38 
3°67 +O'Ol 4°14 
4°20 —O'OI 3°00 
5°40 o'o1 4°60 
4°61* +0°03 7°72 
4°16 +002 6°10 
4°54* 0°03 7°55 (B7 V) (1) 
5°65 +005 8:50 
4°03 O‘II 2°70 
5°89* 0°07 4°89 
5°63 (7°63) 
4°80 ze H lines A3 
4°22* 0°03 4°63 
5°88 (3°18) 
5°93 O12 4°37 
5°68 0°07 6°27 
6°04 O’OlI 5°71 
6°07 (5°27) 
5°77 (8-97) 
5°95 0°02 3°22 
6°12 (2°32) 
3°39 0°05 3°09: 
6°24 o'o 66°74 
6-00 (5°70) 
5°75 +o'o9 4°68  K line A2 
5°46 (3°46) 
4°68 +O°05 0°53 
5°73 —o'o2 6°13 or Bo? 
5°79 cats re K line A3 
5°90 -0°02 5°10 
5°71 oo 6491 
5°41 —o'o2 5°81 
5°86* +002 840 (Br°s5 V:) 
5°69 (5°39) 
—o°79 —=O°83 375 
5°83 —o°02 7°83 
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Name 
80 G Mon 
¢ Men 
101 G Mon 
120 G CMa 
122 G CMa 
6 Mon 
158 G CMa 
d Vol 
96 G Pup 
108 G Pup 
f Pup 
3 1 Pup 
161 G Pup 
o Pup 
P Pup 
b Pup 
D Car 
16 Pup 
e Vol 
280 G Pup 
q Pup 
x Cha 
48Ge Vel 
o Veil 
d Car 
& Vel 
f Car 
ec ‘Cor 
b Car 
1Gw Vel 
x Vol 
G Car 
105 G Vel 
a Car 
t Car 
z Vel 
k Vel 
B Car 
¢ Car 
Kk Vel 
b Vel 
32 7? Hya 
145G Vel 
m Car 
O Vel 
157G Vel 
v CarA 
165G Vel 
166G Vel 
d Vel 


HD 
49043 
50506 
52312 
54031 
54309 


au m 
© 
~— 


85 


Ouuw wv 
ool 
mun 

uw N 
mb OW - 


61330 
62623 
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66591 
67797 
68520 
69253 


70060 
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73634 
74195 
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74956 
75311 
76728 
77002 
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73045 
78791 
79275 
79351 
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79735 
79940 
80007 
80404 
81188 


82434 
$2446 
82984 
83944 
84461 
84816 
85123 
85953 
85980 
86440 
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HD: Ao 


(Bo V: pe) (2) 
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Name 

v* Hya 

a Sex 

183 G Car 
q Vel 

w Car 

J Vel 

i Cee 

oe (er 

p Car 
a1at Car 
6 Car 
228 G Car 
229 G Car 
236G Car 
238 G Car 
52 Cha 

i Vel 

5p Ce 

y Car 
74 Leo 
a Cen 
sy Crt 
3G Cen 
A Cen 
50G Cen 
A Mus 

j Cen 

18 G Mus 
1G Cru 
0 Cru 

62 Cru 

» Cru 
92G Cen 
5 Cen 

p Cen 

56 Cru 

Y Crv 

B Cha 

¢ Cru 
113 G Cen 
25G Cru 
119 G Cen 
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y Mus 
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x Mus 
45 G Cru 
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p Cen 

n Cen 


and luminosities of B, A and F southern stars 
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288 
282 
289 
294 


286 
301 
293 
282 
284 


294 
295 
301 
293 
296 


295 
296 
289 
283 
299 


“Tm 


Tase. [I—(cont.) 


b 
Tr O03 
+05 
rO5 
ri4 
+33 
+56 
+10 
T25 
rol 
+08 


+28 
+ 19 
+318 
ri4 


2 


+27 
+08 
+29 
+18 
+16 


+15 
+o1 
+19 
+o1 
+03 


—18 
+13 
+04 
+19 
+08 


+15 
+-og 
—05 
TII 
+20 


+03 

+28 
+133 
—o6 
—03 


r13 
+13 
+22 


+13 


+10 
+10 
—ol 
—10 
+13 


Sp 


Bs: Vn 


Bs Ve 
B3 IV 
B4 1V 
B2 1V 


At III-IV 


Bs Ill 
A2V 
Bs V 
Bi V 


F2 Ill 
B2 V 

B3 Ve 
Bg Ill 
Bs IV 


B3 IV 
B3 Ill 
A3 IV 
AlV 
Bg V 


m(HP,) 


4°70 
5°32 
4°II 
4°81 
4°56* 
4°66* 
5°35 
2°90 
4°50 
2°41 
4°54* 
3°76* 
3°24 
6-02* 
4°59 
6°04 
5°94 
5‘00* 
3°96 
4°09 
4°61* 
1'02* 
4°47 
6°43 
5°20 


5°07 
4°02 
4°28 
4°51 
4°45 
a9" 
4°01 
3°34 
3°06* 
4°06 


6-09 
5°62 
4°43 
6-o1 
5°27 


3°03* 
3°58" 
5°37 
(4°7) 
5°92 
4°89* 
5°16* 
4°12 
2°98 
5°97* 





Antoinette de Vaucouleurs, Spectral types 


Ey 

—o°or 
~~ 0°02 
Tool 
70°02 
+ 0°02 


—o'06 5:16 


oo 
—0°02 
+ 0°'02 
70°03 


—o'ol 
oo 
o’o 
o"o 

+o°'o! 


+ o°o!l 
TOI! 
ro°ol 

o'o 
—0°02 


+O"1o 
+ 0°03 
+0°02 
+oOrl4 
+ 0°07 


+007 
role) 

+0°21 

—o'06 
0°03 


—0°02 
=>" Ot 
—O°OS 
ro'os 
+O°OlI 


+ 0°09 
+ 0°04 
—o'o2 

°'o 
+0°'09 


°'o 
+ 0°02 
+0°O3 
o'o 
+0°02 


+o°o1 
o'o 
7oO°o!r 
—o'o5 
—0°06 5°97 





yp 





ee a ined Nn at ee 


No. 4, 1957 
Name 
6 Lup 
e Lup 
¢* Lup 
k Lup 
100 G Lup 
35 ¢* Lib 
Y Lup 
90 G Lib 
d Lup 
t Lib 
g Lup 
ys? Lup 
128 G Lup 
B SerA 
14G Nor 
x Lup 
1(b) Sco 
B TrA 
2(A) Sco 
4G Sco 
5G Sco 
p Sco 
67 Sco 
7 Lup 
151 G Lup 
B ScoA 
6 Lup 
229 G Sco 
13 Sco 
v Sco 
o Sco 
52G Nor 
e Nor 
7X Oph 
26G Oct 
71 G Sco 
N Sco 
62G Nor 
pn’ «Sco 
p? Sco 
101 G Sco 
117G Sco 
32 G Ara 
k Sco 
116G Sco 
80G Oph 
7 Sco 
130 G Sco 
46G Ara 
40 € Oph 


and luminosities of B, A and F southern stars 


HD 
136298 
136504 
136664 
137058 
137432 


138485 
138690 
138764 
138769 
139365 


139664 
140008 
140784 
141003 
141318 


141556 
141637 
141891 
142114 
142165 


142184 
142669 
143018 
143118 
143699 


144217 
144294 
144892 
145482 
145502 


147084 
147152 
147971 
148184 
148542 
148688 
148703 
149711 
151890 
151985 


152234 
153613 
153716 
154090 
154368 


154481 
155203 
155450 
156838 
156897 


323 


320 
302 
304 
326 
274 
308 
314 
308 
314 
315 


311 
320 
300 
318 
318 


325 
312 
321 
296 
331 


TaBLe II—(cont.) 


b 


+13 
+10 
+16 
L144 
+16 


+30 
+] 
35 
+o8 
- 20 


+o8 
H15 
-415 
- 47 
-O2 


Las 
+21 
—o8 
20 


r2iI 


+21 
17 


+03 
+o8 
+Or 
+03 
+03 


foie) 
+ oS 
—fI 
+03 
+02 


+03 
—16 
+07 


Sp 
B3 IV 
B3 IV 
B3 IV 
Ao IV 
B4 V 


Ao III-IV 
B3 V 

Fo V 

B3z Vn 

B6 Vn 


B3 Vne? 
B3 IV 
B2 IV 
B3 V 
Bs V 


Bi V 
B3 IV 
F7 V 
B3 Vn 
B3 V 


As Ill: 
B6 V 
B3: V 
B3 V:e 
A3 IV 


B1 Ia 
B2 V 


B4 V 
Br II-IIle 
Og°5 I-II 


Ao ITI-IV 
F2 IIl 

B1 Ill 

B2 V 

F2 V 


m(HP,) 


3°33 
3°95* 
4°66 
4°64 
5°48 
5"90* 
3°48* 
5°16 
4°74 
4°00* 
4°60 
5°38* 
5°57 
3°57 
5°72 
4°69* 
4°74 
2°97 
4°59 
5°36 


5°39 
4°61* 
5°99 
3°74" 


E, 
—0'03 
0'0 
+-0°03 
+0°02 
+902 


oo 
—0°02 

o"o 

o"o 
—0°02 


0°03 
—0°02 
+ 0°04 

roHre) 
O17 


0°04 
7 oO°og 
—0'O2 
+ 0°08 
+O"1o 


+Orrl 
—0o-'02 
+O°02 
—0'03 
role) 


0°63 


Porrl 
+0°40 


0°47 
+ O°OS 
+O'l2 
—0°02 
+o°ol 


+033 
—o'o2 
+o'o6 
rO°43 
Lo-71 


—o°06 
+0°03 
+0°20 
+0°04 
oho) 
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m,-—M Remarks 


6°13 

‘05 
‘37 
18 
"02 


a | 


~sI > 


9°60 
*26 
55 
‘60 


masts wu 


2°01 
6°48 
6°55 
a ae 
9°31 
5°49 
6°47 
oe) 
6°35 


Oo f=) “NI oO 
Nuss O Lo 
AAAs N Oa 


nnwn 
rn Of 
> WwW 


wn 
-— 


6-62 
5°90 


2°84 
6°46 
6°85 
5°34 
(5°05) 
10°98 
7°06 
8°93 
5°74 
6°94 


9°40 
5°53 
7°16 
8°37 
10°04 
6°86 
2°47 
9°62 
8-21 
I°Is 


(Az IV) 


HD: Bo 


(B3 V) 


(B2 V) 
(Br V) 


(Bo-s V) 


(Bz IV ?) 
(As II) (6) 


(B2 IVp) 


(Bx Ia) 


(Bis V) 


(Bo-s Ia) 


(Og°5 lab) (7) 


(B2 IV) 


31* 
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Name 


8 Oph 


167 G 
77G 


13G 


66 


7 


102 


6G 
45G 


9G 
77G 
83 G 


91 G 
94 G 


32 G 
31 G 
A 


134 G 
A 


45G 
47G 
162 G 


186 G 
63 G 

55 
70G 


228 G 


73G 
51 
75G 


Sco 
Ara 
Ara 
Ara 


Sco 
Sco 
Sco 
Oph 
Ara 


Sgr 
Ser 
Oph 
Pav 
Her 
Tel 
Ser 
Ser 
CrA 
Tel 


Sct 
Tel 
Tel 
Sct 
Ser 


Ser 
CrA 
Ser 
Sgr 
Pav 


Sct 
CrA 
Tel 
CrA 
Ser 


Ser 
Aql 


CrA 
Ser 
Ser 
Sgr 
Tel 
Sgr 
Pav 


Sgr 
Tel 
Pav 
Aql 
Pav 


HD 


157056 
157243 
157246 
158094 
158427 


158926 
159358 
160578 
161056 
161783 


163685 
164259 
164284 
165040 
166182 


167128 
167356 
167833 
168905 
169467 


170296 
170465 
170523 
170740 
171034 


171856 
172777 
172910 
173117 
173948 


175156 
175362 
175510 
175813 
176687 


177074 
177756 
178322 
178628 
180885 


181869 
182681 
183806 
186005 
186219 


186500 
186543 
186837 
187532 
187653 
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TaBLe II—(cont.) 


l b Sp m(HP,) E, m-M 
328 +05 B2IVv 3°61* oor 6°88 
313 —0o6 B6V 5°03 0°07 5°92 
302 —12 Br Ill 3°36 +O°12 7°30 
299 —15 B8&V 3°68 —o'03 4°18 
308 —10 B3 Vn 2°82 oo §644 82 
319 04 B2iv 1‘96* oO) «6526 
342 10 B8V 5°50 (6-00) 
319 06 «6«=B2IV 2°62* +0'02 5°86 
346 +10 B3 Vn 6-50* +0°42 7°24 
306 —14 B3V 6:05* +0°05 7°90 
329 —0o4 B3I1V 5°94 -o'10 8-94 
351 08 F3V 4°60 —0'03 1°30 
358 12 B2Vn 4°75 +O°Il 7°02 
299 —21 Ap 4°69* O'1s ia 
15 17 B2V 4°24 (6°84) 
306 19 B3V 5°71* +o'11 6°38 
340 —o2 Az2l 6-10 9°02 10°74 
348 +o1 AsglIill 6°24 (6°24) 
318 —15 B3 Vn 5°34 +oO°'OI 7°31 
316 —17 B3zIIl 3°68 +0°02 7°32 
345 04 A2V 4°60 —O'OI 3°40 
316 17 B6I1V 5°28* +002 6:22 
316 17 Bs IV 5°50* °o'o 67°70 
349 o2 B2V 5°66 0°40 7°06 
329 13 Bz IV 5°34 0°04 8-02 
340 —o8 A7IIIp 5°75 es ree 
325 1 ArV sas O'll 5°62 
327 15 B3V 4°82 oo) §6—66-82 
338 —11 Bs: V 5°76 +o18 6-52 
301 25 Br Ve 4°28 O°3r 7°IS 
347 —0o9 B3 Ill 4°93 +O-3I 7°70 
327 18 B8& IV 5°42 —0°03 7°12 
311 23 Bolll §°23* +0°03 7°14 
327 19 FoV 4°84 0°07 1°63 
334 17 A2V: 2°98* -+o-o1 1°75 
334 —17 AolV 5°53 +006 5°75 
358 -—o7 B7V 3°34 +0705 3°89 
323 23 Bs V 5°80 +008 6-86 
326 —22 B6V 6°48* -+o-10 7:28 
330 —22 B4sIV 5°62 +O°O5 7°97 
325 24 Bo Ill 4°34" —0°04 6°34 
336 -—21 BoV 5°58 +008 5°34 
321 —27 Ap 5°79 oo 
352 —20 Folll 5°07 (4°47) 
289 —31 Am 5°47 ee 
336 -—26 B8 Vn 5°53 +o12 5°67 
308 —31 Asp 5°48 rae ba 
303 —31 Bs V 6°27 OO 67°57 
358 —20 FolV 5°49 (3°49) 
303 —32 A3V 6°27 (4°47) 
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Remarks 
(B2 IV) 
(8) 
(B2-5 V) 
(MKK Std) 


(F3 V) 
(B2 Ve) 
Am? 


(B2 V) 


(Ao Ia) 


(A3 V) 


(B2 V) 


Am? 








el tl ein 


Pi aia hii sees 


No. 4, 1957 
Name 

76G Pav 
79G Pav 
61 Sgr 
259G Ser 
x’ «Ser 

%’ Pav 

v Pav 

x Del 
52G Cap 
64 G Cap 
x Oct 
8G Ind 
x Cap 
23G Ind 
e Mic 

y Pav 

18 Agr 

Y Ind 

6 PsA 

€ Aqr 

5 Peg 
54GB Oct 
A Cap 
20G Gru 
Oo Ind 
66G Oct 
U] PsA 

o Aqr 
25G PsA 
« <Aqr 

c Peg 

A PsA 

p Aqr 
48 y Aqr 
m7 Gru 

6 Tuc 
40G PsA 
72G Ind 
570 Aqr 
7 <Agqr 
49G PsA 
€ PsA 

€ Gru 

52 Peg 

a PsA 
80 G_ Ind 
a Peg 

83 Aqr 

v Gru 


and luminosities of B, A and F southern stars 


HD 
188097 
188164 
188899 
1891938 
193571 


195207 
196519 
196775 
196867 
198174 


199443 
199532 
199623 
201184 
202103 


2026027 
203605 
203705 
203760 
204554 


205767 
205852 
206553 
207052 
208321 


208450 
208741 
209014 
209409 
209522 


209819 
210027 
210934 
211838 
212061 


212132 
212581 
213135 
212728 
213320 


213998 
214484 
214748 
215789 
217232 


217792 
217831 
218045 
218060 
218242 


351 


2gI 
352 
288 

21 


35 
340 
353 
304 

53 


336 
283 
57 
36 
325 


TaBLe II 
b Sp 
2 Am 
2 A3p 
23 A2lIV 
2 A7IIIp 
36 A1V 
37. FoV 
36 6 B8 V 
17 B4 Vn 
17 B8&V 
38 B6V 
—36 As IlIp 
34 ~Fs5 IlIp 
2 F6V 
—41 AoV 
-44 Ag Ill 
—45 A2p 
—41 F5 V 
-41 AgIlIl 
—45 Folll 
—48 A2lIV 
-41 A7IV 
—25 FolV 
—29 A7IV 
—45 A2V 
—53 A3V 
—49 FolV 
—37 F3 lll 
—38 B8&V 
—44 BO: Vn 
—55 Bs Vn 
so 6 B8IV 
-25 Fs V 
57 B8 Ill 
—51 B8&V 
47 AolV 
57 FolVp 
—47 B8&V 
—60 FoV 
—45 A3V 
—54 ArlV 
—49 B8&V 
—62 A2Vp 
—63 B8&V 
—57 A2V 
—43 FoV 
—67 FolIV 
—46 F2IIIp 
—41 BoV 
—59 F2V 
—67 ArV 


(cont.) 


m(HP,.) E. 
577 
6°42 wet 
5'O1 0°05 
5 "90 0°13 
5°57 “0°02 
4°78 —o'ol 
5°26 + 0°02 
5°98 
3°61 0°03 
5°74 rorro 
5°83 
5°17 
5°80 
5°25 70°05 
5°80 
4°73 0°03 
4°22 0°02 
5°54 or 
6°17 + 0°04 
5°95 
4°71 —o'ol 
5 “‘RR* 
6°50 
5°46 
5°50 
5‘29* —0O'O2 
5°34 on 
5°37 o°ol 
4°51 0°09 
5°84 —0'03 
4°84* 0°05 
4°35" oo 
5°31 0°02 
5°31 0°04 
4°07" 0°05 
5°76 0°09 
4°7° 0°07 
5°89 
5°65 ner 
5°09* 0°09 
3°95 o'o! 
6°15* 0'0 
4°26 0°02 
3°63 o'ol 
5°68 
5°07 
5°57 oe 
2°86* 0°03 
5°83* 0°07 
5°54 —0'02 


m,-M_ Remarks 


(7°58) 
4°02 
6°54 

(5°88) 

(4°17) 

(2-00) 
4°80 

80) 


wv 


‘66 
14) 
45 


i) 


_~ 


_~ 


muumun ot 


3°01 
(3°88) 
(4°80) 
(4°26) 
(3°70) 


3°2 
(4°94) 
84 
34 
14 


aU wn 


“I 


39 
‘85 
‘31 


5°69 


“oo 


3°49 

4°99 
(2°89) 
(3°85) 


5"19 


4°42 
4°95 
4°76 
2°40 


(2°68) 


(3°07) 

(4°77) 
2°77 
2°63 
4°74 


“44: 


(Bg V) 


(A7 V) 


(B8 V) 


(Bo V) 
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TaBLe I]—(cont.) 


Name HD l b Sp mHP,) E, m-M_ Remarks 
88 G Gru 218630 316 —66 F6V 5°68 .-. (1°88) 
59 Peg 218918 54 —47 A2V 5‘10 «ee (90) 
%* Agr 219832 37 —63 AoV 5°13 es (4°83) 
107 G Gru 220729 295 —61 F3IV 5°48 ... (2°88) 
Scl 221507 322 —72 Bolll 4°49 —o'05 6°49 
11G Phe = 222095 302 —68 A2V 4°82 +O°OI 5°99 
w? Aqr 222661 39 —71 BoV 4°45 t-o-o5 4°30 @©6©(Bg'5 V) 
106 Aqr 222847 30 —74 B8&V 5°19 +0704 5°67 
5 Scl 223352 352 —78 AoV 4°61 +003 4°22 
27G Phe 224022 305 —73 F8V 591 = 1°81 
35 G Scl 224113 337. —78 Bs IV 6-00 -+0:06 8-02 
42G Tuc 224361 279 —54 A2p 5°98 Sie as 
e Tuc 224686 277 —51 B7V 4°61 +0°05 5°26 
C Sel 224990 345 —80 B4Ill 5‘21* +0°04 8°59 
a Cet 4295132 44 —76 Ao: IV 4°64 —o'o2 5:04 (Bg IV) 
45G Tuc 225253 274 —46 B8 III 5°54 —0'03 8°54 
NOTES 


(1) wm Cet—the classification B7 V given by A. Slettebak and R. F. Howard (Ap. 7., 
I2I, 102, 1955) seems difficult to justify ; the star has sharp lines and is indicated as Bs 
I or II by S. L. Lippincott and D. Hoffleit (A. 7., 58, 167, 1953). 

(2) o Pup—four values of converted B-V colours were found as follows: -—o-o09 
—0'20 —135 —0°025. 

(3) ¢ Vel—Bs I or II by S. L. Lippincott and D. Hoffleit. 

(4) y Car—classification Fo Ia by W. P. Bidelman (P.A.S.P., 66, 249, 1954) ; A 4226 
weak for F2, strong Fe lines 4233 and 4172-78. 

(5) HD 125288=HR 5358, Bs I or II by S. L. Lippincott and D. Hoffleit. 

(6) o Sco—class III doubtful, faint spectrum. On the basis of its proper motion and 
position the star was recently found to be a member of the Scorpio—Centaurus cluster 
(Ap. F., 12%, 557, 1955) with M, = —3°6 and a provisional classification of As II. 

(7) HD 154368—a B-V colour of +0-605 which would give a total absorption of nearly 
3 mag. was rejected. Schilt and Jackson’s converted B-V is +-0°37 giving a total absorption 
of 2:1 mag. in agreement with the value 7 E,=2-1 mag. obtained by P. Th. Oosterhoff 
(B.A.N., 11, 299, 1951). 

(8) y Ara—three values of converted B-V colours were found as follows : 0°17 
—O'2I —0°I0;. 

A referee has called my attention to a recent Cape Mimeogram* (not available at 
Mt Stromlo) which contains a few MK types for some stars in Table II : 

HD 4293 A7 V 
33262 F8 V 
128898 Fo III (Astrontium star which shows a false absolute magnitude 
effect) 
176687 A2 Ill (The Stromlo spectrum was very dark and the lumin- 
osity class is uncertain) ' 
202103 As V 
203608 F8 V 
213135 Fo V 


*Fundamental data for southern stars—First List. 


A stars.—In A stars the helium lines have disappeared and the metallic lines 
and the K line have become progressively stronger. Useful line ratios are 
4481/4385, 4128-32/4030-4 and 4300/4385. Luminosity differences are shown 
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by the intensity of the blend 4171-4179, the Feit line 4233, the appearance of 
the hydrogen lines and the ratio 4417/4481. Among 94 A stars, 7 have been 
classified as metallic-line stars or peculiar. Because of the false absolute magnitude 
effect observed in these stars, a few, possibly two or three, may have been missed 
and classified A5 III or A7 III. In such cases the K line, lost in the granulation, 
may have been too weak to remove the ambiguity. The average difference 
between the two determinations is about one subclass for types and luminosities. 

F stars.—For F stars type-differences are shown mainly by the rapid decrease 
of the hydrogen lines and the increasing strength of the lines 4325, 4226, 4045 
and the G band. The line ratios used for luminosity classifications are 
4417/4481, 4444/4481, 4171/4226, 4063/4077, 4045/4077. The very luminous 
F stars are distinguished by the appearance of the G band, the strength of the 
pair of lines 4172-4179 and the numerous sharp metallic lines. Table II contains 
54 F stars for which the average difference is 0-5 of a subclass both for types and 
luminosities. ; 

Cepheid 8 Dor.—In addition, eight spectra of the cepheid 8 Dor (7, 8), taken 
by S. Gaposchkin, were recorded and classified and the variation with the phase 
was found to be as follows: 


Plate JD 2,430,000 + Phase (in days) Sp 
1477 5740°205 7°79 Fo la 
1496 5769101 ay i: F8_ Ia 
1498 5769°128 1°18 F8 Ia 
1511 §771°222 3°47 Go +1 
1523 5773°140 5°19 G2 | 
1525 5773°177 5°23 G2 I 
1550 5777149 9°17 F6-8 | 
1552 5777°144 9°19 Fo Ia 


Scorpio—Centaurus cluster.—The certain and probable members of the 
Scorpio—Centaurus cluster according to A. Blaauw (9) are listed in ‘Table III. 
Column 1 gives the name of the star, column 2 the HD number, column 3 the 
spectral classification and column 4 the modulus m,— ™ as in Table II. 

The mean moduli are: 


(a) certain members 6:87+0°09 p.e. (n=34), 
(5) probable members 65gto13 pe. (n=38), 
all members 6-72+O011 p.e. (n=72). 


Three emission line stars giving small distance moduli, 5Cen, » Cen and 
a Ara were rejected from the means. Emission line stars are not supposed to 
be intrinsically brighter than normal stars but because of broad diffuse lines 
the luminosity classes may be uncertain. 

The 72 members were grouped according to spectral type as follows: 


Bo, 1, 2 B3, 4 Bs-8 

mo-M 6°68 6°72 6°79 
p.e. +O0°15 +org +018 
(n) (20) (33) (19) 


No systematic effect in relation to the spectral type is evident. 
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Although the cluster contains few giants the members were divided into 
three luminosity classes as follows: 


Name 


Certain members 


a Car 
166 G Vel 
1G Cru 
92 G Cen 
5 Cen 
Cru 
Mus 
Cru 
Cen 
Cen 
p Cen 
3k Cen 
283 G Cen 
4h Cen 
d Cen 

x Cen 

a Cen 

a Lup 

b Cen 

o Lup 

x Cen 

45 G Lup 
5 Lup 

¢*? Lup 

90 G Lib 
%? Lup 

p Sco 

7 Lup 

a Sco 

N Sco 

t Sco 

pe’ Sco 

pe? Sco 

32 G Ara 
ot G Ser 
31 GCrA 


= >>R oO 


IV 

7°02 
+o°1l 

(29) 


Tas_e III 


Scorpio-Centaurus cluster 


IIl 
mo-M 7°58 
p.e. + 0°27 
(n) (6) 
HD Sp my—M 
79351 B2Iiv 7°18 
85980 B3 V 7°57 
103884 Bz V 7°48 
105382 B6III-IV 7-28 
105435 B3z Vne (4:6) 
106983 Bz IV 6°88 
109668 B3 IV 5°59 
112078 Bs: Vn 6-00 
113703 B4IV 7°25 
120307 B2V 6-36 
120324 B3 Ve (5:2) 
120709 Bs IV 6-76 
120908 Bs V 6-91 
120955 Bs IV 7°17 
121743 B2V 6°56 
122980 B2V 7°O1 
125823 Bé III 7°61 
129056 Br III 7°21 
129116 B3zV 6:03 
130807 _B6 III 7°53 
132200 Ba III 7°62 
132955 B4gIV 7°78 
136298 Bz 1V 6°13 
136664 B3z 1V 7°37 
138764 B6IV 7:26 
140008 B6 V 6°48 
142669 Bz 1V 6°77 
143118 B3zV “56 
147165 (Br III) 6-8* 
148703 B2V 7°06 
149438 (Bo V) 6-65t 
151890 B3 Vp 5°74 
151985 B2IV 6°94 
153716 B4V 7°16 
172910 B3zV 6°82 
175362 B8& IV 7°12 
Mean 6°87 
+ 0°09 pe. 


* MK Std E, = +0°26 
t MK Std E,= +005 


Name 


Probable members 


b Car 
i Car 
145 G Vel 
6? Cha 
a Cen 
33 G Cen 
18 G Mus 
p Cen 
8 Cru 
119 G Cen 
Y Mus 
45 G Cru 
B Cru 
pe? Cru 
pe’ Cru 
&* Cen 
202 G Cen 
J Cen 
e Cen 
v’ Cen 
B Lup 
A Lup 
100 G Lup 
Y Lup 
d Lup 
I 


a Sco 
151 G Lup 
B’ Sco 
6 Lup 
vy Sco 
x Oph 
22 Sco 
@ Oph 
a Ara 
9 G Sct 
45 GCrA 


{MK Std E, 


Vv 

6°35 
+O°1s 

(37) 

HD Sp m,y—M 
77002 B3 IV 7°73 
79447 B3 IV 6°83 
82984 B4 Vn 7°23 
93845 B3V 6°44 
98718 B6 Vn 5°84 
99556 Bs IV 7°46 
103079 B4 IV 7°63 
105937 B4V 560 
106490 B2Iv 6°23 
108257 Bs Vn 6-20 
109026 Bs V 5°21 
110956 Bz IV 7°39 
111123. Bo-s IV 5°53 
112091 Bs Ve 6°56 
112092. Bz IV 6°88 
113791 B2Iv 7°80 
115823 Bs III 8°75 
116087. Bs V 5°74 
118716 Br V 5°52 
121790 B3 IV 6°89 
132058 B2V 5°63 
133955 B3 IV 7°66 
137432 Bg \ 7°02 
138690 B3 V 5°48 
138769 ©=6B3 IV 7°55 
1393605 B4V 5 ‘60 
141637 B3zV 6°47 
142114 B3z Vn 6°35 
143018 B2IV 6°76 
143699 Bs V 6:26 
144217 Br V 5°43 
144294 B3 IV 7°04 
145502 B3z V 5*90 
148184 Bz V:e 5°84 
148605 (B2 V) 7°22} 
157056 B2IV 6°88 
158427 Bz Vne (4°82) 
170740 B2vV 7°06 
178322 Bs V 6°86 

Mean 6°59 
+ 0°13 pe. 
+008 
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A systematic effect depending on luminosity is well marked. It is partly due to 
errors in the classification and partly to errors in the calibration of standard 
absolute magnitudes. The number of objects of class III is small but even 
class 1V with a comparable number of stars with that of class V gives a mean 
modulus 0°67 magnitude larger. 

A. Blaauw (10) has recently determined from the proper motions of the 
cluster members (and two other associations), mean visual absolute magnitudes 
per volume of space, for types B1, 2, 3 of class V and Bz class IV. With these 
new values for 27 members of the corresponding spectral types, the following 
mean moduli are obtained : 


B2 IV Br, 2,3 V 
my-M 6°37 6:07 
p.e. +o0°21 + O15 


(n) (6) (21) 


The difference between the mean moduli of class IV and class V is reduced to 
0°3 magnitude. This results mainly from the fact that the Am between Blaauw’s 
absolute magnitudes for classes B2 IV and Bz V is only of 0-3 magnitude while 
it is of o’7 mag. for Morgan and Keenan standard absolute magnitudes (5). 
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Fic. 1.—Relation between absolute magnitudes and equivalent widths of H,, for B stars (MK standards). 


Absolute magnitudes of B stars from Hy line intensities.—Absolute magnitudes 
of B stars may be determined from equivalent widths of Hy and R. M. Petrie (11) 
has obtained a mean relation between luminosity and Hy intensity which shows 
a systematic difference from the absolute magnitudes of Morgan and Keenan 
(Victoria~Yerkes = + 0-7 mag.). 

Among the B stars measured for equivalent widths by Petrie there are only 
19 MK standards, but the extensive list published by S. Giinther (12) contains 
46 MK standards for B stars. In Fig. 1 the equivalent widths of Hy in angstroms 
(EW) for these stars are plotted against the mean visual absolute magnitudes of 













462 Spectral types and luminosities of B, A and F southern stars Vol. 117 


Morgan and Keenan. The supergiants and bright giants, classes I and II, are 
clustered around M=—6 and EW between 1°5 and 3A; the class III giants 
forma distinct sequence from M=— 5 and EW ~3°5 Ato M— =2 and EW10A; 
the mean points for classes IV and V fit closely the mean relation obtained by 
Petrie, except of course where this relation reaches the group of the very luminous 
stars. The sequence of the class III giants is from 1-2 to 1-4 magnitudes above 
Petrie’s relation. 

The equivalent widths published by J. Stock (13) for some MK standards 
for B stars show the same effect. 

This may explain, at least in part, the disagreement between the Victoria 
and Yerkes luminosities. Either Petrie’s relation is valid only for the Main 
Sequence and class IV, or the mean absolute magnitudes of the B stars of class ITI 
in the MK system are too bright. The systematic effect depending on luminosity 
classes in the Scorpio—Centaurus cluster tends to indicate that the latter 
alternative is more probable. But it is equally possible that the mean relation 
between absolute magnitudes and equivalent widths of Hy, being based on only 
one parameter, is approximately a Main Sequence relation for B stars. It does 
not show up the finer distinction of luminosities which can be made by visual 
inspection of the spectra or their tracings. 

Acknowledgments.—The author wishes to express her thanks to Dr Gollnow 
and Dr Buscombe for the use of the observational material in their programme 
and particularly Dr Buscombe who has organized the necessary observations of 
the MK standards, and for some helpful suggestions. 
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THE EFFECT OF DEPARTURES 
FROM THERMODYNAMIC EQUILIBRIUM ON THE 
INTENSITY RATIO OF 2D-‘S DOUBLETS 
IN GASEOUS NEBULAE 


Ali M. Naqvi* and S. P. Talwar 
(Received 1957 January 18)t 


Summary 


A revised calculation of transition probabilities in ground configuration 
of [O11], [S11], and [Atv] is made. Using these transition probabilities and 
allowing for collisional effects, the relative level populations in 2D and 2P terms 
I (D5)2—Ss,2) 
I (D3, 2- Ss 2) 
electron density N, and the electron temperature T,. The theory is made use 
of to estimate N,/7T,'/? for regions of maximum emission of [O11], [S11], and 
[Atv] in nebulae for which the intensity ratio r of *D-‘S doublets are 
observed. It is suggested that the accurate measurements of this intensity 
ratio in a nebula, together with the use of theoretical calculations, give the 
electron density in the nebulae. This may prove to be a very simple and 
powerful method for determining the electron density. 


and the intensity ratio, r , are calculated as a function of the 





t. Introduction.—'The doublets (?D5). 52-45) in [Ou], [Su], and [Arv] 
are observed in numerous emission nebulae. The intensity ratio 
I?D5 2 >4S) 

T@D3.—>*S) 

for [O11] has been calculated by Aller, Ufford and Van Vleck (7) and by 
Garstang (11), assuming the 7D levels to be populated in the statistical weight 
ratio 6:4. The values obtained for [O11] were 0-58 and 0-47 respectively. 

The amount of accurate data on the intensity ratio r for the ions under 
consideration is somewhat limited, yet it can be said that the variation in the 
values of r for an ion from nebula to nebula is real (see Table I). There are 
marked variations for [O11]; the value for NGC 40 is 0°85, whereas for IC 4997 
it is 0:38. For IC 4997, which is a “stellar” nebula, and for which the 
measurements are probably difficult, two independent measurements were made 
by Aller, Ufford and Van Vleck (hereafter abbreviated as AUV) which gave r 
equal to 0-36 and 0-40. 

The value of r for the corresponding [S11] doublet has been calculated by 
Garstang. The theoretical value r = 0°55 is significantly smaller than the observed 
values in gaseous nebulae. 

The assumption of the Boltzmann distribution for level population is not 
justified for most of the nebulae. ‘This assumption cannot account for the 
observed variation of the intensity ratio. As we shall see later, the Boltzmann 


v= 


* Now at University of Sind, Pakistan. 
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value must always be smaller or at most equal to any of the observed values. This 
however, is not the case with the theoretical values of AUV and Garstang. ‘The 
Boltzmann value depends only on the transition probabilities. Earlier, one of 
us (Naqvi (13)) had calculated the transition probabilities of [O 11] lines, including 
the effects of the mutual magnetic interactions. Using Nagqvi’s transition 
probabilities, we recalculated the Boltzmann value of r for [Om] doublet. The 
value thus obtained, r=0-34, is lower than the observed value for any nebula. 
When the Boltzmann distribution does not apply, the value of r depends on the 
electron density, N,, and the electron temperature, 7), of the nebula. In the 
present paper, the departures from the Boltzmann distribution are considered 
in detail. 


TaBLe I 
Observed intensity ratios of the *=D-—‘S lines 

Nebulae r (O11) r (S11) r (A tv) 
NGC 7662 0°62 (24) 0°74 (5) 1°25 (24) 
NGC 2440 0°60 (24) aes 1°25 (24) 
NGC 7027 0°49 Pe 0-72 (4) 0°43 (24) 

0°50 (24) 0°67 (24) 
a re ce (7) cnet (23) 1°33 (24) 
Orion 06 (24) 0°75 (24) 3°3 (24) 
IC 418 0-4 (24) 0°83 (24) an 
NGC 40 0°84 (7) aes 
NGC 2392 0°46 (7) 
NGC 6210 0°47 (7) oe 
BD + 30 3639 0°46 (7) 0°70 (23) 
IC 4593 0°67 (7) am 
NGC 6543 0°50 (7) 100 (23) 
NGC 6826 0°48 (7) nbs 
IC 4997 0°38 (7) oe ie 
NGC 7009 050 (7) sie 1°00 (24) 
NGC 6803 om re 0°87 (1) 
IC 2165 — ae ‘= (24) 

110 (x) 

IC 5217 sas 1°14 (24) 
Humason Planetary 0°42* ee 
Cannor Planetary 0°54* 


* Private communication (Osterbrock, 1956)+. 


2. Transition probabilities—The transitions under consideration are not 
allowed according to the selection rules for the electric dipole radiation but are 
allowed for the electric quadrupole and the magnetic dipole radiations. The 
theory of the transition probabilities is well known (Condon and Shortley (10) ; 
Shortley (21)). Earlier calculations of the transition probabilities were performed 
by Pasternack (15) and by Shortley, Aller, Baker and Menzel (22), hereafter 
abbreviated as SABM. The effects of the mutual magnetic interactions are very 
important for atoms having p* as their ground configuration. In general, the 
effect of the inclusion of the M—M interactions on the transformation coefficients 
is twofold. First, £, My and F, (the spin-orbit, the mutual magnetic and the 
electrostatic integrals respectively) are calculated empirically by equating the 


: oa of us (S.P.T.) is thankful to Professor Osterbrock for letting him know his experimental 
results. 
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theoretical expressions for the energy levels with their observed values, and 
therefore the values of ¢, M, and F, depend upon whether the mutual magnetic 
interaction terms are included in the energy level expressions or not. ‘The 
integral » used by AUV is five times our integral My. Second, the expressions 
for the transformation ‘coefficients (13) contain the M, terms explicitly. The 
empirical determination of ¢, M, for the p* configuration is further complicated, 
due to the fact that the diagonal elements of the spin-orbit interaction for this 
configuration are zero. ‘Therefore { and M, must be determined by interpolation 
between the neighbouring configurations p? and p*. 

For the calculation of the transition probabilities of O11 listed in Table II, 
we have employed the following values of the various integrals, calculated 
according to the above-mentioned procedure: 

€=181 cm". 
M,=0°35 cm". 

F, = 2979 cm“!. 

5, = 0°59ea,,” (Garstang’s value based upon wave functions with exchange). 
For comparison, we have included in Table II the transition probabilities 
calculated by AUV and Garstang. ‘The differences between their values and 
ours arise due to the different choices of ¢, My and s,._ AUV have used values 
of ¢ obtained from the observed energy levels but ignoring the mutual magnetic 
interactions. 

There is a further difference between our procedure and that used by AUV. 
In the formulae for the transformation coefficients, the denominators consist of 
the unperturbed term differences as given by the Slater formulae in terms of the 
electrostatic interaction integral F,. AUV have instead used the observed term 
differences which are affected by the configuration interaction. We feel that 
AUV’s procedure is not justified, and it may not correctly take into account 
the contribution of the configuration interaction to the transformation coefficients. 
In our procedure, the effect of the configuration interaction is taken into account 
at one step only, namely in the determination of F,, but not in the subsequent 
determination of the transition probabilities. 


TaBLe II 
Transition probabilities (in sec™) 
Transition Ou Si A Iv 
3°78 x 10-5 2°59 x 10-* 
D5)2-S3;2 5°41 x1075 (7) 
3°94 X 1075 (x1) 6:31 x 10-4 (11) 
16°51 X 107° 27°02 X 1078 
D3)2-S3)2 < I4°OI X 10 . (7) 
| 13°06 x 10~* (x1) 17°36 x 10-4 (11) 
P3)2-S3)2 00466 o'2711 
Pi)2-Ss2 0'0179 00663 
Ps, 2-Ds5,2 0°1643 0°4262 
12-D 5,2 0°0780 0°1087 
P3)2-Ds3;2 00478 0°1440 
P,,2-D3;3 0°1255 0°4830 


Transition probabilities for A1tv.—The ground configuration of Alv is 39%. 
For atoms having this configuration the effects of the mutual magnetic interactions 








466 Ali M. Naqvi and S. P. Talwar, The effect Vol. 117 


are not as important as for O 11 (ground configuration 2p*). Earlier, we mentioned 
that in general the effect of the inclusion of the mutual magnetic interactions is 
twofold. For atoms having 3p configuration, the first effect alone is significant, 
namely the effect on the calculation of and F,. The M, terms in the expressions 
for the transformation coefficients have negligible magnitude. Hence for atoms 
having 3p* configuration all we need to know is the accurate values of ¢ and F, 
but values of M, are not required. The values of ¢ and F, given by Naqvi (13) 
are: (=1256cm™; F,=2353 cm. 

The wave functions for Atv are not known. Therefore, the integral 5, has 
been determined by the screening constant method, similar to that used by 
Pasternack (15). This method is approximate, and for Alv we find 

Sy=0°974ea,”. 

The transition probabilities are given in ‘Table II. 

Transition probabilities for S11.—The transition probabilities for the 7D-‘S 
transition in S11 have been calculated by Garstang (11). For calculating the 
probabilities of other transitions we have used the following values of { and F, (13) 
and s, (1x): 

¢=497 cm™ 
F,= 1653 cm! 
$,=2'10ea,”. 

The transition probabilities for Sm are also included in Table II. 

3. Collision strengths.—Recently Seaton (16,17) calculated the collision 
strengths (the dimensionless parameter related to cross-section) for excitation 
by electron collisions between different terms of the ground configuration of O11 
and Sir. In our calculations in Section 5 we have needed to know these collision 
strengths (henceforth called as the collision cross-section parameters) from a 
level of one term to another term, and we have assumed them to be proportional 
to the statistical weights of the levels concerned. That is to say, if A and B are 
two terms and A consists of two levels, A; and A,, then the various collision 
cross-section parameters are given by, 


Q(4;,B) Q(A;,B) _ Q(A,B) 


WwW 





(1) 


Further, the collision cross-section parameters for excitation and de-excitation 
between any two states are equal, i.e. 


Q4p=Qp,. 


The collision cross-section parameters between the levels of the same doublet 
term are not known for the atoms considered here. The frequency of collisional 
processes between them is hence ignored. Since we completed this work, we 
have seen an estimation made by Seaton (18) of the cross-section parameters 
between the two levels of the 7D term of O11 and Su, but his procedure and the 
values do not appear to us to be justified. Until proper calculations for the 
cross-section parameters between the ?D levels of p? ions are made, it is considered 
safer to neglect the collisional processes between the doublet levels although 
considerable differences in the electron densities for regions of maximum 
emission of O11, Si and Av (as given in Table III) may result if different 
assumptions were made regarding the cross-section parameters for the 2D levels 
in p* ions. 
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We have used in Section 5 the cross-section parameters calculated according 
to equation (1); for Alv the cross-sections are estimated following a suggestion 
made by Seaton (16). 

4. Departure from thermodynamic equilibrium.—The p* configuration consists 
of the following levels: 4S, *Ds/o, 7D5/o, 2P 1/2, 2P 32. In Om the doublet terms 
are normal, but in Sit and Atv they are inverted. For convenience we shall 
adopt the following abbreviations. 


Term Abbreviations Level Abbreviations 
‘S S "Sse s 
2D. d’ 
2 ) ove 
I D 2D. d 
2p P “E 1/2 p’ 
L *Px p 


The intensity ratio of the D-S lines can be written as 
I(d'—s) N, A(d'—s) he(d'—s) 
"= I(d—s) N,  A(d—s) ‘ h(d—s) (2) 
where v is the frequency of the line A, the transition probability, and N,, N,- 
are the equilibrium populations of the levels concerned. In thermodynamic 
equilibrium, 





_ A(d’ —s) ho(d' —s) 
"th "3 Ad—s) h(d—s) ° 

Substituting the transition probabilities given in Table II, and the frequencies 
of the lines involved, the following values of ry, are obtained 

ry (Ol) =0°34, 
rp, (5 U1) =0°546, 
rq, (AIV) =0°14. 

A comparison with Table I will show that there is no thermodynamic 
equilibrium in gaseous nebulae. 

Bowen (g) showed that the atoms giving rise to the forbidden lines in gaseous 
nebulae are excited due to electron collisions. Radiative excitations are negligible 
on account of the extreme dilution of the radiation from the central star. The 
fact that stellar radiation is extremely diluted forms an essential reason for 
departure from thermodynamic equilibrium in gaseous nebulae. Proton impact 
could also produce transitions between the metastable ground configuration 
terms of ions with outer p? shells (¢=2, 3, 4) but recently Seaton (20) has shown 
that under the conditions prevailing in gaseous nebulae the number of such 
transitions is small compared to that produced by electron impact. 

The radiative and collisional de-excitations compete with each other. When 
the frequency of radiative de-excitations is small compared to that of collisional 
de-excitations, the collisional processes alone are responsible for excitation as 
well as de-excitation and virtually detailed balancing and thermodynamic 
equilibrium is obtained. On the other hand, when the number of radiative 
de-excitations is significant, we cannot expect to find thermodynamic equilibrium. 
Therefore the state of thermodynamic equilibrium can be expected only when 

(1) N, is large; 

(2) T, is small; and 
(3) transition probabilities are small. 
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The large departures of the observed ratio r for Atv from the Boltzmann 
value, as compared to that for O11, are due to the larger transition probabilities 
for AIv. 

5. Relative populations of the levels—The equations of equilibrium have 
been discussed earlier by Hebb and Menzel (12) and more recently by Seaton 
(17, 18) and Aller (2). Ignoring the frequency of transitions between the levels 
of the doublet terms (i.e. ?P,.2°?P,. and ?Dg. 2 *D5.) and neglecting the 
de-excitations from ?P term in the equations of the steady state for d’ and d levels, 
we can write for the ratio of the level populations 








Ny ____ CN,(Qsp+ Qppe-Xor/*e) + wpA(d—s) 

N, ~'5 CN,(Qgp + Qype %0r*™) +a,A(d —s) (3) 
Ny _ og CNA Qep + Qps) + wpA(p —s) + A(p—d) + A(p—d’') 
N,, ~°5 CN{Qpp + Qps) +opA(p’—8) + Alp’ —d) + A(p’—d’) 4) 


where C=8-54 x 10-*/T,1*. These ratios depend upon N, and T,. ‘The curves 
showing the variation of N,/N, with N,/T,)? for the ions under consideration 
are given in Fig. 1. For large values of N,/T,»?, all the curves in Fig. 1 
asymptotically approach the Boltzmann values. 

The variation of N,,/N,, is very small and hence is not included in Fig. 1. 
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Fic. 1.—Theoretical population ratio of *"Ds5;. and *D3;_ as a function of Ne/T¢!". 


6. Intensity ratio I(d'—s)/I(d—s) and the estimation of N,/T,?.—From 
equations (2) and (3), we obtain, 








CN, I 
1+ —(OQgp + Qppe*ne/*7e) Fans | a 
r=I(d'—s)|Id—s)=15 | “p Men 9 5)" 
seine —tpp/tT : 
[: + wp (Qep + Qype Xpp/ Aa as) = | 





he The exponential term in parenthesis is much smaller than the other term. ‘Therefore, 
while calculating Ne/T¢'/? (Table III, cols. 2, 3 and 4) T, in the exponential term only is taken 
equal to 10%. 
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The values of r, computed for various values of N,/T,"2 for Ou, Su and AIv 
are plotted in Fig. 2. All three curves asymptotically approach their thermo- 
dynamic equilibrium values r_),, (0-34 for O11, 0°546 for S11, and 0-14 for Atv). 
The curve for O11 reaches near the asymptotic value at N,/7'¥2~ 103, that for 
Sir at about 1o* and that for Atv at about 10°. 
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Fic. 2.—Theoretical intensity ratio r as a function of Ne T2'*. 


It would be worthwhile to utilize the theoretical curves to estimate the ratio 
N,/T,+? in the nebulae for which the intensity ratios are known observationally. 
The electron temperature 7, is usually computed from the relative intensities 
of the 44363 and the green nebular lines of [O11]. The electron densities N, 
are then estimated using the surface brightness method of Aller and Liller (3). 
However, on account of stratification which exists in most of the nebulae, real 
differences in electron density do occur within the body of the nebula. The 
surface brightness method should give essentially an overall average value of 
the electron density throughout the volume of the nebula, which may be much 
smaller than that in a cloud of maximum emission from an ion. 

More recently Seaton (17) obtained estimates of the electron temperature 7, 
and electron density N,, from the relative intensities of the forbidden lines. 
For ions like Nu, Ou, Orn, and Su, the intensity ratio of the nebular to the 
auroral or transauroral transitions gives a relation between T, and N,. Thus 
simultaneous transcendental equations involving N, and 7, are formed for 
each ion. They are then solved simultaneously to give N, and T,. However, 
with this method, one can expect a unique solution for N, and T, only when the 
ions considered arise in the same strata. Recently Aller and Minkowski (6) 
found that in NGC 7027 the simultaneous solution of N,, 7, equations for 
[Oi] and [Nu] lines gave (after applying corrections for space absorption) 
N.~6°3 x 104 electrons per cm® and T,~14000 deg. K. On the other hand 
employing only the [O11] data together with the photoelectrically measured 
surface brightness corrected for space absorption and the distance derived 

32 
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by Berman (8), Aller (2) found N,~1-7x10* electrons per cm* and 
T,~15600deg.K. These results may imply that [N 11] and [O11] radiations 
originate predominantly in different regions within the nebula, and that 
the forbidden line method (17) gives essentially an average density for the regions 
of maximum emission of the ions considered. 

The electron densities in the clouds of maximum emission of an ion in a 
nebula are not known. We can, however, make use of the theoretical curves 
in Fig. 2, to estimate the values of N,/7,2? in a nebula for the regions of maximum 
concentration of the ions considered here and thus get an idea of how this varies 
within a nebula for clouds of maximum emission of different ions. 

The various atomic constants used and the experimental intensity ratio r 
for O11 are quite accurate and hence the curve for O11 is quite reliable. The 
reliability of the theoretical curve for S11 is somewhat less, and that for Atv 
still less. It is difficult to give an exact measure of the reliability of these curves, 
but for Alv it may give no better than an order of magnitude estimate. This 
determination is most accurate for the steep portion of the curves. 

In Table III, we give the estimates of N,/T,!? for Om, Su and Atv for 
various nebulae. The values for O11 are perhaps quite accurate, whereas for Av 
they are rough estimates only. But still it appears that [S 11] and [A 1v] radiations 
are produced in filaments of much higher density than those from which [O11] 
radiation is more predominant, assuming that the variations in 7, are small. 
Hence a nebula, perhaps, consists of numerous filaments, knots, and tenuous 
regions such that the density may very well range from less than 104 electrons/cm® 
to something like 10° electrons/cm*. 





TaBLe III 
Values of N,/T.}" for O11, S11 and Atv 
Nebulae N,/T,}? NJT3* N,/T,}? NASB) NAFL) 
Ol Su AIv T 2 “Tt 
NGC 7662 80 550 760 51 440 
NGC 2440 92 aa 760 _ 310 
NGC 7027 174 562 1°2 X 10 38 290 
160 1000 ' 29 
NGC 6572 <n { 540 480 88 440 
Orion 92 540 re ala 180 
IC 418 500 340 dec 117 60 
NGC 40 38 Pia os ee 
NGC 2392 230 a as 8 
NGC 6210 230 ae Se 84 
BD +-30 3639 230 740 a iat 
IC 4593 66 we = 33 eae 
NGC 6543 176 160 ue 65 300* 
NGC 6826 190 a ose 21 286* 
IC 4997 750 en aes ap co 
NGC 7009 160 oes 1°99 X103 53 253* 
NGC 6803 re = 2°82 x 10 Sete 
19 X103 
IC 2165 ee ae 4 apa" 
IC 5217 me ve 1*15 x 108 
Humason Planetary 400 Bee 
Cannor Planetary 125 


*Estimated from surface brightness (Seaton, (19)). 
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For the sake of comparison, we have included in Table III, the values of 
N,/T,** computed from the surface brightness method of Aller and Liller (3) 
(abbreviated as SB) and the values obtained by Seaton (17) using the forbidden 
line method (abbreviated as FL). 

Our values of N,/T,’ for regions of maximum [O 11] emission are, in all but 
two cases, about 2 to 5 times greater than Aller’s values. The values of N,/7T,1 
given by the FL method (17) are all larger (with the exception of IC 418, for which 
Wyse’s intensities may not be reliable) than our values for Ou, which can be 
explained if [O11] radiation predominantly originates from less dense portions 
of the nebulae, since the values given by the FL method are essentially average 
values for the regions of maximum emission of the ions considered. 

Very recently Osterbrock (14) observed the 43727 intensity ratio at many 
points in the Orion Nebula. The intensity ratio is observed to vary much at 
different points, thus supporting the theory. All the observed values lie between 
the lower limit 0-34 and the upper limit 1-5. 

It is admitted without reservation that this method of estimating the electron 
density in clouds of maximum emission of p* ions can be effectively applied only 
if the intensity ratio r of the *D-—*S lines is known accurately. Further, this 
determination should be particularly accurate for the steep portion of the curves. 
It is here assumed that there are no large scale variations in the value of electron 
temperature 7’, within a nebula. But, strictly speaking, one should determine 7, 
for the ion concerned from the ratio of the nebular and the auroral or transauroral 
lines. Thus accurate intensity measurements of [D—S], [P-D], and [P-S] lines 
for Ou, Si and Atv in gaseous nebulae are much to be desired. 

We are grateful to Professor D. S. Kothari for helpful discussions. 


Department of Physics, 
University of Delhi: 
1955 April. 
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